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‘jeder dieser Zustände höchstens einmal besetzt sein. 


ANNALEN DER PHYSIK 


Prozesse zwischen leichten Teilchen nach der Diracschen Theorie*) 
Von B. Kockel 


(Mit 1 Abbildung) 


Inhaltsiibersicht 


Solche zwischen Elektronen und untersucht, in 
denen die Zustandsänderungen der Elektronen ohne Beteiligung schwerer Teil- 
chen unter Emission oder Absorption von Lichtquanten geschehen. Die Wahr- 
scheinlichkeit für das Eintreten der einfachsten Prozesse dieser Art (Compton- 
effekt, Erzeugung oder Zerstrahlung eines Elektron-Positron-Paars) sind in Ab- 
schnitt 5 noch einmal in einheitlicher Rechnung bestimmt. Abschnitt 6 gibt die 
Wahrscheinlichkeit der Zerstrahlung eines Paars unter Bildung von drei Licht- 
quanten. Dabei ergibt sich diese ,,Dreifachzerstrahlung,, bei hoher Energie des 
Paars mit derselben oder größerer Wahrscheinlichkeit als die ,,Zweifachzerstrah- 
ling‘, was natürlich nur bedeutet, daß bei solchen Energien die Einteilung der 
Zerstrahlungsprozesse nach der Zahl der erzeugten Lichtquanten unbrauchbar 
ist. Abschnitt 3 zeigt, daß keine Prozesse möglich sind, in denen bei Beteiligung 
beliebig vieler Lichtquanten mehr als ein Elektron seinen Zustand wechselt, z. B. 
ist also die simultane Bildung zweier Paare aus Lichtquanten unmöglich. Die 
Tragweite dieses Ergebnisses und seine Abhängigkeit von den benutzten Eigen- 
funktionen der Elektronen und dem störungstheoretischen Ansatz werden auf 
Seite 281 und 282 diskutiert. 


f 
In der Diracschen relativistischen Theorie des Elektrons findet 1 man wie in. 
der nichtrelativistischen Quantenmechanik, daß den Elektronen bestimmte durch 
ihre Eigenfunktionen gekennzeichnete Zustände zur Verfügung stehen. Durch 
eine geeignete Periodizitätsbedingung kann man immer erreichen, daß die Zu- 
stände eine diskrete Mannigfaltigkeit bilden. Nach dem Pauli-Prinzip kann dann 
Die Zustände negativen 
Energievorzeichens sind nach Dirac fast alle besetzt; die volle Besetzung aller 
Zustände negativer Energie soll keine physikalische Bedeutung haben; unbe- 
setzte Zustände negativer Energie bedeuten Positronen. Jedes Elektron kann 
unter dem Einfluß einer Störung in einen anderen vorher nicht besetzten Zustand 
übergehen. Nach der Diracschen Deutung haben solche Übergänge ja nac :h dem 
Energievorzeichen im Anfangs- und Endzustand des Übergangs eine andere 
physikalische Bedeutung: ++ Elektronensprung, — — Positronensprung, — + Er- 
zeugung eines Elektron-Positron-Paars, + — Vernichtung eines solchen Paares. 


*) Eingereicht zur Erlangung u venia legendi in der Philosophischen Fakultät 
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Die den Übergang auslösende Störung kann verschiedener Art sein. In dieser 
Arbeit sollen solche Prozesse untersucht werden, in denen der Elektronenüber- 
gang durch ein Strahlungsfeld, also durch Absorption oder Emission von Licht- 
quanten geschieht. An solchen Prozessen können schwere Teilchen beteiligt sein, 
sei es in der Form, daß alle oder ein Teil der beteiligten Elektronenzustände an 
einen Kern gebundene Zustände sind (Beispiele: Absorption, Emission, photo- 
elektrischer Effekt), sei es, daß die Mitwirkung eines Kernes notwendig ist für den 
Ablauf des Prozesses (Beispiel: Bremsstrahlung). Weiter gibt es Prozesse, die 
ohne jede Beteiligung eines schweren Teilchens ablaufen; Prozesse dieser Art 
sollen in dieser Arbeit näher untersucht werden. Als einfachste derartige Prozesse 
ergeben sich die vier Prozesse der folgenden Aufstellung, alle dadurch gekenn- 
zeichnet, daß ein Elektron unter dem Einfluß zweier Lichtquanten in einen an- 
deren Zustand gelangt. 

Elektron und Lichtquant — Elektron und Lichtquant (Comptoneffekt) 

Positron und Lichtquant — Positron und Lichtquant (Comptoneffekt am 
Positron). 

Zwei Lichtquanten > Elektron und Positron (Paarerzeugung, Materialisation 
von Lichtquanten). 

Elektron und Positron — zwei Lichtquanten (Paarvernichtung, Zerstrahlung 
eines Elektron-Positron-Paars). 

Die zunächst auch denkbaren noch einfacheren Prozesse, bei denen der Elek- 
tronenübergang unter Beteiligung nur eines Lichtquants geschieht — das wären 
Prozesse, die aus der obigen Aufstellung hervorgehen, wenn man ein Lichtquant 
streicht — treten nicht auf, weil sie nur unter Verstoß gegen Impuls- oder Energie- 
satz möglich wären!). 


2. Prozesse höherer Ordnung 


Kompliziertere als diese einfachsten Prozesse — der Einfachheit halber sämt- 
lich als ,,Prozesse höherer Ordnung“ bezeichnet — sind in zweierlei Weise denkbar. 
Zunächst können nicht ein, sondern zwei oder mehr Elektronen unter Absorption 
und Emission von Lichtquanten in einen anderen Zustand übergehen. Das wären 
also Prozesse, bei denen Anfangs- und Endzustand zusammen mehr als zwei 
materielle Teilchen enthalten. Z. B. gehört hierzu die Erzeugung von zwei oder 
mehr Elektron-Positron-Paaren in einem Strahlungsfeld. Weiter können die 
Prozesse der obigen Aufstellung unter Beteiligung von mehr als zwei Lichtquanten 
auftreten. Ein Beispiel ist der Comptoneffelkt unter Bildung von zwei oder mehr 
Streuquanten. Die Prozesse heißen dann Comptoneffekt, Paarerzeugung oder 
Paarvernichtung ,,in höherer Ordnung“ ?). Schließlich können beide Erweiterun- 
gen der einfachsten Prozesse zugleich auftreten). 


1) Im Feld eines Kerns treten sie auf, weil der Kern überschüssigen Impuls auf- 
nehmen kann. 

2) Hier hat die Bezeichnung ,,héhere Ordnung‘‘ einen wohldefinierten Sinn, denn die 
Diracsche Störungstheorie — siehe die Zitate bei (1) — muß bei Prozessen, in denen 
n Lichtquanten auftreten, in n-ter Ordnung angewendet werden, und die Wahrscheinlich- 
keit dieser Prozesse und die zugehörigen Wirkungsquerschnitte werden proportional e?*, 
wo e die Elektronenladung ist. 

3) Die obige Aufstellung und die besprochenen Erweiterungen unfassen noch nicht 
alle unter Emission und Absorption von Lichtquanten ablaufenden Prozesse. Beispiel: 
Streuung von Licht an Licht, wo im Anfangs- und Endzustand gar keine materiellen 
Teilchen vorhanden sind. H. Euler u. B. Kockel, Naturwiss. 28, 246 (1935); H. Euler, 
Ann. Physik (5) 26, 398 (1936). : 
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Solche Prozesse höherer Ordnung können in der Form ablaufen, daß in einem 
unterhalb der Erkennbarkeit kleinen zeitlichen und räumlichen Abstand zwei 
einfachere Prozesse auftreten. Z. B. könnten so vier Lichtquanten zwei Elektron- 
Positron-Paare bilden, indem jeweils zwei Quanten ein Paar bilden. Gemäß der 
Zusammensetzung aus zwei (oder mehr) einfachen Prozessen müssen schon für 
diese selbst Impuls- und Energiesatz erfüllt sein. Bei dem als Beispiel genannten 
Prozeß, dessen Gesamtbilanz durch: 4 Lichtquanten — 2 Paare gekennzeichnet 
ist, würden also zwei Impulsgleichungen und zwei Energiegleichungen gelten: 


+ Ge P-ı + C9, + pi, + + pin 


(g Lichtquanten-, p_ Elektronen-, p.. Positronenimpulse). Darüber hinaus können 
aber auch Prozesse auftreten, bei denen Impuls- und Energiesatz nur für den ge- 
samten Prozeß erfüllt sind: ,,Eigentliche Prozesse höherer Ordnung‘. Bei der 
Erzeugung zweier Paare aus vier Lichtquanten wären also nur die beiden Glei- 


chungen 
+ Ge + Gs + Ga = + P-2 + + Pte 


=c + + c? + + Vm? c? + pi, + + Pia} 


erfüllt, nicht aber entsprechende Gleichungen, in denen nur einige der acht be- 
teiligten Teilchen auftreten. 


In den folgenden Abschnitten wird für die eigentlichen Prozesse höherer Ord- 
nung gezeigt: 


Dritter Abschnitt: Im Rahmen der benutzten Ansätze für die Eigenfunk- 


tionen der Elektronen und ihre Wechselwirkung mit elektromagnetischer Strah- = 

lung findet man: Es gibt keinen Prozeß höherer Ordnung, in dem unter Absorp- | 

tion und Emission von Lichtquanten mehr als ein Elektron seinen Zustand wechselt. : 
Es gibt also z. B. keine Erzeugung von mehr als einem Paar, keine Zerstrahlung 

von mehr als einem Paar unter Absorption oder Emission von beliebig vielen 
Liehtquanten. Ebenso gibt es nicht den entsprechenden Comptonprozeß höherer , 
Ordnung, also keinen Prozeß, bei dem n (> 2) Elektronen + absorbierte Licht- 

quanten in n andere Elektronen + emittierte Lichtquanten übergehen. pee. 


Man muß sich bei diesem, Prozesse der beschriebenen Art absolut ausschlieBen- 4 a 
den Resultat allerdings darüber im klaren sein, daß eine exaktere Theorie zu Ein- Er 
schränkungen dieses Ergebnisses führen könnte. Die Elektronengesamtheit ist 
in der später folgenden Rechnung unter Überlagern der individuellen Eigen- 
funktionen aller Elektronen beschrieben worden, statt, wie es richtig wäre, durch _ RS i 


eine Eigenfunktion in den unendlich, vielen Variablen aller Elektronen. Damit ist | 


z.B. die gegenseitige Beeinflussung der Elektronen — gleichen oder ungleichen 
Spins — nicht beriicksichtigt. Dariiber hinaus ist die Wechselwirkung der Elektro- has “a 
nen mit den Lic -htquanten nicht so beschrieben worden, daß die Eigenfunktionen _ 2. 4 
der Elektronen im Feld der Lichtquanten exakt berechnet wurden. Die W. echsel- 29S | 


wirkung ist vielmehr durch ein Störungsverfahren genähert bestimmt worden. 35 


kann also sein, daß eine exaktere Theorie doch eine nicht verschwindende Wahr- 3 
19* 
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iss für den gleichzeitigen Niveauwechsel von zwei oder mehr Elektronen 
ergibt. Da aber die in der Rechnung verwendeten genäherten Verfahren in vielen 
Fällen zu gut mit der Erfahrung übereinstimmenden Resultaten führen (z.B, 
Comptoneffekt, Abbremsung von Elektronen in Materie bei Energien bis etwa 
10° eV *)) und deshalb anzunehmen ist, daß eine exaktere Theorie nur zu gering- 
fügigen Modifikationen der bisherigen Regultate wenigstens im Bereich uicht zu 
hoher Energie führen wird, können wir erwarten, daß von unserem Ergebnis 
folgendes richtig bleiben wird: Der gleichzeitige Niveauwechsel von zwei oder 
mehr Elektronen unter dem Einfluß von Lichtquanten ist wenigstens für nicht 
zu hohe Energie wesentlich weniger wahrscheinlich, als nach der Anzahl der be- 
teiligten Lichtquanten zu erwarten wäre. In Formeln: Für die zugehörigen Wir- 
kungsquerschnitte ist bei n beteiligten Lichtquanten 


mc?) \he 

einer dimes isionslosen Kombination von mc und Impulsen, die von der Größen- 
ordnung 1 ist, wenn die Impulse die Größenordnung m c haben, zu erwarten. Statt 
dessen wird ® < eine solche Größe. 

Vierter Abschnitt: Dagegen gibt es alle solche Prozesse höherer Ordnung, in 
denen nur ein Elektron seinen Zustand wechselt unter Beteiligung von mehr als 
der Mindestzahl von zwei Lichtquanten, die bei den entsprechenden einfachen 
Prozessen auftreten. Folgende Prozesse sind also möglich: Comptoneffekt, ein 
Elektron + mindestens ein absorbiertes Lichtquant gehen über in das gestreute 
Elektron + mindestens ein emittiertes Lichtquant; Paarerzeugung, mindestens 
zwei Lichtquanten bilden ein Elektron-Positron-Paar und n >0 emittierte Licht- 
quanten; Zerstrahlung, ein Paar geht unter gleichzeitiger Absorption von n >0 
Lichtquanten über in mindestens zwei Lichtquanten. Dabei sind in dieser Auf- 
stellung die einfachsten Prozesse mit enthalten; man erhält sie, wenn nur die 
jeweils genannte Mindestzahl von Lichtquanten (0 oder 1 oder 2) beteiligt ist. 

Fünfter Abschnitt: In diesem Abschnitt werden Wahrscheinlichkeiten und 
Wirkungsquerschnitte für die — schon in der Einleitung genannten — einfachsten 
Prozesse noch einmal in einheitlicher Rechnung zugleich abgeleitet. 

Sechster Abschnitt: Die Wahrscheinlichkeit eines der möglichen Prozesse 
höherer Ordnung, die Zerstrahlung eines Paars unter Bildung von drei Licht- 
quanten wird berechnet und mit der des entsprechenden einfachen Prozesses ver- 
glichen. Dabei ergibt sich, daß diese „Dreifachzerstrahlung‘“ erst für sehr hohe 
Energie ebenso wahrscheinlich wird wie die „Zweifachzerstrahlung‘, so daß im 
experimentell erreichbaren Energiebereich sowohl die Lebenszeit langsamer 
Positronen wie der Anteil der Zerstrahlung an der Abbremsung schneller durch 
Materie fliegender Positronen weiter durch die Zweifachzerstrahlung bestimmt 
bleiben. In dem Energiebereich, in dem die Wahrscheinlichkeit der Dreifach- 
zerstrahlung ebenso groß wie die der Zweifachzerstrahlung wird, ist daselbe von 
der Vierfach-, Fünffachzerstrahlung usw. zu erwarten; das heißt soviel wie: In 
diesem Energiebereich sind das ganze Verfahren, die Prozesse nach der Anzahl 
der beteiligten Lichtquanten einzuteilen, und die dieser Einteilung angepaßten, 
in der Rechnung verwendeten störungstheoretischen Ansätze samt den abgelei- 
teten Resultaten als illusorisch anzusehen. x 


4) Vgl. W. Heitler, The Quantum Theory of Radiation, 1936. 
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3. Prozesse mit mehr als zwei materiellen Teilchen 


Die Ubergangswahrscheinlichkeit fiir Prozesse zwischen Elektronen, Posi- 
tronen und Lichtquanten ist 5) 


2n Via Vox V, 
W is "93 "36 n—In 1 
Y [chh 9 
k,k+ı = e) OV (X €) (2) 


Dabei bedeuten die Indizes 1, 2,3...n Gesamtzustände, also eine bestimmte 
Konfiguration von Lichtquanten und besetzten Elektronenzuständen. Die mög- 
lichen Elektronenzustände sind auf das Volumen V normierte ebene Wellen 


7 Qyoa, also gekennzeichnet durch die drei Angaben p (Impuls), o (Spin), A (Ener- 


gievorzeichen) und den zugehörigen Diracschen Vierervektor ay,,. Das Matrix- 
element Yz,x+1 der Störungsenergie beschreibt einen Übergang vom Zustand k 
zum Zustand & + 1, der dadurch geschieht, daß das in k vorhandene Elektron 
(p’, o’, A’) unter Aufnahme des Impulses g (Betrag g) eines absorbierten Licht- 
quants oder unter Abgabe des Impulses g eines emittierten Lichtquants übergeht 
in den Zustand (p’ +- g, o’’, 2’’). Die Gesamtzustände k und k + 1 unterscheiden 
sich also dadurch voneinander, daß ein Lichtquant neu entstanden oder ver- 
schwunden ist und ein Elektron unter Einhaltung des Impulssatzes einen Niveau- 
wechsel vorgenommen hat. e ist die Polarisationsrichtung des betreffenden Licht- 
quants, ve=n,e,+a,e,+%,e,, wo die x die Matrizen der Diracgleichung 


W 
—y= Py + P,) y sind. 


Die Nenner FE, — E, stellen die Summen aller vom Anfangszustand 1 bis zum 
Zwischenzustand k freigewordeneft Energien dar: 


E, —E, (E, — E,) +(E, — E;) (E,—E,) (Ex -ı —E,) 
wo : 
E,;_,—E;=ci' + Vm? c? + (p’ + g)? eg 


ist. (Oberes Vorzeichen bei Absorption, unteres: Vorzeichen bei Emission). Für 


den gesamten Prozeß gilt der Energie-Erhaltungssatz E, = E,,. _ bedeutet die 


Dichte der Endzustände auf ihrer Energieskala. Das Summenzeichen & in (1) 
bedeutet ‚Summe über alle Möglichkeiten, vom Zustand 1 unter Emission bzw. 
Absorption aller beteiligten Lichtquanten zum Zustand n zu kommen“. Wie 
man Anzahl und Art dieser Möglichkeiten bestimmt, ist mehrfach in der Literatur 
behandelt worden ®). 

Anfangs- und Endzustand mögen sich nun dadurch voneinander unterscheiden, 
daß zwei Elektronen in einen anderen Zustand gelangt sind. Das eine 
—— kann seinen Endzustand etwa über folgende Zwischenzustände erreicht 

aben. 


Dos Ap, To > Pr, > Po, Ap, 
5) Uber die Herleitung dieser aus dem Ansatz ebener Wellen fiir die Eigenfunktionen 
der materiellen Teilchen, der Quantisierung des Strahlungsfeldes und der Diracschen 
Störungstheorie hervorgehenden Formel vgl. z. B. W. Heitler, 1. c. — H. Euler, |. ce. 
—B. Kockel, Z. Physik 107, 153 (1937). 
“ hy Vel. die weiter unten zitierten Arbeiten und fiir einen komplizierteren Fall: H. 
uler, l.c. 
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Das andere mége folgende Zustandsw echsel durchgeführt haben 


Po, Ao, 06 Pi, Ai, 01 Po, Ag, 0g Pin: 


(4) 


(m kann =n sein). Bei einer Erzeugung zweier Paare wären A, = 2, = —1, 
An = /m= +1, bei einer Zerstrahlung zweier Paare wäre es umgekehrt. Suchen 
wir alle diejenigen Übergangsmöglichkeiten auf, die über die in (3) und (4) ausge- 
wählten Zwischenzustände führen, so finden wir jedenfalls, daß erst die Ubergangs- 
reihe (3) vollständig und danach die Reihe (4) ablaufen kann, weiter, daß zuerst 
ein Teil der Reihe (3), dann ein Teil der Reihe (4), dann der Rest der Reihe (3), 
dann der Rest der Reihe (4) ablaufen kann usw., kurz, die Reihen (3) und (4) 
können beliebig ineinander geschachtelt werden, z. B. in leicht verständlicher Schreib- 
weise 
(01) (12) (23) (34) (1’2’) (2’3°) (45) (3’4’) (56) (67).... 


.(m — 1’,m’) (n—2,n—1) (n—1,7n) 6) 


km 


Für alle auf diese Weise erfaßten Ubergangsméglichkeiten ist das Produkt der 
(kommutativen) Größen V x,x+1 in (1) gleich groß. Die Diskussion der & über alle 
Übergangsmöglichkeiten reduziert sich also auf die Erörterung der Summe über 
die Energiedifferenzprodukte im Nenner von (1). Wir bezeichnen die vom Anfangs- 
zustand bis zum k-ten Zustand freiwerdende Energie (einschließlich der Licht- 
quantenenenergie) für die Reihe (3) mit F,, für die Reihe (4) mit G,, z. B. also 


PF, = ch, Vm? ct + + pi + e|pı 
| ps —pil. 


Auf Grund des Energiesatzes ist F,, + G,, = 0. Jeder der Nenner in (1) läßt sich 
in der Form er 


v=) v=k, 1 v=k m-1 


Dabei sind 0 <k, Sk, Sk,...k, <n Zahlen zwischen 0 und n. Das erste 
und das letzte Produkt sind durch 1 zu ersetzen, wenn der obere Index kleiner 
als der untere ist. Jede Zahlenfolge k,...k,, stellt eine Möglichkeit der Inein- 
anderschachtelung der Übergangsreihen (3) und (4) dar; z. B. bedeutet k, = k, = 

. ky, = 0, daß erst die Reihe (4) vollständig, dann die Reihe (3) abläuft, umge- 
kehrt bedeutet k, = k, = ...k,, =, daß erst die Reihe (3) vollständig, dann die 
Reihe (4) abläuft. Die Zahl der Möglichkeiten ist gleich der Zahl der Kombina- 


tionen von n Elementen zur m-ten Klasse mit Wiederholung = E <a . Fiir die 


Berechnung von (1) ist der Ausdruck 


ky ks n 


alle k-Folgen 
auszuwerten. Es ist 


IT F, G, II’ (F,+6,)° 


?) IT’ heißt außer dem Fall, daß zugleich v = n, u = m ist, also außer F,, + @„- 


also 


abhi 


ode 


| | 
\ 
Ny, 
{7 — 
= 
nh 
1 
= 
nac 
unc 
sch 
aus 
un 
in 
3 abs 
na 


IT (Fe +@,) IT (Fy +6,) 
(4) 0 ke+1 
(F, +G,,-1) II (F,+@,); 
zuerst kh=0h= 
he (3), Wir werten zuerst die Summen über k, aus (unter Weglassen der von &, nicht 
reib- kr n ke 
= IF IT (F, +6)= „I, F, i=, ii +6,) 
(5) = pF, li (F, +4G,). 
td tae 
Entsprechend erhalten wir für die Summe über k, 
er alle 
IT (F, +@) (Fy +@) IT Fr. 
Im weiteren ie der Rechnung entstehen die RT +@)), 
t sich 4 
II(F, +@,) usw. Der letzte Schritt der heißt 
1 
G,,,). 
1 
erste n k 
leiner » +G,-1) > II F, ire (F, + G.,,) 
'nein- 1 Okm+1 tet, TALS 
n n 
imge- A=— II (F, (F,+@,).. TI (F, [] (F, +G,,) 
bins- ater schlieBlich 6, 
Fy +6 
A= “u=1bism F, + Gn 
Bir, I Fee TT Ge (Pe + Gy) Fr. II Gu 
1bisn 1bism 1bisn bisn 1bism = 
(7) a=1bism 
nach dem Energiesatz. Für jede Kombination zweier solcher Übergangsreihen (3) 
und (4) ist also der Beitrag zur Summe in der Formel für die Übergangswahr- 
scheinlichkeit (1) gleich Null. Da die Gesamtheit aller Übergangsmöglichkeiten 
aus solchen Kombinationen zusammengesetzt ist, ist W=0, wenn Anfangs- 
und Endzustand sich dadurch voneinander unterscheiden, daß zwei Elektronen 
in einen anderen Zustand gelangt und beliebig viele Lichtquanten emittiert oder 
absorbiert worden sind. Ein solcher Prozeß tritt also nicht ein, auch wenn er 
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Die Rechnung der letzten beiden Seiten läßt sich erweitern für den Fall, daß 
N >2Elektronen in einen anderen Zustand gelangen. Man muß dann N 
Übergangsreihen (3) kombinieren, d. h. nach Art von (5) ineinander verschachteln, 
Für ihren Beitrag zur Übergangswahrscheinlichkeit hat man einen (7) entspre- 
chenden Ausdruck auszuwerten, für den man 


3 


erhält, weil der Zähler die gesamte freiwerdende Energie darstellt, die nach dem 
Energiesatz gleich Null sein muß. Also gilt folgender Satz allgemein: Die Über- 
gangswahrscheinlichkeit ist gleich Null für solche Prozesse höherer Ordnung 
zwischen Elektronen und Lichtquanten, bei denen sich Anfangs- und Endzustand 
dadurch voneinander unterscheiden, daß zwei oder mehr Elektronen in einen ande-, 
ren Zustand gelangt sind ®) ®). i 


er; 4. Comptoneffekt, Paarerzeugung und Zerstrahlung bei Beteiligung von 
4 . beliebig vielen Liehtquanten 


Der vorige Abschnitt sagt nichts aus über Prozesse, in denen nur ein Elektron 
in einen anderen Zustand gelangt. Solche Prozesse bleiben also möglich. Nach 
dem Energievorzeichen des einen Elektrons, das bei diesen Prozessen von einem 


Anfangsniveau Pp, do, A, nach dem Endniveau p,, o,, A, gelangt, hat man . 
A= +1 A,= +1 Comptoneffekt am Elektron, 
44 =—1 A, =—1 Comptoneffekt am Positron, 


Ag = +1 A,=—1 Zerstrahlung, 
Ag =—1 A, = +1 Paarerzeugung. 


Impuls- und Energiesatz verlangen, wenn g,, g, die Impulse absorbierter und 
emittierter Lichtquanten bedeuten 


Sind A, und A, gleich (Comptoneffekt), so kann man ohne Einschränkung 
der Allgemeinheit in ein Bezugssystem übergehen, in dem 9, = — py ist. Das 
neue System muß dabei gegen das ursprüngliche Bezugssystem die gleichförmige 
Geschwindigkeit ¢ — 
m2 + pz + a+ pe 


am Je 


- haben. Im neuen System gilt dann 


Das ist sicher méglich, wenn beide Summen mindestens ein Lichtquant ent- 
halten, also muß mindestens ein ‚Lichtquant absorbiert und mindestens eins 
emittiert werden. 

Sind), und} ungleich, so kann man ohne Einschränkung der Allgemeinheit 
in ein Bezugssystem übergehen, in dem p} = pj ist. Das neue System hat dabei 
Po — Pn 


gegen das ursprüngliche die Geschwindigkeit ce - - 
8 Vm pi + Vm? + pe 


8) An diesem Ergebnis ändert sich auch nichts, wenn man, um der fast vollständigen 
Besetzung der Zustände negativer Energie Rechnung zu tragen, die Jordan-Wigner- 
schen Vorzeichenfunktionen einführt. 

9) Also können z. B. nicht drei Lichtquanten zwei Paare bilden, was nach Energie- 
und Impulssatz durchaus möglich wäre. B. Kockel, Z. Physik 107, 153 (1937). 
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Im neuen System gilt ikea 
= Ga = 2% (8) 
\ 
(9) 


Ist A, —Ap = + 2 (Paarerzeugung), so müssen mindestens zwei absorbierte Licht- 
n 5 

quanten vorhanden sein; denn wenn nur eins vorhanden wire, so wiirden aus (8) 

und (9) die einander widersprechenden Beziehungen 


“sake 


folgen. Ist umgekehrt A, — A, = — 2 (Zerstrahlung), so müssen mindestens zwei 


Lichtquanten emittiert werden; würde nur eins emittiert, so wäre zugleich 

was nicht möglich tst. Im folgenden fünften Abschnitt werden die einfachsten 
möglichen Prozesse, in denen also gerade zwei Lichtquanten beteiligt sind, in 
einheitlicher Rechung behandelt. Sie sind von „zweiter Ordnung“, weil infolge 
des Auftretens zweier Lichtquanten in (1) das Produkt zweier Matrixelemente 
V,., .+1 auftritt. Im sechsten Abschnitt werden dann die Prozesse „dritter Ord- 
nung“, also die der Prozesse, an denen drei Lichtquanten beteiligt sind, unter- 
sucht. 


5. Comptoneffekt, Paarerzeugung und Zerstrahlung in einheitlicher Behandlung 

Wir können die Berechnung der Wahrscheinlichkeit dieser Prozesse einheitlich 
durchführen, wenn wir gemäß dem Schema von 8.286 Energievorzeichen A der 
Elektronenzustände und für Impulse und Energien emittierter Lichtquanten eine 
Vorzeichenumkehr einführen. g und g (eine dann negative Größe!!) bedeuten für 
ein emittiertes Lichtquant, daß es den Impuls —g und die Energie —cg besitzt. 
Danach ergibt sich folgendes Schema für alle diese einfachen Prozesse. 


Anfangszustand des Elektrons: p,, A», 0% Energie = Agc & 
Zwischenzustand entweder: p’=p,+ g9',4,0’ Energie = 2’ ce’ 
oder: p’=p,+9",4",0"  Energie= A’ ce” 
Endzustand : p=p. + 9",4,0 Energie=Ace 


2 


q’, 9’ seien die beiden beteiligten Lichtquanten, e’, e’’ ihre Polarisation. Die RE 


Energiedifferenz zwischen Anfangs- und Zwischenzustand ist entweder 
+g’—HeE) oder ce (Aue A” &’’), wo die positive Vm: ce + po 
usw. ist. Nach den Gleichungen (1) und (2) erhalten wir 


2a h? [a¥ (x e’) a’] [a’* (x e’ Jal 
222 Che tg 


RAE | 
Die Summe in (1) über alle Übergangsmöglichkeiten wird zu der Summe über 
die beiden Energievorzeichen und Spinzustände der Zwischenzustände und zur 
Summe über die beiden Fälle ,,entweder und ‚oder‘ der obigen Aufstellung, was 
ausgedrückt werden kann durch = Summe iiber die beiden durch Permutation 


W = 


der Lichtquanten entstehenden Ausdriicke. Im allgemeinen interessiert man sich 
nicht dafiir, in welchem Spin- oder Polarisationszustand die zu Anfang und Ende 
des Prozesses vorhandenen Teilchen sich befinden. Man fragt statt dessen, mit 
welcher mittleren Wahrscheinlichkeit der untersuchte Prozeß geschieht, wenn in 
dem Strom der einfallenden Teilchen keine Spin- oder Polarisationsrichtung be- 
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vorzugt ist. Um diese vereinfachte Frage zu beantworten, hat man in (10) über 
die Spins oder Polarisationen der vor dem Prozeß vorhandenen Teilchen zu mitteln 
und über die Spins und Polarisationen der nach dem Prozeß vorhandenen Teilchen 


zu summieren. Das ergibt 
2? Qn ei ec? h? h? 1 [at (a e’) [a’* (x e”’) a] |? 
W = RAE 122 222 C (Ape +9’ —A’e’) |” (11) 


Sie Dummation iiber ‘ie o’,2’, o’, 2’ kann man in bekannter Weise durch- 
führen!) unter Benutzung der Operatoren h = ap +B mc usw. Man erhält 


49 — (Po 9’) &0 9° — (Po 9’) 


Summierung über die Spinzustände o,, o nach dem gleichen Verfahren ergibt 


eth? h 1 
od 32 4E V?|g’g” ay 4 3 
24 


(Age +9’ + h’) (a e” h) 
5 sur > Jo 9’ — (Po 9’) (de 


+ h’) (xe’) 
h 
Jo 9’ — (Po 9°) (Ao + a}. h 


Die weitere Rechnung wird weitgehend vereinfacht, wenn mas ein Bezugssystem 
benutzt, in dem p, = 0 und damit ¢, = mc ist. Man erhält dann 


n ehh 


‘ 

a , [77 
Ir (xe’) (Agme +g’ +h’) 
€ Spur | = (Ae +h) 


g 


Wegen der beiden zweigliedrigen Summen über die Permutation besteht die Spur 
aus vier Teilen. Man erhält unter eval Summierung iiber die Polari- 


= & Spur...= 32 (me +49 — 
— 16 cthg’ + Ag”) 
bar 32 (m c T ho g’ Ao — ’ sin? 


10) Vgl. sämtliche später genannten Arbeiten. _ 
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, über wo # der Winkel zwischen g’ und g’’ ist. Für W ergibt sich en 


itteln wa 2 eh? h 1 
Ichen AE Ve|g’g’ | 


{2 me—mcesin?# + Ag" 
(12) 7 
1 
Unter Berücksichtigung der aus dem Impulssatz p = g’ + g’’ und dem Energie- 
satz Ae=A,mc +g’ +g” folgenden Gleichung 


It 9") =9' 9" +Aome(g’ +9”) 
kann man (12) auch umschreiben in 
an 1 1 1 1 
(12”) 
Die weitere Berechnung spaltet sich von hier aus für die vier verschiedenen 
* Prozesse auf. 
Comptoneffekt: Hier ist A = A, = + 1 (oder beide = — 1 für den Compton- 
effekt am Positron, was an den Formeln nichts ändert). g’’ ist das emittierte Streu- 
quant des Effekts, also führen wir von jetzt an wieder g’’ in — 9”, g’’ in —g’”’ über. 


t W => AE V2\ 9’ | ) ho emc 


d ” 1 (n’ ’ 
ay? Io + Vm:« +(8-9) 

was unter Berücksichtigung von (13) zu ; 

FAS a tle 173 a % 
if AE h® d2 m c? g’ 


führt, wobei dQ” der Raumwinkel ist, in den das Streuquant fliegt. Setzt man in 
e 


Fiir die Dichtefunktion (s. S. 283) erhält man 


1 [73 
ystem 


(12’) ein, so ergibt sich unter Einführung des ‚„Elektronenradius“ 7, = u 
ä = 


Spur Ist n’ die Anzahl der pro sec und cm? einfallenden Lichtquanten, so erhält man 
olari- hieraus für die Zahl dn’’ der pro sec und pro streuendem Elektron in den Raum- 
winkel dQ” in der Streurichtung 9 gestreuten Quanten oe 7 


dn” =n! W— =n — 7 dQ” — +— — sind}. 
Ist J’ = n' cg’ die Intensität des pro sec und cm? einfallenden Lichts, so erhält 
man für die im Abstand R in der Streurichtung # pro sec durch den cm? gehende 
Intensität des gestreuten Lichts @ 
~ Rea’ 


J” =J Ri sin? #). 


v 
J = 


[73 
eg’ n 
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(11) 
urch- 
| 
1 
a 
14) 


Berücksichtigt man, daß nach (13) sich für den Impuls g’”’ und die Frequenz y” 
des Streulichts (»’ Frequenz des einfallenden Lichts) 


, h , 
gmc v v 


I = me+g’(1— =1+y(1—cos8) 


mc 
ergibt, so erhält man aus (14) die Formel von Klein und Nishina!!) 
a n 1 + cos?# y? (1 — cos 9)? 


Durch eine RER Integration läßt sich die Gesamtzahl n’’ der Streu- 
quanten pro sec mer zu 


(15) 


+18% +4y Y—2y—2 
n’ + log (1 + 2y)]. (16) 
Zerstrahlung = Paarvernichtung: Hier ist AA= +1, A=—1. Beide | 


Lichtquanten werden emittiert, also ist in (12) g’ in —g’, g’ in—g”, g’ in—q’, 


q’’ in—g” zu überführen. Für die Dichtefunktion AE erhält man, wenn das Licht- 


quant g’ in den Raumwinkel dQ’ emittiert wird und p, der Impuls des Positrons ist 


1 1 
Mit Hilfe von (13) läßt sich das umrechnen in te en 
AT 
mc? (g’ + g’’) 


womit sich aus (12’) ergibt (ce; = Energie des ran er 
’2 
W = — r2 dQ’ 7 in* 8). 
Ist 2, die Anzahl der pro sec und cm? en Positronen und n* die An- 
zahl der Elektronen in der Volumeneinheit, so erhält man hieraus für die Zahl dn’ 
der pro sec aus der Volumeneinheit in den Raumwinkel dQ’ gestreuten Quanten g’ 
V FR 1 g? /g’ g” 
dn’ n* W = — dQ’ +7 sin?) 

Ist x der Winkel zwischen g’ und p,, so folgen aus Energie- und Impulssatz 
me + me) (€4 + me) — py cos x) 

&, + me — py + mc — p, cos « 


gy = sind = sin x. 


Für die ‚Gesamtzahl n, der pro sec in der Volumeneinheit auftretenden Zer- 
strahlungsprozesse hat man zu bilden 


In’ 
ye n, = 2a | sin « dx, (17) 


wobei der Faktor 4 einzusetzen ist, weil die beiden durch Vertauschen von g’ 
und gw ” auseinander hervorgehenden Fälle identisch sind und nur einmal gezählt 


11), 0), Klein u. Y. Nishina, 2. Physik & 52, u. BER Je. REN = Physik 62, 
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werden diirfen. Die Integration ergibt mit y = = * 
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Durch eine Transformation in das Bezugssystem, in dem ae und Positron 
entgegengesetzt gleichen Impuls vom Betrag P haben, kann man bei einer (pro 
sec und em?) einfallenden Zahl von N_ Elektronen und N, Positronen die in 1 sec 
und 1 cm? auftretende Zahl N, von angeben zu?) 


m? 73 — pa 
| (19) 
Da, a N für P—0 divergiert, gibt mam für kleine P die Zahl N, der Zer- 


strahlungsprozesse besser in Abhängigkeit von der Teilchendichte (Anzahl pro 
em?) N* und N* im Elektronen- und Positronenstrahl'*) 

m? — pa 1+ß gg ‘ 
log — 2(2 — . (20) 
Durch Grenzübergang (P—0, @—mc) erhält man hieraus die in Materie 


von N , Atomen pro cm? zu je Z Elektronen je Sekunde und je (langsamem!) Posi- 
tron auftretende Zahl von Zerstrahlungsprozessen zu 


r ave | 
N,=cN*N* 


N,=N,Zare. (21) 
Paarerzeugung: Hier ist umgekehrt A, = — 1,4 = + 1 zu setzen. Da beide 


Liehtquanten absorbiert werden, ist an den g-Vorzeichen nichts zu ändern. Da 
zur Vereinfachung der Rechnung (s. 8. 288) p, = 0, &, = mc gesetzt worden war, 
erhält man zunächst die Wahrscheinlichkeit dafür, daß zwei Lichtquanten g’ 
und q’’ ein Paar mit ruhendem Positron und einem Elektron vom Impuls p 
im Raumwinkel do bilden. Für die Dichtefunktion erhält man 


1 Vep 
den - =  — dw 
AE (me 4 + + p) Wee 
und aus (12’) - 
emcp 2 
W = ro do - V9’ 9’ (5 8). 


Sind »’ und n” die Zahl der pro’ sec und cm? einfallenden Lichtquanten, so 
ergibt sich für die Zahl dn der pro sec und em? erzeugten Paare mit ruhendem 


Positron und in den Raumwjnkel dw fliegendem Elektron FE 
nV W 1 mep/g 
dn = . ‘= n'n — P + sin?d). 


12) P. A. M. Dirac, Proc. Cambridge Philos. Soc. 26, 361 (1980). 

13) Große Buchstaben N, P, G, beziehen sich stets auf das Schwerpunktssystem der 
einfallenden Teilchen; n, N bedeuten Teilchenströme (pro sec und em?), n*, N* Teilchen- 
dichten (pro em). 

4) Der entsprechenden also 


2 


s. W. Heitler, The atin tt of ss 1936, S. 207 (8), wo e? im Nenner fehlt. 
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Wir gehen iiber in das Bezugssystem, in dem die Impulse der einfallenden Licht- 
quanten und wegen des Impulssatzes auch die Impulse der erzeugten Paare ent- 
gegengesetzt gleich groß sind. Man erhält (neues System in großen Buchstaben, 


& = Winkel zwischen Licht- und Teilchenimpuls im neuen System, G= = des 
einfallenden Lichts, P =/G? — m? c?) 


= G4 — P cost x + 2 m? c? P2 sin? « 


dN = N’ re (@ — P? eos? x)? 


Die Gesamtzahl N aller pro sec und cm? gebildeten Paare erhält man hieraus 
durch Integration über den Raumwinkel zu”) 


_ 1 am „mcPr3— pt 1+Bß 7 
N=—NN an log — 2 (2—B*) | 


B 
6. Prozesse dritter Ordnung 


a) Die Übergangswahrscheinlichkeit 
Die Prozesse dritter Ordnung unterscheiden sich von den vier einfachen Pro- 
zessen zweiter Ordnung (Comptoneffekt am Elektron oder Positron, Paarer- 
zeugung, Zerstrahlung) durch das Auftreten eines dritten Lichtquants.: Da dieses 
entweder absorbiert oder emitt’ert werden kann, gibt es acht verschiedene Pro- 
zesse dritter Ordnung. Bei vier von ihnen (Prozesse A) sind die einfachen Prozesse 
dadurch erweitert, daß ein zusätzliches Lichtquant emittiert wird, bei den rest- 
f lichen vier (Prozesse B) findet der einfache Prozeß unter der zusätzlichen Anregung 
Bei. eines beim Prozeß absorbierten Lichtquants statt. 

¢ Die Ausgangsformeln für die Wahrscheinlichkeit dieser Prozesse lassen sich 
wieder einheitlich gewinnen, wenn man wie in Abschnitt 5 Energievorzeichen i 
einführt und verabredet, daß für Impuls und Energie emittierter Lichtquanten 
—q und —cg geschrieben wird. Man hat dann folgendes Schema 


P 


Anfangszustand des Elektrons 1% ,Ao, do 
4 : p’=p,+q, 4, 0’ 
p’ = Po + G1 + Ge, A”, 
P= Po + Gi + Ge + Gs; A, 0. 


Statt der hier benutzten Reihenfolge g,,Q2,q3 der Lichtquanten kann auch 
eine andere Permutation q,, Q3, Jz bis Gs, Ge, 9 auftreten. _S bedeutet in den fol- 


genden Formeln die Summe iiber diese sche Permutationen. Analog (10) erhält 


man 
RAE V® | 929s | : 
as x, a’) (a’* x, a’’) (a’’* x,a) | * 


(x; bedeutet x e,). Die zweite Klammer im Nenner kann auf Grund des Energie- 
satzes auch Ae —g, — 4’ e’’ geschrieben werden. Im allgemeinen interessiert man 
sich wieder nur dafür, mit welcher Wahrscheinlichkeit der Prozeß beim Aufein- 
15) G, u. A. Physic. 46, 1087 
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andertreffen von Teilchenstrémen geschieht, in denen kein Spin- oder Polarisa- 
tionszustand bevorzugt auftritt. Das erfordert Mittelung iiber diese Richtungen 
im Anfangszustand, Summierung iiber diese Richtungen im Endzustand oder — 
Summierung über alle diese Richtungen und Zufügen eines Faktors f = } für %s 
die vier Prozesse A, f = } fiir die vier Prozesse B. Also \ 


2 6,3 h3 53 


a AE | 919293! 


(a? x, 0’) (a’* | 

TT +9, — (Ae —g3—2" 


Sämtliche Summierungen über Spinrichtungen und Energievorzeichen lassen 


sich wie in Abschnitt 5 durchführen unter Benutzung der Operatoren A= ap _ 
+ßmc. Man erhält dabei 
Te 


& + (de 2}  2he 


„y_% (Ae — gs +h”) +91 th’) 4 0 


_ Berücksichtigen wir Energie- und Impulssatz, so ist 
(Ao + 9)? — = 2 (Ao €0.91 — Po G1) 
(A = =2 (Aegs—P Qs)- 


Setzen wir das in (23) ein und ziehen den Faktor 2*c*/ A, ¢€ heraus, so wird 


W=— es h3 h3 f Vw 


ea es 


9 + 3.04 h 
pur >, (Ao £0 91 — Po 91) (AE Gs — P Gs) 


(Ao & +91 +h’) 
‘(Ag eo +h 
(Ao £0 9ı — Po 91) (AE 93 — P Gs) (Ao + ho) 


Aus dieser Gleichung läßt sich wegen des Auftretens von e® entnehmen, daß sich 
die Wirkungsquerschnitte zu 


— (<) (dimensionslose Kombinationen von m c, 91; 92; 93; Po» P) 
(25) 
ergeben werden. Sie sind also wegen des Auftretens des Faktors =— =. im 
allgemeinen kleiner als die Wirkungsquerschnitte der entsprechenden Prozesse 
zweiter Ordnung. Es kann nun eintreten und tritt auch ein, daß für bestimmte 
Größenordnungen der Energien und Impulse diese Verringerung auf den 137. Teil 
durch den dritten Faktor in (25) wieder ausgeglichen wird. Weiter kann man aus 
(24) ablesen, daß eine allgemeine für alle Größen der Teilchen- und Lichtquanten- 
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impulse gültige Berechnung wegen des Auftretens von 36 Spuren ($ mit sechs 

Summanden - $ mit wieder sechs Summanden) nur schwer möglich ist. Auch die 

® 

im Abschnitt 5 benutzte Vereinfachung p, = 0, & = mc macht die Berechnung 

noch nicht wesentlich leichter; sie gelingt erst, wenn man noch voraussetzt, daß 

zum mindesten ein Teil der auftretenden Impulse groß gegen m c ist. 


b) Bereits berechnete Wirkungsquerschnitte und Übergangswahrscheinlichkeiten 


Bisher ist die Wahrscheinlichkeit ihres Auftretens berechnet worden für zwei 
der Prozesse A, bei denen ein zusätzliches Lichtquant emittiert wird, und zwar: 

a) A, =A= +1, die „Doppel-Comptonstreuung“, d.h. ein Comptoneffekt, 
bei dem das einfallende Lichtquant zwei Streuquanten bildet !®). 

b) A, =—1, A= +1, d.h. die Erzeugung eines Elektron-Positron-Paares 
unter Absorption zweier Lichtquanten und gleichzeitiger Emission eines dritten 
Lichtquants!?). Noch nicht berechnet ist der Fall 

die „Dreifachzerstrahlung“, d.h. die Zerstrahlung 
eines Paares unter Bildung von drei Lichtquanten. Für diesen Prozeß soll im fol- 
genden die Wahrscheinlichkeit angegeben und mit der der Zerstrahlung in zwei 
Lichtquanten verglichen werden. 


c) Die Dreifachzerstrahlung 
Wir gehen aus von (24) und fiihren die Rechnung in dem Bezugssystem durch, 
in dem das Paar entgegengesetzt gleichen Impuls hat. Wir haben also 


Po= VP. 


=m ne Zugleich berücksichtigen wir, daß alle drei Lichtquanten emittiert werden, 


dom also in (24) g, g in —q, —g über. Wir erhalten so (mit f=}, A, = +1, 

i=—]) 


(€ 9, — P G1) (E93 + P Gs) 

Wr (€ = 93 — (e—g +h’) oy 
“3” (€ 91 — P G1) (€ 93 + P Gs) 
woh=ap+ßmeh"=ap—ag +ßme,h’"=ap+ag +Pßme ist. 

Der Gesamtimpuls des Paares ist gleich Null, die Energie ist 2ce, also muß 

91 + G2 + 93 = 9, 9; +92 +93 = 2e sein. Die Impulse der Lichtquanten können 

also in einer Ebene beliebiger Normalenrichtung liegen; sie miissen sich nur zu 

einem Impulsdreieck schließen (Impulssatz), und der Umfang dieses Dreiecks muß 

2¢ sein (Energiesatz). Für eins der Lichtquanten, etwa g,, sind Impulsrichtung 

und -größe frei wählbar, für ein zweites BERNER, etwa g,, ist noch die Richtung 


‘(e +h), 


frei wählbar. Für die Dichtefunktion —, erhält man 


ae 
16) W. Heitler u. L. Nordheim, Physica 1, 1059 (1934). a on pay, 

B. Kockel, 1. ce. 
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sechs $# also unter Benutzung von Energie- und Impulssatz 
h die 91 9 
= m 94919540, (€ — 9s) 
nung Man erhält an 
1 3 3 
aten Sind N_, N, die Zahl der pro sec und cm? einfallenden Elektronen und Posi- 
zwei # tronen, so läßt sich die Zahl dN derjenigen pro sec im cm? eintretenden Dreifach- 
war: | zerstrahlungen, bei denen eins der entstehenden Lichtquanten einen Impuls 


fekt, zwischen g, und 9, +dg, hat und in den Raumwinkel dQ, um g,, ein wei- 
teres Lichtquant in den Raumwinkel dQ, fliegt, angeben zu 


ares 1 Ar 1 e? 2 e2 91 93 m? c? g, < i: 
dl = —N- Pr (5) D> Spur... 


itten + he p® 2e(e 

lung Die Gesamtzahl von Dreifachzerstrahlungen pro &ec und cm? bekommt man — 

fol- hieraus durch Integration 

zwei 


| 102422 (me na) “pF Beleg) 2 2 
Der größte und kleinste mögliche Wert für den Impuls eines Lichtquants sind 0 
und e. Der dem Faktor $ in Gleichung (17) entsprechende Faktor 4 ist beizu- 
fügen, weil die Fälle Emission von g,, 2, 9, von Qi = 91, 92 = Gs» 95 = Je usw. 
nur einmal zu zählen sind. 

len, Für die Berechnung der Gesamtzahl von Dreifachzerstrahlungen ist es nun 
1, vollständig gleichgültig, welche Lichtquanten wir bei der Raumwinkel- und Im- 
pulsintegration bevorzugen. Wir können also auch schreiben unter Ziehen einer 


rch, 


der $ bis vor das Integral =, 
e: 3 (27) 


(€ 9, —P 91) (€ 93 + P Gs) (e 3 (€ 9; —P Gi) (E93 + P Gs) .+m- 


Die erste = gibt einfach einen Faktor 6, also + Pr 
uB 
len 19, meg 
uB (28) 
ng d) Die Dreifachzerstrahlung bei hoher Energie der Teilchen 


Die Berechnung von (28) ist nur schwerlich allgemein möglich. Sind jedoch 
die Energien des Elektrons und des Positrons groß gegen m c?, so ergibt sich eine 
Möglichkeit zu einer genäherten Berechnung auf Grund einer Eigenschaft der 
Nenner in (27). Jeder oe 


egtpg a))} fiir &p>me 


Ann. Physik. 6. Folge, Bd. 4 ve 


i 
> 


* hat die Eigenschaft, wenn g antiparallel bzw. parallel zu p ist, sehr klein zu werden. 
Man erhält also eine brauchbare Näherung, wenn man nur für die Nachbarschaft 
solcher Lagen die Spur in (28) exakt berechnet. In Anhang I wird gezeigt, daß 
_ man dabei nur zwei von den sechs Spuren zu berücksichtigen braucht. In Anhang II 
und III werden diese Spuren berechnet. Addiert man ihre Beiträge, so erhält 
.man nach (27), (28), (39) und (42) 


1 gi m Is 
3 
1 (e — g,) sin? (p gs) 
(p + + p cos (p 9,))? (28) 
p sin d, sin d, cos (9, — 93) } 
2 gs (€ — 9s) (€ + pcos (P Gs) (€ + pcos (P 9,)) 


(8, Bs, Po, Ys siehe Anhang III). Der stärkste Beitrag von der Nachbar- 
schaft der Lage: g,, gs antiparallel zu p,g, parallel zu p, (9, & p we) her. Nur über 
die Nachbarschaft dieser Lage wird integriert. 2. die in ie eckigen Klammer 


stehenden Ausdrücke ist das einfach. Mit = = ergibt sich 


(sna 1 fede _ 4a 


— N 
P. q 


Entsprechend ergibt sich 
[ 1 p* (© — gs) sin? (p gs) 
+ pcos (P gs) (+ 

Ne “Also hat man zunächst als Beitrag zu (29) 


Raumwinkelintegral noch Anteils (auBerhalb der eckigen Klam- 
A ist in Anhang IV berechnet und ergibt als Beitrag zu (29) = 


Also ist schließlich 


rein mick P_ [9349s P(pP—9s) _P 
N c N N: ic p log* mc fo +4 5 
0 


LP . 
was auch vereinfacht werden kann zu?8) 


4 mie Pp. Gs , ple-g) 
N= c (<4) hc p? log* + 93 93 (€ — 95) } (30) 


= 18) Entsprechned dem Rechenverfahren fehlen in (30) Ausdrücke, die log pro- 
portional sind; mit anderen erg in (30) zu noch ein Zahlenfaktor a der Größen- 


ordnung 1 in der Form log? statt log? — auftreten. 
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Bilden wir hieraus die Gesamtzahl der Dreifachzerstrahlungen, so finden wir, 
daß das Integral über dg, divergiert sowohl an seiner unteren wie an seiner oberen 
Grenze. Diese Divergenzen sind nicht etwa fiktiv, also nicht nur durch die Ver- 


nachlässigung des Rechenverfahrens entstanden). Beide Divergenzen haben ihre ° 


Ursache darin, daß Lichtqhanten beliebig kleinen Impulses auftreten können. 
An der unteren Grenze des Integrals kann g,—0 gehen; an der oberen Grenze, 
wo 9, % eist, kann einer der beiden anderen Lichtquantenimpulse beliebig klein 
werden. Divergenzen dieser Art sind, als ,,Infrarotkatastrophe“ bezeichnet, schon 
öfter aufgetreten, so z.B. bei der Behandlung eines anderen Prozesses dritter 
Ordnung, der Paarerzeugung unter gleichzeitiger Emission eines Lichtquants 
(Prozess A, = +1, A=—1 von Seite 294) oder der Behandlung der Brems- 
strahlung**). Die Auflösung dieser Schwierigkeit liegt in folgendem: Die genannten 
Prozesse können geschehen unter Emission weiterer und zwar beliebig vieler 
Lichtquanten*!), und nur die Intensität des insgesamt emittierten ‚Lichts bleibt 
endlich **). Im Grenzübergang g, — 0 bleibt also die Intensität endlich®), die Anzahl 
der Lichtquanten geht jedoch gegen unendlich. Damit wird die übliche Methode, 
die Prozesse nach ihrer Ordnung, d.h. der Anzahl der beteiligten Lichtquanten 
einzuteilen, illusorisch. Berechnungen, die dieses Verfahren einschlagen, können 
also die Wahrscheinlichkeit dieser Prozesse nur für den Fall angeben, daß die be- 
teiligten Lichtquanten einen hinreichend von 0 verschiedenen Impuls besitzen. 
Die Gleichung (30) kann also entsprechend nur aussagen, wie groß die Zahl der 
Dreifachzerstrahlungen unter Emission hinreichend kurzwelliger Lichtquanten ist. 

Solange noch keine Theorie vorliegt, die es ermöglicht, Prozesse mit einer be- 
liebigen Anzahl beteiligter Lichtquanten ohne die Entwicklung nach dieser Anzahl 
zu berechnen, bleibt nichts anderes übrig, als (30) mit einem Abschneideverfahren 
zu behandeln. Wir integrieren über dg, also nur von mc bis e— mc und erhalten 


| 
e) Vergleich: Zweifachzerstrahlung-Dreifachzerstrahlung 

Die (31) entsprechende Zahl von Zweifachzerstrahlungen, d.h. die Anzahl 
von Zweifachzerstrahlungen im cm® bei einer einfallenden Zahl von N_, N, Elek- 
tronen bzw. Positronen pro sec und cm? erhält man aus (19), wenn man dort ß—1 


gehen läßt und nur die höchste, hier lineare Potenz von log Pe aufschreibt zu 


P= e? \2 m? c? p 
N =—N_N,(-5) log (32) 

19) Rechnet man die Spuren und die Integration über die Raumwinkel exakt fiir 
9,» 0 und g, » € aus — im letzteren Fall führt man am besten elliptische Koordinaten 
mit g, als Verbindungslinie der Brennpunkte ein — so findet man wieder die Proportionali- 
tät zu oder 

20) N. F. Mott, Proc. Cambridge philos. Soc. 27, 255 (1931); A. Sommerfeld, Ann. 
Physik 11, 257 (1931); H. Bethe u. W. Heitler, Proc. Roy. Soc. London A 146, 83 (1934). 

2!) Der Energiesatz muß natürlich gewahrt bleiben. 

2) F. Bloch u. A. Nordsieck, Physic. Rev. 52, 54 (1937). 

23) Die Intensität J (v) kann auch unendlich werden, wenn die Frequenz » — 0 geht. 
fJ (v) dv bleibt natürlich endlich. Das ist z. B. der Fall bei der Bremsstrahlung, die ein 
Elektron beim Durchgang durch ein reines Coulombfeld aussendet. 


20* 
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e2 
Ke’ also 137mal 


Die Dreifachzerstrahlung ist also zunächst um den Faktor 
seltener. Für sehr hohe Energien kann dieser Faktor jedoch durch die höhere 
Potenz des log > ausgeglichen werden. Wegen der Unsicherheit in der Berech- 


nung von (31) und der geringen Variabilität des Logarithmus kann man aber die 
Energie, für die das eintritt, nicht aus (31) und (32) entnehmen. Die folgende Ab- 
bildung zeigt die Wahrscheinlichkeit der Zweifachzerstrahlung (Kurve A nach 


Er: . Formel (19), Kurve B nach Formel (20)) und der Dreifachzerstrahlung (Kurve ( 
> nach Formel (31)), wobei im letzteren Fall eine Unsicherheit um den Faktor yı 
Br. nach oben oder unten angenommen wurde. Die Abbildung ist weit über den in- 
Es teressierenden Bereich der Teilchenenergie von etwa 0,5 bis 1000 MeV hinaus- 


gezeichnet, damit die bei sehr hohen Energien auftretende größere Wahrschein- 
lichkeit der Dreifachzerstrahlung sichtbar wird. (Über die Gültigkeit der Ergeb- 
‚nisse in dem Gebiet gleicher oder größerer Wahrscheinlichkeit der Dreifachzer- 


Br f strahlung vgl. Seite 282). Da in dem experimentell erreichbaren Energiebereich 
ar die Zweifachzerstrahlung häufiger eintritt als die Dreifachzerstrahlung, ist an 
"a den Folgerungen, die man nach Berechnung der Wahrscheinlichkeit der Zweifel- 


zerstrahlung zog, nichts zu ändern. Sowohl die Lebenszeit langsamer Posi- 
tronen wie der Anteil der Zerstrahlung an der Abbremsung schneller durch Materie 
fliegender Positronen bleiben weiter durch die Zweifachzerstrahlung bestimmt. 


A 
An. für Zweifachzerstrahlung (rechte Orinaten) 
» » (linke Ondinaten) 
=. 
| Se 
we 
: 


Anhang I. Die Spuren bei der Dreifachzerstrallung 

Die Nenner dieser Spuren sind nach (27) von der Form IE 
91939. 9, (€—P (p g,)) (€ + pcos (p ga)) (e+p cos(p (e—peos(pg,)) (3) 
wo 1, j, k eine der Permutationen von 1, 2, 3 ist. Jeder der Klammerfaktoren 
wird wegen p> mc, e = p fast = 0, wenn g, oder g, fast parallel zu p, g, oder & 
fast antiparallel zu p werden. Für die Nachbarschaft dieser speziellen Lagen kann 


man mit = ‚q<.1 schreiben 
e + pcos (pg) = 4 (g? + sin? (p g)). 
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Die Zähler Z haben die Form 


Spura,Ba;A =X e;; 


| 

(35) 


A=(e+ap+fmc —g,; (e +ap + Bme —g, 

(36) 

Wir untersuchen A weiter. Unter Benutzung von (e+xp+Pßme) 
-(e—a« p—B mc) = 0 findet man 


A= {x; (9; §:) — 2 (p e,)} (e +ap+Pßme) { (91 91) — 2 (p e;)}- (37) 


Wenn g, fast parallel zu p ist, so ist (p e,) proportional zu sin (p 9,), 9, —« 9; = 
gi (1 — «,) + zu sin (pg,) proportionalen Termen (z-Richtung sei die Richtung 
von p). Füre+ap-+ßmc kann man wegen p > me schreiben p {(1 + «,) + 
Bqa+4q*}. Da (1+ .«,)(1—a,) = 0 ist, folgt, daß das Produkt der 2. und 3. 
Klammer in (37) sich zusammensetzt aus Ausdrücken, die zu 8 g, 9? oder sin (p g;) 
proportional sind, wenn ¢ — p cos (p g,) fast = 0 wird. 

Wenn g, fast parallel zu p und damit der Nennerfaktor e — p cos (p g,) fast = 0 ir 
wird, gilt das Entsprechende. 

Werden sowohl g, wie 9, fast parallel zu p, so findet man auf die gleiche Weise, - 
daß dann A sich zusammensetzt aus mindestens zu sin (p g,) sin (p 9,), 8 q sin (p 9,), 7 
Bqsin (p g,) oder g? proportionalen Ausdrücken. ' 

Fast man diese Ergebnisse zusammen und berücksichtigt, daß für B dasselbe \ 
gilt, so findet man, daß für jeden fast = 0 werdenden Nenner im Zähler eine Pro- 
portionalität zu fq oder sin (p qo), oe = 1, 2 oder 3 auftritt. Prüft man daraufhin 
die Spuren in der Hinsicht durch, ob sie bei der Winkelintegration über d2, dQ, 


einen zu log! proportionalen Anteil geben können, so fipdet man, daß das nur 


auftreten kann, wenn alle vier Nennerfaktoren fast = 0 werden. Das tritt ein, 
wenn zugleich g, und g, parallel oder fast parallel zu p, g, und g, antiparallel oder 
fast antiparallel zu p sind. Das wiederum ist nur möglich für die symmetrische 
Spur mit k = 3, 7 = 2,7 = 1 und für die Spur mit k = 2,7 = 3,1 = 1. 


Anhang II. Die symmetrische Spur 


In (37) setzen wir 2 = 1 und werten diesen Ausdruck noch exakt aus zu 


A= 91) (G1 + Gi — 2 (P er) + 9, +ap+ fame). 


Da die beiden Faktoren ¢g,— pg, und (pe,)? zu sin? (pg,) proportional sind, 
setzen wir für die anderen Faktoren die Werte ein, die sie für die exakte Parallel- 
lage: g, ||p annehmen und erhalten unter Berücksichtigung von ex p 


A = len + 4 (p ex)? —g1)} (1 +0). 


Für B erhält man entsprechend 


B= {2 (eg, +492) go +4 (pes)* (I 
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Die Spur ergibt sich unter Berücksichtigung der in der Parallellage eintre- 
tenden Antikommutativität von x, und «, zu 


= Spur a, Ba, A = 82 {(eg,—pgi)gi. + 2(p e:)? —g:)} 
+P Gs) Js +2 (p es)? (e—9)}- 


Da in der Parallellage auf Grund des Energiesatzes g, = & sein muß, verein- 
facht sich das zu 


Spur a, B,A=32e(eg, —Pp9ı) {(e gs +P9)9 + 2 (p es)? (e — Ys) }- 
Anbringung der zugehörigen Nennerfaktoren und Summierung über e,, &,, & 
ergibt schließlich 

€9, —PG + PGs (€ 93 + P 93)? 


oder 


1 1 (e —g,) sin? (p gs) 
- 256 — | —_______ 

| — pcos (p 91) Le + pcos (p gs) gs (€ + pcos (p (39) 
N Anhang IH. Die Spur k=2, j=3, i=1 

Wir gehen von (34) bis (38) aus. Für A erhält man denselben Wert wie in An- 


hang II, also nach (38) unter Berücksichtigung von & x g, . 


für B erhält man entsprechend (37) 
B= {2 (pes) + (9s +. Qs) x3} (e p —B mc) {X (92 +o 92) +2(pe,)}. (41) 


Wir setzen, um die Verhältnisse in der Nachbarschaft der speziellen Lage 9; 
und g, antiparallel zu p zu erfassen 


9s + 9s, = 963 E—P= z-Richtung = p-Richtung 


WO €2,, €3, die zur z-Richtung senkrechten Komponenten von e,, e, sind. Dabei 


sind h,, /; proportional zu sin (p g,), 6; proportional zu sin? (p g,) usw. Weiter kann 


man in (41) 8 mc streichen, da es keinen Beitrag zur Spur geben kann, nachdem 
A (s. (40)) B-frei geworden ist. Man erhält unter Streichen von zu sin® (p g,) oder 
sin? (p g,) proportionalen Ausdrücken 
B= {2 (pes) + (1 + fs) +60, Ca) +04 (0 +, fs)} 

{2 (p €2) + 92 (x Cos + fe) (1 + (x C25) &, + (x Ces + fe) He} 


oder, wieder unter Weglassen von allen mehr als in zweiter Ordnung kleinen 
Ausdrücken und Berücksichtigung der Antikommutativität von a, mit & eg 


Fi 


Spur 


D 
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sichti 
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jetzt : 


und ] 
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B ist jetzt in zweiter Ordnung klein, also brauchen wir bei der Bildung von 
& B xs nur 9 = & % = & Cs, zu setzen und erhalten 


0, Bx, = [2 p Js % Cog X XK + 4(p C3 —Gs fs) (P fe) x Cog 
+ 2 (pes fs) x Hex + 2 (P C2 fe) C2, Oy 
Für «, B a; A ergibt sich unter gleichzeitiger Spurbildung ta 
Spur x, Ba, A = 16 (eg: —P (2 (€25 83)? — 1) 
' 
+ 4 (p es —gs fs) — 92 fe) (Cas Css) 
dj + 2 (p es —gs fs) (Be + 2 —Ge fe) (Bs Cos) + (Be Hs) (Es C35)? 
+ (Gs Cas) + ess) (Bs ss) — (Bo Cas) (Bs Cas) as) 
— (He 25) (Bs a5) (C25 es.) | . 


Die Summation über die Polarisationsrichtungen ergibt (#2, Js, G2, 9, polare 
und azimutale Winkel der g,, g;-Richtung — p-Richtung als Aghse) 


X Spur a, Ba; A= 128 p? (eg, — p gi) sin sin cos (9 — gs). 


Nach Anbringen der zugehörigen Nennerfaktoren ergibt sich, wenn wir berück- 
sichtigen, daß in der Nachbarschaft der speziellen Lage ¢ © p & 91.92 & € — 9; ist 


222% Spur, Ba; A 


® 
Gi 92 Js (€ — p cos (P g,))*(e + p cos (Pp + p cos (P gs) 
128 1 sin d, sin cos — 95) 


(€ — G3) € — p COs (P 91) (€ + P cos (P Ge)) (€ + COs (P gs)) 
pr 


Anhang IV. Die Raumwinkelintegration für die Spur k=2, j=2, i= 1 
Nach (29) ist zu bilden 


1 p° sin 9, sin 9, cos (Pa — 95) 

J= > dQ, dQ. 2 
2 gs — gs) (€ — p cos (p gi) (€ + Pp cos (P (€ + pcos (Pgs) 

Wir drücken den Integranden durch #,, #3, 9, und g, aus und zählen dabei 


jetzt 9, und 9, so, daß sie in der speziellen Lage selbst = 0 sind. Aus Energie- 


und Impulssatz folgt für die Nachbarschaft dieser Lage mit g = — 


p 
2 (€—gs)* 


{q (e —gs)* + p? + 0,9; cos — 


€ + pcos (p gs) = 


1 
sind, sind, cos —Ys3) Is (9s ds; cos (91 —9))- 


2, 6 | 
(40) q 
(41) 4 
x, fe 3 
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Außerdem schreiben wir natürlich wieder 
3 


| 
e— pcos (pg,) = 4 +9) e+pcos(p gs) =F(¢+%) 


und bekommen 


J= 2 93(€ €— Js 9s) 
of 9; 03 — pd, 9; cos — Ys) 


oder nach leichter Umformung und Durchführung der Integration über die azi- 
 mutalen Winkel 


wl 

enc 
oes ‘pe p — (e — 9)? Dif 
V tg? (© — 95)" + p* + gi — 4 85 daf 


u 
wobei die Integration nur iiber die Nachbarschaft der speziellen Lage durch die geg 


Integrationsgrenzen 0 bis ungefähr 1 ausgedrückt sei. Die eckige Klammer ist 80 * 
for 

[.--] für pd, < 9395. (43) lei 


Damit ergibt sich vol 


(44) 
mA, 
„Ir 
[ dr | dy (7 + x) + y) 93 | tats 
kei 
= -f2 log a + Terme, die nur einem log | Prop. sind +.. 1 ; 


_ Die Fehler, die die Vereinfachung (43) mit sich bringt, sind ohne Einfluß auf 
das Ergebnis (44), weil sie nur in einem Gebiet 49,49, der Größenordnung ¢ wo 


auftreten. Es ist also unter Vernachlassigung aller nicht zu log? - proportionalen 


Jog? 


mc” 


p 

Termeı 

Be. Leipzig, Institut für theoretische Physik der Universität. i die 
ean (Bei der Redaktion eingegangtn am 28. Dezember 1948.) lie 
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Über die Berechnung des Schallfeldes von kreisförmigen on 
Membranen in starrer Wand 


Von H. Stenzel 
(Mit 3 Abbildungen) 


Inhaltsübersicht 


Zur Bestimmung der für den Schalldruck im Aufpunkt (r, y) charakteristischen 
Funktion R (r,y) werden die speziellen Funktionen R (00,y) (Aufpunkt im Un- 
endlichen) und R (r, 0) (Aufpunkt auf der Mittelachse) eingeführt. Mit Hilfe der 
Differentialgleichung für die Legendreschen Kugelfunktionen P,, (z) wird gezeigt, 
daß die Koeffizienten der nach fallenden Potenzen von r fortschreitenden Reihe 
ven A (7, y) durch eine auf die Funktion R (00, y) ausgeübte Differentialoperation 
gegeben sind. Ist insbesondere R (00, y) als Potenzreihe von z (z = cos y) gegeben, Bi 
so wird i (r, y) ebenfalls durch eine Reihe dargestellt, die nach Potenzen von cos y aa 
fortschreitet. 

In analoger Weise ergibt sich mit Hilfe der Differentialgleichung der Ray- 
leighschen Funktion f, (¢kr) die Bestimmung der Koeffizienten der nach Potenzen 
von sin y fortschreitenden Reihe von  (r,y) durch eine auf die Funktion R (r, 0) 
ausgeübte Differentialoperation. 

Die Tragweite der abgeleiteten Formeln wird an Hand von Beispielen für die | 


1. Die Formulierung der Aufgabe 


Ausgehend von der Rayleigh schen Fundamentalformel für das Geschwindig- 
keitspotential y einer in starrer Wand schwingenden Membran: 


wollen wir mit: ze w (o) ei”! den im Aufpunkt A BEE Schalldruck 


= icox® in der Form darstellen: 


2acoF w,, 
“r Win 


-i(w t+2/2-xr) 1 (r-R) 


w (0) dF. 


Dabei ist Rotationssymmetrie fiir die schwingende Membran vorausgesetzt, so 

daß wir ohne Beschränkung der Allgemeinheit den Aufpunkt in der YZ-Ebene 

des Koordinatensystems annehmen kénnen, dessen Mittelpunkt O mit dem Mem- 

branmittelpunkt zusammenfällt, während die -Membranebene in XY- ‘Ebene 


F 


h die 
er ist 
(43) 
= 
= | 


Ferner bedeutet: 
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der Schallwiderstand des Mediums 


die zur Wellenlänge A gehörende Wellenzahl 
der Zeitfaktor der fortschreitenden Welle 
die Koordinaten des Strahlerelementes dF inder XY -Ebene 


die Membranschnelle für das Strahlerelement dF a 
die mittlere Membranschnelle 


der Abstand des Aufpunktes A vom Strahlerelement dF 
die Koordinaten des Aufpunktes in der YZ-Ebene 


die zu beliebigem w (9) gehörende charakteristische Funktion 
des Aufpunktes (r, y) 


die zu w (0) = (1— 0*/a?)P-! gehörende charakteristische 


Funktion des Aufpunktes (r, y). P 


‘Wir definieren die relativen (dimensions- 
ee losen) Druckkomponenten p, und p,, durch: 


+t Pm) 


und erhalten dann: 


Abb. 1. Zur Berechnung des Schall- 
der kreisférmigen Membran rechnufk von § (r, y) bzw. Pa +t Pm er 


Nun ergibt sich aus der allgemeinen Theorie!), daß für Aufpunkte, deren Ab- 
stand vom Membranmittelpunkt gréBer als der Membranradius ist, — nur solche 
Aufpunkte sollen im Folgenden betrachtet werden — das Geschwindigkeitspoten- 
tial durch die nach Kugelfunktionen fortschreitende, absolut und gleichmäßig 
konvergente, Reihe dargestellt werden kann: 


1) H. Backhaus, Ann. Physik 5, 3 (1930). Da wir nach Sommerfeld e** als 
Zeitfaktor verwenden im Gegensatz zu e**®* bei Backhaus und Rayleigh, müssen wir 
; (an Stelle ae ; bei Backhaus) schreiben. Aus demselben Grund ist auf 
der rechten Seite von (6) in der Rayleighschen Definition + i xr durch —i xr ersetzt. 


Pa +i Pm 
x 


mit: 


F 
Unsere Aufgabe soll es dann sein, für R 
Reihenentwicklungen anzugeben, die insbe- 
sondere für große xr eine praktische Be- 


möglichen. 


co 
N" tan Pon (cosy). 


n=0 


H.St 


Ans 


die v 


und 


| 
| 
F=na! die strahlende Kreisfläche mit dem Radius a, 
x = 2n/i 
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An Stelle der Hankelschen Funktion HP. 3(* 7) führen wir durch die Beziehung?) 


Has (kr)=(- 0) mtleixr (ixr) =S,, (xr) —iC,, r) (5) 


die von Rayleigh®) benutzte Funktion ein: 
- 
PR . m iad + 1) je (m — 1) m (m + 1) (m + 2) 7 ur, 
fm =1 +i +i | 
jm = 


und erhalten dann für R (r,y) die folgende Darstellung: 


oo 
(r, y) (x a) fon x r) Pon (cos y). (7) 


Bekanntlich *) läßt sich für die Kolbenmembran (w (9) = 1) die Funktion R (r, y), 
wenn y = 0 oder r = 00 ist (d.h. der Aufpunkt auf der Mittelachse oder im Un- 
endlichen liegt), sehr einfach durch direkte Auswertung des Integrals in (3) finden. 
Das gilt auch für allgemeinere Schwingungstypen w (0) = (1— o?/a?)P-!, wie wir im 
nächsten Abschnitt erläutern werden., Wir können deshalb die entsprechenden 
Funktionen R (r,o) und R (r,y) als bekannt ansehen. Für die Kolbenmembran 


92 inv 
ergibt sich für großes r wegen R (co, y) = = oes = ”) aus (2) die bekannte 
Naherung: 


+ipt = — (za)? 2J,(xasiny) 


' 

Es liegt dann nahe, nach einer allgemeinen Reihenentwicklung von i (r,y) zu 
fragen, bei der das erste Glied durch R (r,o) bzw. durch i (00, y) gegeben ist. 
Wir werden damit auf Reihen fiir R (r,y) geführt, die nach steigenden Potenzen 
von sin y oder von cosy bzw. nach fallenden Potenzen von r fortschreiten, und 
stellen uns die Aufgabe, die Koeffizienten dieser Reihen zu bestimmen >). 

Das Hauptergebnis der Arbeit besteht dann in dem Nachweis, daß die gesuchten 
Koeffizienten sich unmittelbar durch Differentialoperationen, die auf R (00, y) 


*) W. Magnus u. F. Oberhettinger, Formeln und Sätze..., 2. Aufl. Berlin, Sprin- 
ger 1943, S. 32. 

3) Lord Rayleigh, Theory of Sound, London, 1929, $ 323. 

4) H. Backhausu.F. Trendelenburg, 2. techn. Physik 7 , 630 (1926); H. Stenzel, 
E.N.T. 4, 251 (1927). 

5) In ‘ähnlicher Weise wird in der Arbeit: A. O. Williamsand L. W. Labaw, Acoustic 
Intensity and Distribution from a Piston Source [Acoust. Soc. of America, Vol. 16, 231 
(1945)] eine Näherungsformel für Aufpunkte, die in der Nähe der Mittelachse ‘und in groBer 
Entfernung liegen, abgeleitet, doch ist das Ergebnis fiir den absoluten Schalldruck: 


2 
| Po +i Pm |= 28in (1 


wenig befriedigend,da es wegen r/z= tg ~siny sich praktisch von der bekannten obigen 
Näherung: 


(za)? 2 J, (za siny) zxa® ( ) 
| Pa + Pm | = xasiny  2r * 
kaum unterscheidet. ® 


| | 
tion 
sche q 
. 
yns- 
ch: | 
all 
(3) 4 
R 
be- 
Be- 
er- 
{ 
she 
en- 
Big = 
5; 
wir 
ruf : 
zt. 
— - 
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bzw. auf Ri (r,o) auszuführen sind, bestimmen lassen. Bevor wir diesen Nachweis 
führen, sollen im folgenden Abschnitt die für R(co,y) und R (r,o) geltenden 
Entwicklungen zusammengestellt werden. 


2. Die Darstellung der Funktionen R (co, y) und R (r, o) 
Mit Hilfe der Besselschen Funktion nullter Ordnung: 


22 
; 


1 
ixosin y cos 
| € dp 


ur 

geht für r = oo die Formel (3) über in: | 


R (co, y) = J, (zo sin y) - w (0) odo. (8) 


Setzen wir für J,(xosiny) die bekannte Reihe ein, welche nach Bessel- 
schen Funktionen höherer Ordnung J,, (x 0) und nach Potenzen von z (2= cosy) 
fortschreitet ®): 


2 
Jo sin y) = Jo (x9) + (%0) N) 


so folgt: 
KR (co, y) = + 
wobei B,, durch: 
22 2 
m! 2m J (#2 0) w (0) odo 
ö 


definiert ist. 
Bezeichnen wir in (3) die zu w (o) = (1 — 9?/a*)?~1 gehörende Funktion mit 
RP (r, y), so ist: 


J, (xasiny) 7 
RIP] (co, y) = 2 p! (12) 


Daraus folgen die beiden Darstellungen: 


2? z z4 ( P 
WP (00,9) = ly, (za) + Sper (za) +I (xa) +4 (13) 
. (x a)? sin? y (x a)* sinty 
Für y = 0 geht die Formel (3) über in: er 
R (r, 0): Wy, F J Vr + p? Pr (e) de 


0 


6) Watson, Theory of Bessel Functions, Cambridge 1929, S. 141. 
7) H. Stenzel, Ann. Physik 7, 972 (1930). : 
8) G,N. Watson, Theory of Bessel Functions, TE 1929, S. 141, dabei ist 
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ziel 
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nden 


(14) 


(15) 


i ist 


3. Die Darstellung von R (r, y) mit Hilfe von R (00, y) 


Für die Kolbenmembran w (6) = 1 folgt daraus: 


Qixr 
(xa)? 


Ku (r, 0) = — 


Für beliebiges w (0) benutzen wir zur Auswertung des Integrales (15) die Be- 
ziehung ?°): 


re'*r > i" (x 0)" f, (i xr) (17) 
Vr? + a=0 
Setzen wir dies in (15) ein, so folgt: 
0) — 1+ “BESO (18) 


2xr (2 zr)? 


» 

- 5; n! (x w (0) odo. (19) 

It speziell w (0) = (1 — o?/a?)P-1, so ergibt sich: ¢, = 


und es folgt: 
RiP] (r,o) = 1+ 


(p +1)(p+2)...(p tn) 


tf, (i xr) (ka 


(20) 


nach fallenden Potenzen von x r 
a) Die Darstellung der Entwicklungskoeffizienten durch einen Differentialoperator 


Setzen wir für fa, (i 17) die sich aus (6) ergebende Reihe in (7) ein, so können 
wir R (r, y) in der Form schreiben: 


= Sol) + (21) 
wi: (2n + m)! 

on (2 —— P..(z Bs 29 

= m (z) (xa) (2 n— m)! F 2n (z) 22 


und insbesondere: 


Se (2) = don (ea) QB) 


= Asn (xa) 2n (2 + 1) Pen (2) (24) 


r= 


Wir zeigen zunächst, wie wir den zweiten Koeffizienten in (21): 
g 


®) H. Backhaus u. F. Trendelenburg, Z. techn. Physik 7, 630 (1926). 
10) H. Stenzel, Ann. Physik 41, 253 (1942). Dies ergibt sich aus der Formel (43), 
wenn man darin x o durch x r und — (x a)? sin? y durch (x 0)? ersetzt und die Beziehung: 


Cp (#1) +18,, (% (— i)" 
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durch ©, (z) bestimmen kénnen. Dazu benutzen wir die ee Differential- 
gleichung der Kugelfunktionen: 
dz 
Den auf der linken Seite von (25) ci Ausdruck fassen wir als Operator D, 
von P,, auf und schreiben: 


2-1) n(m+1)P,@). (25) 


dP,\ 
dz 


Wenden wir dann die Operation: 
‘Dy [Pan] = 2n 22+ 1) Pan 
auf beide Seiten der Gleichung (23) an, so folgt: 
D, [So (2)] = Ei (za) 2n (2n 4-1) Poy (2). 
Das heißt, es ist nach (24): 
Legen wir z. B. die Kolbenmembran -. w (o) = 1) zugrunde, so folgt aus: 


AN (00,y) = (2) = 2 1-2?) 


(26) 


durch Differentiation 
(z2—1) = —2zJ, (za yı — 22) 
und daraus mit der Abkürzung zay 1 — 2= xa siny = u EN Me 
3% 
2(1 + 27) J, (u) —2 22 uJ, (u) I 
Um &, (z), &; (z) usw. zu finden, definieren wir die ji ratoren: bd 
D,[P)=D,[D, [Pal] — 


D; | [D. — [Pa] 


. 
Im [P,]= D, ilP,]] m (m — 1) 1 


Dann ist, wegen n(n + 1)—1-2= (n—1) (n + 2), ur 
n(n+1)—2-3= (m—2)(n +3) 

a (n + 2)! 

D, {P,] = (n—1)n(n +1) P, 
(n — 2)! 
D, [P,'] = (n — 2) (n— 1) n(n +1) (n + 2) (n + 8) P, = 


und wir können G, (z), &; (z) usw. analog wie ©, (z) bestimmen. Um dazu die 
Operation D, [P,,] auszuführen, setzen wir: : ; 


Da 


Al: 


Zu 


un 
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D: 


un 


di 


i - : - 
~ 
da 
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(27) 
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(25) 


or D, 


die 
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dann ist: 
(22—1)g, = 229,4+ 92 


= 


(28) 
+9- 


Damit ergibt sich: : 
d ” 7 
{®-Da}-2g= —29,} = gf’. 
Also ist: 


D, [P„]= le 1)? 


Wir behaupten, es gilt allgemein: 


de” 
Zum Beweis setzen wir: | 
dann ist: 


1) gm- 1 = 2(m—1)Z2G9m-1 + Ym- 


Daraus folgt durch (m — 1) maliges Differenzieren die zu (28) analoge Gleichung: ia 


(30) 


(m-2) (m-1) 


= (m—1) + 


und schließlich durch nochmaliges Differenzieren : 
d d 
= m + 9. 


Nehmen wir nun an, daß (29) für m — 1 bewiesen ist, d.h. daß: 
Dm-1 [Pa] = 9m-ı 
TAR, 


ist, so folgt wegen (27) aus (30): 


Da (29) fiir m = 1 und m = 2 bewiesen ist, gilt (29) fiir beliebiges m. 
Wir erhalten also: 


11) Diese Gleichung kénnen wir wegen (29) auch in der Form schreiben: 
de” 
und aus einer Beziehung der zugeordneten Kugelfunktionen: 
m (n +m)! _ m 
direkt ableiten (s. W. Magnus u. F. Oberhettinger, Formeln u. Sätze, S. 75). 


| 
(26) | 
als: 
(27) 
= 
3 
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co 
b) Die Darstellung von R (r,y) mit Hilfe der Potenzreihe R (co, y) = u B, 2 


Wir gehen von der Reihenentwicklung (10) für R (co, y) aus und bekommen unter 
Anwendung von (32) die folgenden Gleichungen: 


se. 


;2 

2!2-4(xr)? 112-4-6(z2r® 


1779 | 
-4-6-8(xr)if 


tay" 


Führen wir die Differentiation aus, so können wir unter Benutzung von (6)  (r, 7) 
in folgender Weise iibersichtlich darstellen: 


(ix n| 


3-4 6 


Für den allgemeinen Beweis von (33) begnügen wir, uns damit, die Überein- 
stimmung an dem zu B, gehörenden Faktor zu zeigen, da sich daraus unschwer 


Sy (2) = B, + B,2* + B, +4 
~ 2! 6! , 
4! ‘ m , 6! m, 8! 
S; (2) = B, [z(2*— + B, [28 — + BI B, [2° — 1)?” + 
Dann folgt aus (21): A 


ein 


‘, H. 
die 
na 
4! 
B B, 
Pol 
ıB 
Yon 
. erse 
aus 
nacl 
Pa + 
t 1 
+5 
(3 
3 
m 
= (2 
(33) 
19) 
; 
7 
+5 
¥ ® Aı 


„ın 


inter 


(r, y) 


24 — z 

312 (xr) 2xr Dur 

[22 (22 — 1)27’’ ‚4-32. 3-4 
1%] _;4:32 3 +, 4! 

212-4 (xr)? 2-4(xr)® 2xr 2xr 2-4 (xr)® 
112-4-6(xr)® -4-6(xr)? 2xr 2° 2 2xr 
— 14) 8! 6! 4! 
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die allgemeine Richtigkeit erkennen läßt. Wir schreiben dazu die einzelnen Glieder 
nach Ausführung der Differentiation in geeigneter Form untereinander: 


4! 


Dabei ergibt sich durch Addition der einzelnen Kolonnen die behauptete Über- 
einstimmung. 

Wir haben damit das allgemeine Ergebnis: 

Die durch (7) definierte Funktion R (r, y) läßt sich unmittelbar aus der Potenz- 
reihe von ©, (z) ableiten, indem jede Potenz z* durch ein von x r und z abhängiges 
Polynom: 


1) 2 
Yon (z) = 2?" fenlixr)—i ~ fon-1 (34) 
(2n — 3) (2m —1) 22" (2 n)! 
ersetzt wı 


Gehen wir von der Kolbenmembran aus und benutzen dementsprechend die 
aus (13) für p = 1 sich ergebende Reihe für ©, (2): 


so können wir nach Potenzen von z (z= cosy) ordnen und beeen aus (2) 
nach Einführung der Funktionen S,, und C,, aus (5): 


Pati Pm (za) [Sp (xr) — Cy + J, (xa) [S, Ger) — 60, 
(za) [S, (xr) (xr)]+- 
112 \r 
7 (xa) [8,(zr) —ıC, war} 
cost 7 


32 


J, (ea) [Sp (en) Js ea) [83 (er) 


J; (x a) [S, (xr) —iC, (2) Gea) [8, er) iC, 
+35 (4) J, (xa) [S, (xr) —iC,(xr)]+-- 


_ (xa cos’ y {a 
93 Ir 


J, (za) [S, (xr) —iC,(xr)] (xr) (xr)] 


+3 
+: 
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Setzen wir cos y = 0 (d.h. der Aufpunkt liegt in der Membranebene), so ver- 
schwinden in (35) alle Glieder außer der ersten Zeile. Dies Ergebnis steht in 
Übereinstimmung mit einer früher'?) vom Verfasser abgeleiteten Formel. 


Prinzipiell ist die Entwicklung von R(r,y) durch (34) für beliebiges w (o) 
gegeben. Dabei kann w(o) z.B. als Ergebnis von Messungen durch eine Kurve 
bestimmt sein. Man hat dann nur nötig, die Koeffizienten B,, nach (11) durch 
graphische Integration zu berechnen. 


c) Die Darstellung der Entwicklungskoeffizienten durch Besselsche Funktionen 
Wir en nm uns nun die Aufgabe für w ee aus Sl?! (2) 
J, (xa y1- - 2?) 


(xa yı — 22)? 
führen der Differentiationen nach (32) zu berechnen. Wir setzen zur Abkürzung: 


= 2? p! die Entwicklungskoeffizienten © ol! z) durch direktes Aus- 


‚ J 
u=xa Vi-2, 9,=7 

= — (xa)? — und 
ei sab es folgt nach den Regeln der Differentialrechnung aus (32): ay ia 

3 
( -1) M\, 2m) 
Die höheren Ableitungen von p, ergeben sich aus: 
, 
4 
indem man beide Seiten wiederholt differenziert: er be, 2 fos 
9» 
(xa)? + (xa z) Pp+s 4 
a. 
15 4A 2G p43 + 10 (za 2)’ + (442) - 
12) H. Stenzel, Ann. Physik 41, 250 (1942). > 
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(36) 


So erhalten wir allgemein : eo 


xa 

x az)” 

2 
=1 (Zn 2 1) [a Z)Pp+n+i 2n Pp+n+2 
6 (38 


Setzen wir (37) und (38) in (36) ein, so kénnen wir nach Ausfiihrung der Dif- 
ferentiationen [(1 — 22)™\™,[ (1 — 22)” (m-2) usw. die Entwicklungskoeffizienten 


Gm (2) berechnen. So ergibt eine einfache Rechnung für die 4 ersten Koeffi- 
zienten die folgenden Formeln, wobei die Abkürzungen: u=xasiny und 


I»] 
= (u); 
o °*p: 
el» 
+ (6 — J,,,(u) 4+ hut? (u); 
(— 54 + 432 1? — 144 J,,, (u) 
+ (15 — 378 1? + 504 1 — 48 1%) u? Js,» (u) 
+ (45 2 — 198 1° + 720) wt-? Ju, ,(u) + (15 — 18 8) (u 
+ (—720 + 13824 1? — 18432 t + 2304 1%) w-? J,,, (u) 


+ (105 = 5520 12 + 18864 0 — 8832 1° + 3841) 
+ (420 #2 — 4080 t* + 4896 1% — 786 1?) u? Js,» (u) 
+ (210 — 720 + 288 #8) J,» (u) + (281° — 32 u’? (v) 
+ BuS-P Js, (u). 

Für z= 1, d.h. wenn der Aufpunkt auf der Mittelachse liegt, versagen die 


Formeln (39). Für z = 1 finden wir aber ze) — (fürn <m, und ZU” = 2-4 :6--- 
2m. Damit ergibt sich dann aus (36): Va A 
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Und wir erhalten wegen ¢,,,._ „= —— aus (37) und (38) fiir den Aufpunkt auf 
m: 


2 


3 der Mittelachse 4 
he @ n)!(2n)! | (xa)?” ri n (za)?"*? n(n— 1) (zuj?"+* 

n! p+n +1!" (41) 

( )= Dit 31(p Fn 

er: (p +n + 3)! 


Wir hätten (41) und (42) auch unmittelbar aus (20) unter Berücksichtigung von 
== (6) ableiten können. 


4. Die Darstellung von R(r,y) mit Hilfe von R(r,o) 


nach steigenden Potenzen von sin y 
a) Die Darstellung der Entwicklungskoeffizienten durch einen Differentialoperator 
ios Im 3. Abschnitt waren wir mit Hilfe der nach fallenden Potenzen von xr 


fortschreitenden Reihe für /, (i x7) von einer entsprechenden Reibe fiir R (r, y) 

- ausgegangen und konnten die auftretenden Entwicklungskoeffizienten ©,, mit 
Hilfe der Differentialgleichung der Kugelfunktionen bestimmer. Jetzt werden wir 

4 ganz analog mit Hilfe der nach steigenden Poztenzen von sin y fortschreitenden 
£ Reihe für die Kugelfunktionen von einer entsprechenden Reihe für R (r,y) aus- 
, gehen und die auftretenden Entwicklungskoeffizienten &,, mit Hilfe einer 
Differentialgleichung für f, (¢x7) bestimmen. Wir benutzen die Darstellung der 
Kugelfunktionen durch eine hypergeometrische Funktion"): 


1; sin? y). 


Daraus folgt: 


2n(2n +1) (2m —2)2n(2n + 1) (2n +3) 


P;a (cosy) = 1— e+ sint y — ++ +1), 
(43) . 
Setzen wir dies in (7) ein, so können wir R (r, y) in der Form schreiben: 
| 
R (r, y) = Xo (r) ER z, (r) T, (r) . (44) 


13) E.Hobson, Spherical and Ellipsoidal Harmonics, Cambridge 1931, S. 219. Dabei 
ist in der dort angegeben Formel m = 0 gesetzt. 

14) An Stelle (43) hätten wir auch die bekanntere Darstellung: 

2 ai 9 9 9 
2n(2n + 1) sin? »/2 + (2 —1)2n(2 hs 1) (2 n + 2) 
1? 1?- 

(s. Rayleigh, Theory of Sound § 336) benutzen können. Doch werden baad dann sich er- 
gebenden Koeffizienten T,, weniger übersichtlich. Ey 


Pn (cos y)= 1 — —— sin‘ y/2 - 


Da 


erse 


Es 


fizi 


die 


Dar 


Ord 


‘ 
= 
H.; 
= 
= 
ty 
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bac) 


auf 


41) 


42) 
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Xo (r) = Agy (x a) fan (a xr) = = R (r, 0) 
(r) = 2n(2n + 1) Agy (za) fon (Üxr) 


(r) = = (2m — 2) 2n(2n + 1) (2n + 3) Agy (za) fen (6% 7) 


usw. 


Wir wollen dann zur Vereinfachung auf der rechten Seite von (6) 7/2 x r durch g 
ersetzen und die entsprechende Funktion betrachten: 


Es macht dann keine Schwierigkeit, die beiden folgenden Behauptungen zu veri- 
fizie ren : 


1. f, (g) genügt der Differentialgleichung: ole 
6) 

2. Es besteht die Rekursionsformel: 

Ähnlich wie früher definieren wir nun die Operatoren: X ER 


= 


Dn fnl= D1 [Dm-1 m — 2) (2m — 1)® m-1 In]. (50) 
Damit ergeben sich aus: 


n(n + 1) — (2m — 2) (2m — 1) = (n — (2 m — 2)) (n + 2m— 1) 
die Beziehungen : 


DU + VW fy usw. 
Dann ist: 


(r) = [To], T(r) = [To], usw. (51) 


Eine direkte Darstellung von ®,, [T,] durch Differentialquotienten höherer 
Ordnung in Analogie zu der Formel (29) scheint hier nicht zu existieren. en 5 


| =a 
3 
ator 
xr 
A y) aa 
mit 
wir a 
den ¥ 
us- 
ner 
der 
= 
43) 
| 44 ) 
: | 
| | 


D,. fn] für m = 1, 2,... zu untersuchen. 
Um ®, [g” /„] zu finden, benutzen wir die Rekursionsformel (47). 
Damit folgt wegen (46) und (48): 


2 (2m + + 1) (m+ 2) nf 
=(4n+2 + (m —) (farı 
Und unter Berücksichtigung von (47) ergibt sich: 


hatm- 


2, [g" ha = 
Dann folgt aus (49) und (54): 


Und wegen (52): E 


Dz fn] = (rn — 


So ergibt sich schließlich allgemein: 


Din [q” 09 (n<1 


und (55): 


DO, (9 fal = + (m—1) fy. 


Analog leiten wir die weiteren Beziehungen ab: „ig: 
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b) Die Darstellung von R (r,y) mit Hilfe der 


n 


~ J.) 


Dr fnl=3¢ Inst + (n — 3) — 2) 


fn] =m (nm —1) 1) fram 


mn). 


(r)=1+ie, +- 7? & +1 +: 
+14-5e +: 
(i xr) fg (i xr) 


= 1-2 Quer 


xr)? 


Um den Operator D,, allgemein (7, 0) = anwenden zu können, 
ist es in Rücksicht auf die Reihenentwicklung (18) zweckmäßig, die Operation 


(55) 


Liegt dann ® (7,0) in der durch (18) angegebenen Form vor, so folgt aus (51) 


(56) 


(p 


_ fiir 


dab 


>16 H. 
all 
21 
(x 
ay gi 
| | 
fl 
(x 
= 
I 
Fü 
aus 
4 . ry ie 
wy 
+ 
E (52) 
(53) 
ae 
- 
A 
| 
+ 
. 
= 
— | 
‚(ixr) 


(52) 
(53) 
(54) 


(56) 


1 x 2 (za)? fs (i #7) 3(za)t f, (ixr) 
2 +e 31 (ure ++ 
G(r) 1-2 . 2-3 (xa)? fg(ixr) , .. B-4(xa)* f, (i xr) 
(57) 
(ry 2-3- (xa)? hier) 4-5(za)t f,(ixr) 
fo (ver) + in 5! 6! (2 xr)? 


Führen wir noch die Funktionen S,, und C,, entsprechend (5) ein, so erhalten wir 
aus (57), (2) und (44): 


Pa +i Pm = (ver) — Cy — 


a)* 


(xr)] 


[S, (xr) — C, 


(x a)® 
23-31 (xr) 


[S, (x r) (x r)] — 


(ka? 
1!1!2? (2-2! (xr) 


[S, (xr) 


(z a sin y)* {1 (za) 2-3 (za) 
3-4(xa)s 
+ er) — 
313196 15: S, (xr) Cg (xr)] — [S, (xr) (zr)] 


3-4-5 y 


+... 


= 5. Konvergenzbetrachtungen und Berechnung praktischer Beispiele 


Für das Folgende beschränken wir uns auf den Fall der Kolbenmembran 
(p= 1). Zunächst wollen wir (35) in etwas anderer Form schreiben, indem wir 
_ für die Besselfunktionen J, (xa) die Darstellung benutzen: 


a)" 
J, 
dabei ist: 


(x a/2)? (x a/2)* 


1!n 2) 


1 
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| und insbesondere nach (20) für die Kolbenmebran (p = 1 og 
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Dann folgt aus (35): 


(za)? 8, (x r) —iC,(xr) (za)! S,(zr)—iC,(xr) 
Pa ti Pm = 5 “r Xı 92.91 (x r)? Xq 


(x a)® S, (x r) —iC, 1/ -3/ 4 
4“3 


S, (xr) (x7) 2 (xa)* S, (xr) (xr) , 
1111282 | 2- Oe + 92.31 (xr 
A 
3(xa)* 8 r)—iC, (xr 
(xa cos {1-2 (xa)? S, (xr) (xr) -3(xa)*S; (ixr) —iCs (ixr) 


(xa cosy)® (1-2 -3(xa)? S,(zr) (Kr) S, (x r) — iC, (x (x r) r) 


313128 | 2-4! er += jp (x — 
3-4-5 (xa)® S, (xr) — iC, (xr) A 
usw. 


Wir sehen, daß die einzelnen Glieder der Sinus-Entwicklung (35) sehr ähnlich 
den entsprechenden Gliedern der Cosinus-Entwicklung (59) sind. Da die in (59) 
hinzutretenden Faktoren «,, 4/2 X; usw. unter einer festen (von n un- 
abhängigen) Schranke liegen, genügt es die Konvergenz von (58) zu beweisen. 

Daß der reelle Teil von (58) absolut für alle xa und xr konvergiert, folgt ohne 

reiteres aus s der Ungleic hung: 


Für den imaginären Teil benutzen wir die Ungleichung: 


(x r)” \ 


und betrachten in (58) etwa die unendliche Reihe in der geschw eiften Klammer der“ 
4. Zeile mit dem m-ten Glied: 


_m(m + 1) (m + 2) (x a)?” 
an = m Om+s m" 
2 (m + 3)! (x r) 


Unter FREIEN von (60) bekommen wir dann die Vergleichsreihe mit dem all- 
gemeinen Glied | bp, | > | am | . Dabei ist: 
b m(m + 1)(m + 2)1-°3-5---2m+9 (x a)?” 
Dan 2 2” (m + 3)! (x r) (x r)” 


15) Die Ungleichungen ergeben sich aus der Reihenentwicklung für S, und (,. 
Vgl. E.N.T. 12, 19 (1935). 
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Hier strebt der Quotient zweier aufeinanderfolgender Glieder 


Om (m + 3) (2m + 11) a? a? 
On 2m (m + 4) r? 


Die Vergleichsreihe ist daher für a/r < 1 konvergent und beschränkt. Daraus er- 
gibt sich, daß der imaginäre Teil von (35) für r >a konvergent ist. 

Wir wollen nun die Wirksamkeit der 3 Entwicklungen (58), (35) und (21) (39) 
für den Fall xa = 10, “T= 100 untersuchen. Mit Hilfe von (58) ee wir - 


Tabelle 1. 
Tabelle 1 Tabelle 2 
xr= 100 7 


xr— 100 
sin y | Pa +1 Pm 


Pa +4 Dm 


0 | —0,13716—0,47416 i —0,00330—0,00309 i 1 

0,1 —0,12700—0,41606 i —0,00396—0,00422 i 3 

0,2 —0,09846—0,26925 i 2 | —0,00588—0,00752 i 
0,3 | —0,05821—0,10029 i —0,00878—0,01364 i 
0,4 —0,01794+0,02294 i | —0,01180—0,02156 
0,5 +0,00912+0,06699 i | —0,01293—0,02881 i 
0,6 +0,06542+0,04500 i | —0,00893—0,02806i | 10 
0,7 +0,0047 —0,0003 i +0,0027 —0,0058 i | 10 
08 —003 i 


Dabei geben die Ziffern in der letzten Kolonne an, wieviel Zeilen in der Formel 
erforderlich waren, um 5 Stellen hinter dem Komma-zu berechnen. Die Ziffer 10 
bedeutet also, daß die Reihe (58) bis (k a sin y)"® ausgewertet wurde. Für siny > 0,7 
reicht auch dies nicht mehr aus, die Konvergenz wird dann immer schlechter. Es 
ist dann besser, die Cosinus-Entwicklung nach (35) zu verwenden. Hier erhalten 
wir die Tabelle 2. 

Wie wir aus der oben bemerkten Ähnlichkeit der Entwicklung von (58) und 
(35) erwarten konnten, ergibt sich hier eine gute Konvergenz so lange cosy S 0,7 ist. 

Schließlich berechnen wir p, + ? p,, nach der r-Entwicklung (21) (39) für » = 1 
und erhalten die Tabelle 3. 

Zunächst überzeugen wir uns, daß der Wert für y = 0, den wir direkt durch: 


Pa +4 Pm = — — 0,13716 — 0,47415 


berechnen können, bis auf 5 Stellen genau wiedergegeben wird. Das war nach dem 
Verlauf der einzelnen Summanden Go, Sıyaxr usw. auch zu erwarten. 
Aus diesem Verlauf werden wir aber auch für sin y = 0,1, 0,2, 0,3 usw. folgern, 
daß sämtliche Werte von p, + i p,, in Tabelle 3 bis auf mindestens 4 Stellen genau 
wiedergegeben werden. Wir können dies für die Werte sin y = 0,1 bis sin y = 0,7 
aus dem Vergleich mit der Tabelle 1 beweisen und ferner aus Tabelle 2 den Wert 


cos y = 0,6. d. h. sin y = 0,8 zum Vergleich benutzen. an 
Für einen Vergleich mit der bekannten Annäherungsformel: Ate ye 


—ie*" (xa)? 2J, (xasiny) 
xr 2 (xasiny) +’ 


= 


welche aus (21) hervorgeht, wenn wir uns auf die erste Näherung ©, (2) beschränken, 

ist in den letzten Kolonnen von Tabelle 3 dieser Wert und | p¥ +9 | ange- 

a Wir sehen daraus, | pe +ips | einigermaßen g genaue Werte gibt. Eine 


= 
= 


4 
11% 
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Tabelle 


“xr 


©, (z)/2-4 (x r)2 8 Ss (z)/2-4-6(ser)* it S,(z)/2-4-6-8( 


—0,04416 
—0,03450 
—0,01269 
+0,00613 
+0,01179 
+0,00625 
—0,00104 
—0,00324 
—0,00120 
+0,00046 
—0,00002 


| 

—0,00646i | +0,00085 
—0,00472 i +0,00057 
—0,00107 i +0,00004 
+0,00152 ¢ —0,00025 
+0,00158 ¢ —0,00015 
+0,00025 i +0,00005 
—0,00057 i +0,00008 
—0,00036 i +0,00001 
+0,00007 i —0,00002 
+0,00008 7 0 

—0,00002 i 0 


Ausnahme machen die 


Werte fiir sin y = 0,4 und 
sin y = 0,7. Das ist ver- 


ständlich, wenn man be- 


denkt, daß hier die 


Werte kasiny = 4 und 


xasiny = 7 in unmittel- 
barer Nahe der Nullstellen 


von J, (xasiny) liegen, 
da die ersten beiden Null- 


stellen von J, (x) durch 
x, = 3,832 und x, = 7,016 


gegeben sind. Es ist klar, 
daß dann die zweite Nähe- 


rung ©, (z) von wesent- 


a? 


lichem Einfluß sein wird. 
Dagegen können wir 


nicht erwarten, daß der 
von y abhängige Phasen- 


wer 


gang vong*=arctg p* /p* 
den wabren Verhaltnissen 


gerecht wird, da der mity 
veränderliche Faktor nur 


positive oder negative 


- 
- 
- 


Werte annimmt, d. h. die 
Phase bleibt zwischen 


zwei Nullstellen von 
J, (zasiny) konstant und 


360° 0° 


320° 


springt an den Nullstellen 


Abb. 2. Phasenfunktionen 


y = are tg p,,/p, für xr = 20, 
xr= 40, xr= 60, xr= 80, 


xr— 100 und die entspre- 
chenden treppenförmigenersten 
Näherungen 


80° 


selbs 
hau 


A 
halte 


(x a) 


am 1 
Be 0 1 +0,25000 i 
0,1 0,88010 +0,20699 
0,2 0,57672 +0,10394 7 
0,3 0,22604 —0,00031 7 
0,4 |—0,03302| —0,05574 i 
0,5 |—0,13103 | —0,05240i 
0,6 |—0,09223|  —0,01845i 
0,7 |—0,00134| +0,00895 i 
0,8 | 0,05866|  +0,01296i 
09 | 0,05451 +0,00305 i 
1 0,00869 —0,00255 
Ver a 
= SBASHERE 
\t | | | | | | 
0,3 
=: 
0,6 
| 
| |: 
ler 
| ¥ 
20 Ts, und 
3 2w° der 3 


die 
ind 
’er- 
be- 
die 
und 
tel- 
llen 
ren, 
ull- 
rch 
016 
lar, 
he- 
nt- 


pre- 
sten 


= 100 


x 
Pa *Pm 


* * 
Par 


0,13716—0,47416 7 
—0,12703—0,41605 
—0,09846—0,26925 
—0(,05820—0, 10030 i 
—0,01794 +0,02293 i 
+0,00910 +0,06698 i 
+0,01539 +0,04499 i 
+0,00486—0,00020 i 
—0,00892—0,02806 7 
—0,01257—0,02449 i 
—0,00330—0,00309 7 


|! Pa t Dm | 
| 


0,494 
0,435 
0,287 
0,116 
0,029 
0,068 
0,048 
0,005 
0,029 
0,028 
0,005 


—0,25318—0,43116 i 
—0,22283—0,37946 7 
—0,14602—0,24866 i 
—0,05723—0,09746 7 
+0,00836 +0,01424 7 
+0,03318 +-0,05650 7 
+0,02335 +0,03977 i 
+0,00034 +0,00058 
—0,01485—0,02529 7 
—0,01380—0,02350 7 


0,500 


0,029 
0,027 
0,004 


selbst um 180°. Hierauf wurde schon in der erwähnten Arbeit von H. Back- 


haus hingewiesen (s. Abb. 2). 


Als zweites Beispiel berechnen wir 9, + 7 
halten die entsprechende Tabellen 1’ und 2’. 


Tabelle 1’ 


—0,00220—0,00375 i 


für za = 10, xr = 20 und er- 


m 


Tabelle 2’ 


20 xr= 20 
sin Y Pa + 0 Pm cos y Pa +t Pm 
0 1,34058 + 1,27412 7 1 0 0,00262—0,04220 i 
0,1 1,27929 + 1,07688 7 4 0,1 0,01206—0,04667 i 
0,2 1,08048 +0,57024 7 6 0,2 0,04244—0,05593 i 
0,3 | 0,72939—0,02240 7 8 0,3 0,09552—0,05610 i 
0,4 | 0,27585—0,43539 i 9 0,4 0,15872—0,02429 i 
0,5 | —0,13090—0,49733 7 10 0,5 0,18579 +-0,05382 7 
0,6 | —0,3018 —0,2502 7 10 0,6 0,09921 +0,13746 i 
0,7 | —0,16 +0,05 i 10 0,7 —0,1309 +0,0944 7 


Hier ergibt sich wieder, daß die Sinus-Entwicklung fiir sin y < 0,7 und die 
Cosinus-Entwicklung fiir cosy < 0,7 zweckmäßig ist. Dagegen können wir für 
die r-Entwicklung keine brauchbaren Ergebnisse (bei Anwendung von 4 Glieder) 
erwarten, da das letzte Glied nicht mehr klein gegen 1 wird. Allgemein ist für 
größeres x a das ausschlaggebende Glied in ©, durch die höchste Potenz, d.h. 


d 3 
urch 


(x a)?/2xr = 2,5 ist dieser Wert ~ '/;, so daß wir nicht einmal die erste Stelle 
nach dem Komma genau bekommen werden. 


Setzen wir dagegen voraus, daß (x a)?/2 x r < 1 ist, so wird der obige Ausdruck 
<"/;9 und wir werden p, +75, dann auf 3 Stellen hinter dem Komma richtig 
finden. Das wird auch durch das erste Beispiel bestätigt, wo (za)?/2xr = !/, 
ist. Demgemäß wurden als weitere Beispiele xa = 10, x r = 40 nach der Sinus- 
und Cosinus-Entwicklung und xa = 10, xr = 60 sowie xa = 10, xr = 80 nach 
der r-Entwicklung berechnet (s. Abb. 2 und die Tabellen 4, 5,6). 


(x a)® 1 


2-14 
Do gegeben. Für xa = 10 und xr= 20, 


L2xr 


1 


ig 
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Tabelle 4 
xr= 40 
sin y Pa +t Pm cos y Pa +t Pm 
0 0,25784 +-1,12563 7 1 0 0,01349+0,0010 i 1 
0,1 0,26367 +0,99005 7 5 0,1 0,01616 +0,00337 i 4 
0,2 0,26331 +0,64252 i 6°] 0,2 0,02383 +0,01118 ¢ 5 
0,3 0,22254 +0,23165 i 8 0,3 0,03483 +0,02578 i 7 
0,4 0,13098—0,07861 i 8 0,4 0,04446 40,04745 i 9 
0,5 0,02418—0,19272 i 9 0,5 0,04322 +0,07023 i 10 
0,6 —0,0387 —0,1303 7 10 0,6 | 0,0190 +0,0729 ¢ 10 
0,7 —0,03 —0,003 i 10 0,7 | —0,027 +0,014 i 10 
Tabelle 5 
xr— 60 
| 2 ~ 4s 
2xr 2-4(xr)? | 2-4-6(xr)? Pe 
0 4 0,416677 | — 0,12269 | —0,02990: 0,00653 |—0,5324 + 0,6027% 
0,1 0,88010 0,344997 | —0,09584 | —0,021877 0,00443 |—0,4568 + 0,5439: 
0,2 0,57672 0,173237 | —0,03525 | —0,00493 i 0,00032 |—0,2712 + 0,3873: 
0,3 0,22604 | —0,000527| + 0,01702 | +0,0070371} —0,00191 |—0,0664 + 0,1897% 
0,4 | —0,03302 | —0,092917| + 0,03275 | +0,007337| —0,00116 |-+ 0,0683 + 0,0206 
0,5 | —0,13103 | —0,08733: | + 0,01735 | +0,001187| +0,00037 |+0,0972—0,0680i 
0,6 | —0,09223 | —0,030757 | —0,00290 | —0,00265: | +0,00065 | + 0,0505—0,0665% 
0,7 | —0,00134 | + 0,014917} —0,00901 | —0,001667| +0,00005 |—0,0079—0,0115s 
0,8 0,05866 | + 0,02160: | —0,00334 | + 0,00032¢| —0,00018 |—0,0314 + 0,0382: 
0,9 0,05451 | +0,005087} + 0,00127 | + 0,000387} —0,000003 |—0,0185 + 0,0429; 
1,0 | 0,00869 | —0,00424% | —0,00005 | —0,000077 | —O, 000003 |+ 0,0012 + 0,0080% 
sin y So(z (2) Ss (2) S, (2) i 
o(2) 2xr 2-4(xr)? |2-4- 6(xr)? | Pet 
0 | 1 0,312507 | —0,06901 | —0,012617| + 0,00207 |_0,0008 0.1814 
0,1 0,88010 0,258747 | —0,05391 | —0,009237| + 0,00140 |—0 5313 — 0,0979 
0,2 0,57672 0,129927 | —0,01983 | —0,002087| +0,00010 |—0,3548—0,0410: 
0,3 0,22604 | —0,000397 | + 0,00957 | + 0,002977| —0,00060 |—0,1462 + 0,0146i 
0,4 | —0,03302 | —0,069687| + 0,01842 | + 0,003097| —-0,00037 |+ 0,0139 + 0,0403% 
0,5 | —0,13103 | —0,065507| + 0,00976 | + 0,00050: | + 0,00012 | + 0,0797 + 0,0320: 
0,6 | —0,09223 | —0,023067 | —0,00163 | —0,001127; + 0,00021 |+ 0,0598 + 0,0086 
0,7 | —0,00134 | + 0,01118: | —0,00507 | —0,00070:| +0,00002 |+ 0,0032 —0,0070: 
0,8 0,05866 | + 0,016207 | —0,00188 | + 0,00013? | —0,00006 |—0,0364—0,0062% 
0,9 0,05451 | + 0,003817} + 0,00071 | + 0,00016: 0 — 0,0346 + 0,0013 
1,0 0,00869 | —0,00318: | —0,00003 | —0,00003: | 0 —0,0052 + 0,0026: 


Wir wollen das allgemeine Ergebnis folgendermaßen zusammenfassen: 

Um den Schalldruck der kreisförmigen Kolbenmembran mit dem Radius a 
im Aufpunkt r, y nach Amplitude und Phase zu berechnen, wobei r > a vorausge- 
setzt wird, unterscheiden wir 3 Bereiche: 
la<r<azaA, ys 45° 
2a<r<aatA, y> 45° 
3. <r < 00, y beliebig 
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und haben im Bereich 1 die Sinus-Entwicklung (58), im Bereich 2 die Cosinus- 
Entwicklung (35) und im Bereich 3 die r-Entwicklung (21) anzuwenden. 

Die abgeleiteten Formeln zeigen nun, daß auch für Schwingungsformen w (9) = 
(1 — o?/a?)P-!, wenn p=2, 3,...gesetzt wird, das entsprechende Schallfeld 
ohne größere Schwierigkeit berechnet werden kann. Durch Superposition sind 
wir dann in der Lage, auch allgemeine Schwingungsformen der Gestalt '%): 

w (0) = + a, (1 — 9?/a*) + a, (1 — ++ a, (1 — 0*/a?)" 
der Berechnung zugrunde zu legen. 

Schließlich ist noch darauf hinzuweisen, daß die für die Anwendung des Ver- 
fahrens geforderte Symmetrie des Strahlungsfeldes auch fiir eine Strahlergruppe 
erfüllt ist, deren rotationssymmetrische Einzelstrahler (z. B. vertikale Dipole) 
vertikal iibereinander angeordnet sind. 

Als einfachstes Beipsiel betrachten wir eine Z A 
Gruppe von zwei akustischen Strahlern (klein 
zur Wellenlänge), die im Abstand a auf der Z-Achse 7 ae 
angeordnet sind (s. Abb. 3). Dann haben wir die 3 dig y > 

4 


Aufgabe, die Funktion: 7 5 4 | 
nach fallenden Potenzen von r zu entwickeln. In 5) ae 
der üblichen Weise würde man r/r, und r/r, aus: Pr 
j d d?\-"/s 
cosy + 
4 nr (63) Abb.3. Zur Berechnung des 
und (i „s Y =) Schallfeldes der Strahlergruppe 
2 


nach dem binomischen Lehrsatz entwickeln, müßte dann die Exponentialfunktionen 
durch ihre Reihenentwicklungen ersetzen, und die auftretenden Potenzen (r — r,)* 
bzw. (r— r,)” durch die sich aus (63) ergebenden Reihen ersetzen und mit den 
Reihen aus (63) multiplizieren. Das bedeutet eine ziemlich langwierige Rechnung. 


Gehen wir dagegen nach unserem Verfahren von: 
ad 
Hi (co, y) = cos Cosy) = cos az 


aus, so finden wir die Entwicklungskoeffizienten mit der Abkürzung = 


== a, 
cos y = z unmittelbar aus: 
(2) =—a [(@—]) sinaz] = a*sin*y cos (a cosy)—2acosysin(acosy) 
&, (z) = — a? [(2*— 1)? cosa = a* cos (x cos y) sint'y — 12 cos*y +4] 
— cos y sin* y sin (& cosy) 


16) Vel. H. Stenzel, Ann. Physik 7, 972 (1930) u. N. W. McLachlan, Ann. Physik 
15, 440 u. ff. (1932), wo die charakteristischen Funktionen R (00, . für r derartige Schwin- 
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Wir beschränken uns auf y = 90°, dann ergibt sich: &, (2) = 1, ©, (z) = at, 
©, (z) = a? (4 + a2), Sy (z) = at (18 + x?) und damit die Entwicklung: 


. nd? d? ad? . nd a? d? 
Rr, W)= 1 - ‘Gaia? + 


4ir 


Dieselbe Entwicklung ergibt sich noch einfacher, wenn wir von 


4 


(— i)" (2n)! 


ausgehen und nach (34) 2?” durch 
2” n! (x r)” 


f„ ersetzen: 


- ix 


Nach den bisher bekannten Formeln macht die Berechnung des Schallfeldes 
einer in starrer Wand schwingenden Membran bei stärkerer Bündelung fiir Auf- 
punkte in größerer Entfernung erhebliche Schwierigkeiten. Zur Behebung dieser 
Schwierigkeiten werden in der vorliegenden Arbeit 3 verschiedene Reihenent- 
wicklungen abgeleitet, wobei jede Reihe für einen bestimmten Bereich des Auf- 
punktes vorteilhaft ist. 
Das dazu benutzte Verfahren, welches die Entwicklungskoeffizienten der für 
den Schalldruck charakteristischen Funktion R (r,y) durch eine auf R (co, y) 
bzw. #(r,o) ausgeübte Differentialoperation zu berechnen gestattet, ist dabei 
von wesentlicher Bedeutung und dürfte sich auch für die Berechnung allgemeinerer 
_ Strahlungsfelder (akustischer oder elektrischer Art) als nützlich erweisen. 


a Wolfenbüttel, Post- und Fernmeldetechnisches Zentralamt 
des vereinigten Wirtschaftsgebiets. = 
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Über eine besondere Form der Zustandsgleichung 


Kü 
= Von C. Kämmereı 

( Mit 2 Abbildungen) 


Es wird eine Zustandsgleichung angegeben, die den kritischen Koeffizien- 

7 ten enthalt und daher fiir verschiedene Stoffe auch in dimensionsloser Form 

verschieden lautet. Dementsprechend ergeben sich aus i hr (auch bei dimensions- 

loser Auftragung) Dampfspannungskirven mit verschiedener Anfangsneigung, 

die der wirklichen nahekommt. 

Ides Die ursprüngliche Form!) der betreffenden Zustandsgleichung: 
Auf- a RT 
leser P+=e-b= (1) 
Auf- läßt sich durch Hinzufügen eines Exponenten 6 von v im Klammerausdruck, wie 
folgt, ergänzen: 
a 
9,7) 
abei Setzt man hierin (n—ß) = 2 und ab = h, so erhält man: Pe, 
lerer RT a h = 


p in kg/m? den äußeren Druck. nd ae 
v in m?/kg das spezifische oder kg Volumen. 
T die abs. Temperatur in °K. 
R die spezielle Gaskonstante in m/Grad. 
- a, h, n aus den kritischen Daten. und der Gaskonstante zu bestimmende 
Festwerte. 
p,. [kg/m?], v, [m?/kg] und 7’, [°K] die kritischen Daten. 


Zur Berechnung der 3 Konstanten a, A und » werden, wie üblich 3 von der 
kritischen Isotherme (7 = 7;,) im kritischen Punkt (p = p,, v = v,) zu erfüllende 
Bedingungen herangezogen und zwar: 


1. die Zustandsgleichung rel 
2. Die Beziehung: (2), = 0 (horizontale Tangente im kritischen Punkt) SEHR 


3. Die Bedingung: (=) = 0 (kritischer Punkt = Wendepunkt). es 


1) Akademie der Wissenschaften Wien, Akad. Anzeiger Nr. 11, Nov. 1943. Pr 27 
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(Wie später gezeigt wird, ist t die Einhaltung der 3. Bedingung nicht unbedingt 
erforderlich). 
Nach Einführung der dimensionslosen Kenngröße x = ET. (Kehrwert des 


kuitischen Koeffizienten K) erhält man die Festwerte von Gl. (3): 


a= ore 4 
a h = p, v;, (1 — 2 x)?/(4 2? — x) 


| n= 1/(1— 22) 


und nach Einführung der reduzierten oder dimensionslosen Drucke, Volumina 
und Temperaturen, p, = p/p,, v, = v/v, und T, = T/T, die dimensionslose Form: 


(1-22) 
P, = rv, (3a) 
Während sich aus bekannten Zustandsgleichungen nur ein einziger, für alle 
Stoffe gültiger Wert von x (z. B. ®/; = 0,375 bei van der Waals, und */,, = 0,267 


bei A. Wohl) ergibt, ist in Gl. (3a) die Kenngröße x mit ihrem wirklichen von 
Stoff zu Stoff verschiedenen Wert einzusetzen. Wegen a 


x — Pete Pate contd 
RT, 848-7, 
ist neben den kritischen Daten auch das Molekulargewicht m [kg/kmol] als Be- 
stimmungsgröße in Gl. (3a) enthalten. 

Auch bei dimensionsloser Auftragung ergeben sich sonach für verschiedene 
Stoffe verschieden verlaufende Isothermennetze, so daß das „Gesetz der 
übereinstimmenden Zustände“ weder für hohe, noch für niedrige Drücke zu 
Recht besteht. Für große Volumina v, geht os (3a) in die dimensionslose 


Gasgleichung: 


über. 


Gl. (3a) lautet z. B. für Wasserstoff aay’ mit x = 0,306 " = 3 BR 


_ 4, 2, 
= 3,27 - : 


oder für Kohlensäure (CO,) mit 2 = 0,276 le == 3 md eae 


9, ‚62 7 
Pz = 3,624 - 


Pr 
> 


rape 


Für Stoffe mit x < 0,25 wird (4 2—1) negativ und die Glieder 2 und 3 der 


rechten Seite von Gl. (3a) vertauschen ihre Vorzeic hen. 
So lautet z.B. Gl. (3a): 


für Wasser mit x = 0,225 (n — 1,82) 


13,46 


1,52 


2. = 4,06 v v vi 
oder fiir Ammoniak (NH,) mit x = 0,243 (n = 1,945) 
35,7 38,82 
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Die Kenngröße x bewegt sich etwa in den Grenzen von 0,225 (H,O) bis 0,327 
(Helium). Der nach van der Waals sich ergebende Wert von 0,375 liegt also zu 
hoch, woraus sich die starken Abweichungen dieser Gleichung besonders für 
Wasser erklären lassen. 

Gl. (3a) stimmt zwar besser mit der Erfahrung überein; doch ist auch bei ıhr 
die Lage des Isothermennetzes z. B. für Wasser nicht befriedigend. So verläuft 
z.B. die kritische Isotherme für H,O und höhere Drucke erheblich unterhalb der 
wirklichen. 

Um die Abweichung zu beheben, kann man im kritischen Punkt eine Unstetig- 
keit der 2. Ableitung zulassen, so daß nur 2 Bedingungen, nämlich: 


1. die Zustandsgleichung 


9 . ep 
2. die Bedingung: (2 “ 
Exponent n von v fiir v, > 1 und v, < 1 verschieden (jedoch > 2) gewählt werden 
kann. Die Waagerechte durch den kritischen Punkt bleibt trotzdem eine Wende- 


= 0 (horizontale Tangente) zu erfüllen sind und der 


tangente. 
seht man wieder von Gl. (3) aus, so erhält man die dimensionslose Zustands- 
gleichung: 
T, - (3b) 
Ps rv, x(n — 2) x(n — 2) ) 


Mit n = 2?/, = 2,666 für v,, > 1 und n = 21), = 2,1666 für v 


„< 1 (Flüssigkeit) 
ergibt sich z. B. für H,O mit x = 0,225: 


fiirv, >1: p, = 4,45-—-— (30) 

T, 18,12 14,¢ 
2 — + (31) 


2,666 für v, > 1: Ä 


pP, = 3,625 - 


r z 


Die Übereinstimmung der so berechneten Drucke mit den VDI-Wasserdampf- 
tafeln 1937 für H,O und einem 7'vS-Diagramm für CO, (Os tert a ER ‚prozesse, 
II. Auflage 1933) ist aus Zahlentafel 1 ersichtlich. 


. 


Kritische Isotherme von H,O (T, 647,4° K, r, 0,60304 m?/kg) x = 0,225 
vm?/kg 0,056145 | 0,02521| 0,01434 0,00304 | 0,002209| 0,00182 
Vz 18,46 8,3 4,725 1 | 0,727 0,6 

Pp, berechnet 0,22166 0,4461 0,6828 1 | 1,125 1,32 

p berechnet [ata | 49,8 100,2 153,5 225,06 253 1297 

p nach Tafel [ata | 50,0 100 150 225,6 250 {300 


Kritische lsotherme von CO, (7), = 304° K, v; = 0,002156) x = 0,276 


vm/kge 0,03 0,02 0,015 0,008 0,005 0,002156 
=0/v, 13,93 9,29 6,97 3,71 2,308 1 
p, berechnet 0.2362 | 0.33703 | 0.4295 | 0.68445 0878| 1 
p berechnet [ata] 17,75 (25,3 32,2 51,3 660 | 7% 
p nach Tafel [ata | 18,00 26,0 33,0 5,2 660 | 


Ann. Physik. 6. Folge, Bd. 4 
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Das Ieitigche Volumen v, wird allerdings verschieden angegeben. So liegt dem 
für H,O in Zahlentafel 1 benutzten Wert von x = 0,225 ein krit. Volumen v, = 
0,00304 zugrunde, (VDI-Wasserdampftafel 1937). Mit v, = 0,0031 (Hütte, 
26. Auflage) wird hingegen x = 0,2294 und für v, > 1: 

6,8 3,54 
p, = 4,36 — (32) 

Nach dieser Gleichung erhält man z. B. für die kritische Isotherme des Wasser- 

dampfes und v, = 4,725: 


p,, = 0,6615 oder 150 ata, also genaue Ubereinstimmung mit dem Dampftafel- 
wert. . 

Bei Assoziierung nimmt das (scheinbare) Molekulargewicht m zu; gleichzeitig 
ändern sich aber auch die kritischen Daten, so daß sich zunächst nicht 
angeben läßt, ob durch die Assoziierung auch die Kenngröße x = u bl 


848 T, 
einflußt wird oder nicht. 


Wahrscheinlich ist in Gl. (3a) und (3b) auch das 2. und 3. Glied rechts tem- 
peraturabhängig. Vgl. die Zustandsgleichungen von Berthelot, A. Wohl, 
Beattie u. Bridgman. Diese Abhängigkeit ist jedoch keineswegs einfach und 
bildet ein Grundproblem der Zustandsgleichung. Dividiert man zunächst das 
2. und 3. Glied rechts in Gl. (3b) durch eine Temperaturfunktion f, = f(T), 
ur. für T,—= 1 ebenfalls = 1 wird), so ergibt sich: 


ass T 1 /fn—1—xn 1 


Nun muß bekanntlich im 9,7 ,-Diagramm bei v, = 1 und 7’, = 1 die Tangente 
an die Dampfspannungskurve gleichzeitig auch die v, = konstant-Linie oder Iso- 
volume (v, = 1) tangieren. Der aus der Zustandsgleichung sich ergebende Wert 
des Differentialquotienten: 


(a ) ge =1 tgx muß also gleich a 


sein, wenn mit p,, und 7',, die dimensionslosen oder reduzierten Sattdampfdrucke 
und -Temperaturen bezeichnet werden. 

Die Anfangsneigung (tgx) der Dampfdruckkurve im p, T,-Diagramm läßt 
sich daher durch Differenzieren aus der Zustandsgleichung berechnen. 

Setzt man nun in Gl. (3c) für den kritischen Punkt als erste Näherung: 


hie fer 
eT, = fa => Pr x(n —2) x(n —2) 

da (Feat „= List. Auf tg « hat also in diesem Fall die Wahl des Exponenten » 


keinen Einfluß. (K = kritischer Koeffizient = 1/x). 
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Die wirkliche Anfangsneigung der Dampfdruckkurve gegen die Abszissen- 
achse kann man andererseits aus einem Schaubild Abb. 1 (nach W. Schüle, 
Technische Thermodynamik, II. Auflage, II. Bd. 8. 80 Fig. 34) entnehmen, in 
welchem log 1/p,,, über (1/7',,—1) als Abszisse aufgetragen ist. 


Nach der bekannten 4 


Dampfdruckformel: 
N 
Fed gm P 
erhält man sonach, da M für 
Festwert bleibt, durch den, _ 
Ursprung gehende Gerade oder _ W4 


Strahlen, deren Neigungs- 


winkel gegen die Abszissen- HT gs 


achse durch die Beziehung: 7 7, Porn ber 
— M gegeben ist. Nach noch 
Gl. (5) ist: 05 G A 
1 9 © 
und Abb. 11. Dampfdruckkurven nach Schüle 
dp, ‘ 9£ ‘ 9£ 
( = tea = 2,3026. 2,3026 tg 6. 
In Zahlentafel 2 sind die nach Gl. (4) aus der Zustandsgleichung berechneten ee 


Werte von tg x den aus Abb. 1 entnommenen gegeniibergestellt . 


Zahlentafel 2 zu Abbildung 2 


| | 0,00304 | 0,0031 | 
co, | NH, | | H,O | C,H,OH 
—1_ | 976 943 | 99% 0,22% | 
x K~ 848 0,306 0,276 0,243 | 0,225 94 0,244 2 
K=1/x 3,27 3,620 4,12 4,45 35 
Atomzahl 2 3 4 | 3 % ae 
tga berechn. nach Gl.(4)} 5,54 6,248 7,235 | 7,9 | | 
tgaaus Abb. 1 6,11 6,4 6,91 | 7,4 
Für van der Waalssche Stoffe mit x = 0,375 wäre tga nach Abb. 1 = 3,86, 4 
während sich nach Gl. (4) ein höherer Wert, nämlich 4,299 ergeben würde. 
Die errechneten Werte von tgx stimmen mit den gemessenen nur annähernd ® 


überein ; doch ist der Gang beider Reihen der gleiche (siehe Abb. 2) und der Vorteil 
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nen, den die Einführung der wirklichen von Stoff zu Stoff verschiedenen 
Kennzahlen x in die Zustandsgleichung mit sich bringt. 

Aus Abb. 2 geht ferner hervor, daß die Neigung der Dampfdruckkurve zwar 
im wesentlichen von der Kennzahl x, jedoch außerdem noch von einem zweiten 
Parameter (wie etwa der 


Be Atomzahl) abhä B 
OH omzahl) abhängen muß, 
da 2 Stoffe mit gleichem z 
ey (Alkohol und Ammoniak) 
, verschieden geneigte Dampf- 
8} 

Ben ) Werte von tga Natiirlich kann man 

nach Abb. 1 

© nach GL. (4) berech- durch Verwendung anderer 
7} Temperaturfunktionen (z. B. 
T:,a + usw.) auch ge- 
Why naue Ubereinstimmung von 
Be tga (berechnet) mit tga 
st (gemessen) erzielen. In 
ae erster Linie sind aber f,, 
Reece und f,, so zu bestimmen, 
ae X= daß das resultierende Iso- 
thermennetz für 7, > 1 mit 
dem wirklichen möglichst 
x=225 genau zusammenfällt. f,, 
Boe und f,, miissen ferner wegen 
a | des Nernstschen Wärme- 
#0 ? 04 satzes bei 7’, = Oeine waage- 
GHsOH Ap rich Waals rechte Tangente haben. Da 
a NH { x=0375 der Sattdampfdruck für alle 
3 Stoffe bei einer gewissen 
Temperatur (z. B. der 
Anfangsneigung einiger Schmelztemperatur (T,= 
TR 0,422) bei H,O) verschwin- 
RR al det, d. h. konstant wird, 


? 


2 


ist anzunehmen, daß sich auch f,, und f,, unterhalb dieser Temperatur nicht 
mehr wesentlich ändern. Auf die weitere Bestimmung der Temperaturfunktionen 


in Gl. (3c) kann aber hier nicht eingegangen werden. 


Wien 14/89, Penzingerstr. 45.8 4 
(Bei der Redaktion am 17. Juni 1944.) 
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klang gebracht werden kann und auf die wichtigsten Folgerungen hingewiesen, 


& 


Von Erich Kretschmann 


Inhaltsiibersicht® 


Es wird gezeigt, wie die einfache Annahme streng BEER TONER Elementar- 
ladungen mit den Erhaltungssätzen von Energieimpuls und Drehimpuls in Ein- 


die sich aus der damit vorgeschlagenen Abänderung der Elektrönentheorie xg 
ergeben. 


Die klassische Elektronentheorie leidet ah der inneren Unstim migkeit, daß einer- 
seits die tatsächliche Unteilbarkeit aller Elementarladungen vorausgesetzt, 
andererseits aber jede dieser Ladungen auf ein endliches Gebiet verteilt gedacht 
wird und Wechselwirkungskräfte zwischen ihren Teilen als bestehend angenommen 
werden !). ; 
Die vorliegende Bemerkung soll kurz darauf hinweisen, daß man diesen Miß- 
stand und die mit ihm verbundenen Scheinprobleme ohne Verletzung der Erhal- 
tungssätze von Energie und Impuls durch eine einfache und naheliegende Änderung 
der bisher für diese Größen als gültig angehommene Ausdrücke beseitigen kann. Be: 
Es genügt, die Felder zweier Elementarladungen e,, h, und e,, h,, und das aus “ 


ihnen nach dem Superpositionsprinzip zusammengesetzte Gesamtfeld e = e, + ¢2, ‘a. 
h = §, + zu betrachten. 
Die Feldgleichungen sind linear in allen Feldkomponenten und gelten einzeln Sam 


fiir das Feld (1) wie fiir das Feld (2) wie fiir das Gesamtfeld bei beliebig vor- 
gegebener, etwa durch irgendwelche nicht elektromagnetische Kräfte bewirkter, 
Bewegung und Verteilung der Ladungen. - SS 
Aus den Feldgleichungen folgen die Erhaitungssätze von Energie und Impuls 
und gelten demnach ebenfalls bei beliebiger Bewegung und Verteilung der Ladungen 
einzeln fiir die erste wie fiir die zweite Ladung, jede mit ihrem Felde, wie fiir beide ony 
Ladungen zusammen mit dem Gesamtfelde. 
Die Erhaltungssätze haben alle die Form ce 


ot 


1) Zur Zeit der Entstehung der Elektronentheorie, nach der makroskopischen Elektro- 
dynamik und vor der Einsteinschen Relativitätstheorie und der Bohrschen Theorie 
der Atomspektren, mußte natürlich eine endliche Ausdehnung der Elementarladungen mit \ 
elektromagnetischer Rückwirkung auf sich selbst als unvermeidbar erse :heinen. Denn ~ 1 
nur so konnte scheinbar eine unendliche Eigenenergie der Elektronen vermieden und te 
zugleich die beobachtete Abhängigkeit ihrer trägen Masse von der Geschwindigkeit N 
erklärt werden. S. M. Abraham: Theorie der Elektrizität II, 3.Kap., $ 16, S. 139 (1905). 
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des Energiesatzes und jedes Glied in ihnen ist vom zweiten Grade in den Feld- 
komponenten und der neben diesen in den Feldgleichungen linear auftretenden 
Ladungs- und Stromdichte. So ist im Beispiel der Energiegleichung 


S= is [e,5, w= {ee + und g=(e,ev) mit 479 = dive. 

Jedes Glied der Gleichungen zerfällt infolge der Zweiteilung : 
e +e, H=h+h, 


in drei Teile, von denen der erste nur Feld und Ladungsdichte des ersten Teilchens, 
der zweite dieselben Daten des zweiten, und der dritte nach dem Schema (1) - (2) 
+ (2)-(1) die beider Teilchen gnthalt. 

Die Glieder der ersten zwei Teile erfüllen nach dem Gesagten die Erhaltungs- 
gleichungen für sich allein, und folglich tun die dritten Teilglieder, die Wechsel- 
wirkungsglieder, wie ich sie nennen will, dasselbe. Der Energiesatz gilt also auch 
in der Form: 


div| Bel + Lee. 6:3} {ler ee) +91, + (es, 0200) + 
(2) 
Dieser, gewiß nicht neue, Schluß läßt sich ohne weiteres von zwei auf n > 2- 
Elementarladungen übertragen. Denn streicht man aus den für das Ganze auf- 
gestellten Erhaltungsgleichungen alle aus Feid, Ladungs- und Stromdichte je 
eines Teilchens aufgebauten Glieder, die Eigenglieder, so erfüllen wie diese auch 
die zurückbleibenden, bei großem n viel zahlreicheren ,,Wechselwirkungsglieder“ 
die Erhaltungssätze für sich allein. Dabei treten als Dichte des Energiestromes, 
der Energie und der Arbeitsleistung der Wechselwirkung statt der Ausdrücke in 
(2) auf: ‘ 


> i,k [e; == > Lk q= > lk (C1; Ox (3) 


Ganz entsprechend sind in der Erhaltungsgleichung z. B. der z-Komponente des 
Impulses der Wechselwirkung 


> 


da ‘Spannungskomponenten und die Kraftdichte aus den bekannten Ausdriicken 
dieser Größen für ein einheitliches Feld zu entnehmen: 


1 1 


s während g = S/c* die Impulsdichte ist. 
Man kann die Doppelsummen hier leicht durch einfache Summen ersetzen, 
wenn man das „äußere‘‘ Gesamtfeld einführt, das auf die |-te Elementarladung 


wirkt: 
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das wäre also das Feld schlechthin, das man mit der /-ten Elementarladung als 
Probekörper mißt. Mit diesem wird: 


(er, Ho) = zz D> (Ew, bil. {er, Em) + (hr, Ho)} 


q= >! (Ey, a%), 12" En 4 & dy 6) 
1 
= {1 Eo. +5. Tey = >! fer, Eo, + bx 


Man sieht an diesen Ausdriicken besonders leicht, daB sie da, wo es auf das Feld 
ainer einzelnen Elementarladung nicht ankommt, oder 


En = & — >! e; und Da = D = >! hy 


gesetzt werden darf, in die altgewohnten Formen iibergehen. 

Aus der differentiellen Form (la) ds Impulssatzes folgt, da der Spannungs- 
tensor symmetrisch hist, der Drehimpulssatz in Integralform ?). 

Die angegebenen und hiermit zur Annahme vorgeschlagenen Ausdrücke (3) 
und (5) erfüllen also alle Erhaltungssätze von Energie, Impuls und Drehimpuls, 
und zwar in der Weise, daß alle Beiträge zu diesen Größen, die inneren Kräften 
zwischen gedachten Teilen je einer Elementarladung zugeschrieben werden müß- 
ten, aus dieser Bilanz der elektromagnetischen Energie und Bewegungsgröße ein- 
fach gestrichen sind. Es gibt nach dieser Auffassung keine Rückwirkung des 
elektrischen oder magnetischen Feldes einer Elementarladung auf diese selbst, 
keine elektromagnetische ,,Eigenenergie“ oder ‚eigene‘ Bewegungsgröße. 

Alles Gesagte gilt unabhängig von Gestalt und Größe der Elementarladungen. 
Für die Elektronentheorie als solche ist es bei Annahme der Gleichungen (3) oder 
(5) das einfachste und methodisch allein Richtige, die unteilbaren Elementar- 
ladungen als streng punktförmig zu behandeln und ihnen rein mechanische Massen 
zuzuschreiben ®). 

Dabei kann man für das Feld der einzelnen, /-ten, Elementarladung die gewohn- « 


ten Feldgleichungen beibehalten. 
\ 
1 ce 4a 
rot h, = = + — ad: dive, = 4270, 
(6) 
pte; = div h, = 0 
und braucht nur die Gleichung 
9, O(t—t)) ‚ (6a) 


ftir die Ladungsdichte der am Orte r, befindlichen Elementarladung hinzuzufügen. 
ö(t — r,) ist die bekannte, uneigentliche Dirac-Funktion *). Dieselben gen 
gelten dann natürlich auch für das Gesamtfeld E= Ye, § = Yh, mit 
0b = So,v, wobei » = v, an den Ladungsorten (0, + 0) ist. 


*) M. Abraham, Theorie der Elektrizität II, § 5, S. 34 (1905). 

3) Integriert man die Energiedichte u nach (5) über den von einer Elementarladung 
eingenommenen Raum, so erhält man eine vom äußeren Felde abhängige Größe, die für 
punktförmige Ladung gleich Null ist, wie man sofort sieht. 

*) Das heißt fdr (x — r,) ist gleich 1 oder 0, je nachdem das 
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ee Felder, die von keiner Elementarladung herrühren, sind hierbei ausgeschlossen 5), 
+ Zum Schluß sei ganz kurz auf zwei ungewohnte Folgerungen aus den neuen 
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_ Ausdrücken (3) und (5) hingeweisen. 

Erstens kann nach (3) die Energiedichte des Feldes ebensogut negativ wie 
positiv sein, statt definitiv positiv bisher. Das kommt daher, daß das Raumintegral 
der neuen Feldenergie nur die Arbeit darstellt, die nötig ist, um alle Elementar- 
ladungen an ihre Plätze und in ihren Bewegungszustand zu bringen, während in 
dem gewohnten Ausdruck der Feldenergie noch die zur Zusammensetzung der 
einzelnen Elementarladung aus infinitesimalen Teilen "erforderliche Arbeit ent- 
halten ist, die physikalisch unmöglich ist. 

Zweitens gehören nach (3) zur Aussendung von Strahlungsenergie stets min- 
destens zwei Elementarladungen, und zwar müssen beide beschleunigt sein. Denn 
ruht die eine, e,, im mitbewegten Bezugssystem dauernd, so fließt der Energie- 
strom, den ihr rein elektrisches und radiales Feld mit dem einer zweiten Ladung e, 
erzeugt, nach (3) nur längs der um e, beschriebenen Kugelflächen. Bei Bewegung 
von e, z. B. auf einer Ellipse um den festgehaltenen Brennpunkt e, würde also e, 
nur aus dem Gebiete zwischen den durch Perihel und Aphel gelegten Kugeln 
Energie abwechselnd aufnehmen, und wieder dahin abgeben, und zwar im Maße der 
von der Anziehungskraft bzw. gegen sie geleisteten Arbeit. 

Die Bewegung eines Elektrons auf einer Bohrschen Bahn ist hiernach bei fest- 
gehaltenem Kern strahlungsfrei ®). 

Überhaupt fällt nach (3) alle inkohärente Eigenstrahlung einzelner Elementar- 
ladungen fort. Dieses Erscheinungsgebiet ist damit von, vornherein aus der 
klassischen Elektronentheorie ausgeschlossen und der Quantentheorie zugewiesen. 

Dagegen bleibt die Arbeitsleistung bestehen, die das Wellenfeld jeder schwin- 
genden Elementarladung an jeder anderen leistet, die es in Bewegung setzt. Den 
Ausgleich für diese Arbeit bildet die (negative) Interferenzenergie der Felder 
beider Ladungen. Das ist auch in der gewohnten Darstellung nicht anders. Nur 
wird in dieser der Vorgang wieder verdeckt durch die Eigenenergien der Ladungen, 
die an ihm nicht mitwirken, aber durch ihre bloße Überlagerung über die Inter- 
ferenzenergie die Feldenergie überall positiv machen. 

Makroskopische Ladungen und Ströme behalten, wie schon bei Gl. (5) gesagt, 
praktisch ihre alten Feldenergiewerte. Im Geltungsbereich der klassischen Elektro- 
dynamik wird also durch die neuen Ausdrücke nichts geändert. 

5) Dies käme noch ausgeprägter zur Darstellung, wenn man als Grundgesetz die 
Lienard-Wiechertschen Gleichungen für die elektromagnetischen Potentiale der 
Elementarladungen und die Definitionen der Feldgrößen durch sie aufstellte. 

6) Auch nach der bisherigen Theorie strömt, wenn das Elektron auf einer elliptischen 
Spiralbahn den ruhenden Kern umkreist, die Wechselwirkungsenergie des Feldes gemäß 
3. abwechselnd in das Elektron und aus ihm heraus, nur jetzt überwiegend hinein, genau 
im Maße der überwiegend positiven Arbeitsleistung des Kernfeldes an dem Elektron. 
Zugleich leistet das hemmende eigene Feld an dem Elektron negative Arbeit und das 
Aquivalent dieser Arbeit wird als Eigenstrahlungsenergie ausgesandt. Dieser zweite 
Vorgang ist aber unabhängig von dem ersten. Denn wenn man mittels mechanischer 
Kraft das Elektron unbeschleunigt dem Kerne nähert, so fällt die Strahlung fort, und 


es bleibt nur das Einströmen von Feldenergie in das Elektron gemäß 3. übrig nebst ihrer 
Umwandlung in mechanische Arbeit. 


bn Halle (Saale), Institut fiir theoretische Physik. iu 


(Bei der Redaktion eingegangen am 10. Dezember 1948.) 
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Me Von F. Möglich und R. Rompe RER 
(Mit 2 Abbildungen) 


Inhaltsübersicht 


Es wird ein Zusammenhang aufgezeigt zwischen der Eindringtiefe magneti- 
scher Felder bei der Supraleitung und der Grenzwellenlänge eines Absorptions- 
gebietes des betreffenden Metalles im kurzwelligen Teil des Sichtbaren. Ferner 
wird gezeigt, daß die Stärke der Absorptior dieses Gebietes von der linearen Aus- 
dehnung der Metalle sehr stark abhängt, und zwar mit geringer werdender Aus- 
dehnung zunimmt. Am experimentellen Material, bei Silber, Quecksilber und 
anderen Metallen wird der. Nachweis geführt, daß ein solches Absorptionsgebiet 
wirklich existiert und die Eigenschaft der Abhängigkeit von der Schichtdicke 
usw. besitzt. Mit Hilfe der neuesten Messungen über die Eindringtiefe magnetischer 
Felder in Supraleiter wird gezeigt, daß dieses Absorptionsgebiet gerade an der 
Stelle liegt, wo die Plasmaschwingungen nach der Theorie zu erwarten sind. 


1. Einleitung 
Nach der phänomenologischen Theorie der Supraleitung Londons und von 
Laues!) steht die Eindringtiefe magnetischer Felder in Supraleiter in folgender 
Beziehung zur Elektronendichte: 


In dieser Gleichung ist n die Konzentration der freien Elektronen, m die Elek- 
tronenmasse und e die Elementarladung. Man hat also die Möglichkeit, aus der 
Elektronenkonzentration auf die Eindringtiefe magnetischer Felder zu schließen. 

Über die Eindringtiefe magnetischer Felder bei der Supraleitung liegen aus 
letzter Zeit zuverlässige Messungen vor, die nach verschiedenen Verfahren durch- 
geführt wurden. Untersuchungen von Shoenberg?) basieren auf der experimen- 
tellen Bestimmung der magnetischen Suszeptibilität von Suspensionen supra- 
leitender Metalle, die zum mindesten für Quecksilber Vertrauen verdienen. Aus 
ihnen geht zunächst hervor, daß die Eindringtiefe bei Annäherung an den Sprung- 
punkt über alle Grenzen wächst. Bei tieferen Temperaturen strebt der Wert der 
Eindringtiefe sehr schnell einem Grenzwert zu, den man als den Wert der Eindring- 
tiefe bei 0° K ansehen kann. Er ist bei Sn und Hg bei 2,3° K praktisch schon 
erreicht. Shoenberg findet Werte von ungefähr 7 - 10-* cm, während Pippard?) 


m! 


2en!n: 


1) Vgl.z.B. M.v. Laue, Theorie der Supraleitung, Berlin und Göttingen, 1947. 
2) D. Shoenberg, Nature 148, 434 (1939); Proc. Roy. Soc. London (A) 175, 49 (1940); 
M. Desirant, D. og te Nature 159, 201 (1947). 


3) A. B. Pippard, Nature 159, 434 (1947); Proc. Roy. Soc. London (A) 191, 399 
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5) M. Steenbeck, Z. Physik 76, 260 (1930). 
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aus N Ba des Skineffektes für Zinn und Quecksilber den Wert 
7,2-10-* cm herleitet, welcher demnach in guter Übereinstimmung mit dem 
von Shoenberg gefundenen steht. Setzt man nun in die Formel für die Eindring- 
tiefe plausible Werte für die Elektronendichte in Metallen ein, so erhält man 
Ergebnisse, die etwa um den Faktor 3,5 kleiner sind als die eben zitierten experi- 
mentellen Werte. 


Man könnte sich auf den Standpunkt stellen, daß bei dem heutigen Stand der 
Entwicklung der Theorie der Supraleitung eine bessere Übereinstimmung zwi- 
schen Theorie und Experiment nicht zu erwarten sei. Wir möchten jedoch ver- 
suchen zu zeigen, daß wir eine Diskrepanz genau des gleichen zahlenmäßigen Be- 
trages vorweisen können an einer Erscheinung, welche normalerweise nicht mit 
dem Phänomen der Supraleitung in Verbindung gebracht wird, nämlich der Metall- 
absorption im kurzwelligen Sichtbaren und langwelligen Teil des ultravioletten 
Spektrums. Hierbei greifen wir zurück auf Vorstellungen über das Zustande- 
kommen der Supraleitung, die wir in zwei Arbeiten *) (im Folgenden mit I und II 
bezeichnet) bereits an gleicher Stelle veröffentlicht haben, und aus denen sich 
ein einfacher Zusammenhang zwischeh der Eindringtiefe magnetischer Felder bei 
der Supraleitung und der Lage eines bestimmten Teiles der Absorption des be- 
treffenden Metalles herleiten läßt. 

Es wird allgemein die Auffassung vertreten, daß die Absorption des kompakten 
Metalles durch zwei wesensverschiedene Mechanismen bestimmt wird, die zu zwei 
spektral aneinander anschließenden, aber doch immerhin gut abgesetzten Spektral- 
gebieten Anlaß geben. Im Rahmen der Bändervorstellung wird das langwelligere 
Absorptionsgebiet gedeutet als Übergang der Elektronen im gleichen Band, wobei 
die Auswahlregeln für die Impulse durch Wechselwirkung mit den Schallquanten 
erfüllt werden. Dies ist auch das Gebiet der „klassischen“, mit der elektrischen 
Leitfähigkeit eng zusammenhängenden Absorption, deren wesentliches Kennzeichen 
darin besteht, daß sie ebenso wie der elektrische Widerstand mit abnehmender 
Temperatur verschwindet. Kurzwelliger liegt die Absorption, die dem Übergang 
eines Elektrons aus dem Leitfähigkeitsbande in ein höheres Energieband zuzu- 
schreiben ist. Diese ist grundsätzlich nicht temperaturabhängig, da sie äuch ohne 
Wechselwirkung mit Schallquanten möglich ist, Sie wird im Folgenden kurz 
Band-Band-Absorption genannt. 

Es ist zuerst wohl von Steenbeck®) darauf hingewiesen worden, daß die 
Ultraviclettabsorption der meisten Metalle etwa an der Stelle liegt, wo unter 
Zugrundelegung der Langmuirschen Formel für die Plasmaschwingungen eines 
Elektronengases die kleinste Eigenfrequenz der Leitungselektronen eines Metalles 
liegen könnte. Wir haben schon früher®) die Ansicht geäußert, daß damit ein dritter, 
von den beiden obengenannten unabhängiger Mechanismus der Absorption durch 
die Leitungselektronen eines Metalles gegeben ist. Hierbei blieb zunächst einmal 
die Frage offen, ob die experimentellen Untersuchungen Raum lassen für den 
Einfluß eines solchen dritten Mechanismus und ob überhaupt die Oszillatoren- 
stärken der Plasmaabsorption nennenswerte Beträge erreichen können. 

In sehr sorgfältigen Untersuchungen hat Wolter?) auf Grund von Messungen 
von Goos zeigen können, daß im kurzwelligen sichtbaren Gebiet eine Zunahme 


= 
4) F. Möglich, R. Rompe, Ann. Physik (6) 1, 27 (1947); 3, 322 (1948). 


6) F. Möglich, R. Rompe, l.c. 
7) H. Wolter, Z. Physik 105, 269 (1937). 
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der Absorption mit kleiner werdender Schichtdicke des Metalles auftritt. Dieser 
Absorption hat man offenbar bisher keine wesentliche Bedeutung beigemessen, 
um so mehr wohl, weil bei mittleren Schichtdicken sich auch im ultraroten Gebiet 
eine Zunahme der Absorption feststellen läßt, die allerdings bei abnehmender 
Schichtdicke wieder zu verschwinden scheint. In jüngster Zeit hat auf Anregung 
der Verfasser in diesem Zusammenhang F. Havemann die Absorption verschie- 
dener Metalle in kolloidalen Lösungen untersucht. Die Metalle zeigen überein- 
timmend eine selektive Absorption im kurzwelligen sichtbaren Gebiet, die in 
voller Übereinstimmung mit den Ergebnissen von Wolter steht. Bei Ag und Hg 
beginnt bei unfefähr 3600 Ä ein Abfall der Absorption, der bei Ag nach einem 
deutlichen Minimum bei 3500 Ä in einen Wiederanstieg des Absorptionskoeffizienten 
übergeht. Dieser Anstieg ist ziemlich steil. Im Gegensatz hierzu verhält sich 
Quecksilber anders. Der Abfall der Absorption, der etwas kurzwelliger liegt, ist 
sehr steil und führt dazu, daß unterhalb von 3000 Ä das Quecksilbersol vollständig 
durchlässig wird. Von einem Wiederanstieg der Absorption ist bis zu 2700 Ä 
nichts zu merken. 

Unserer Auffassung nach sprechen diese Befunde beinahe zwingend dafür, den 
Mechanismus der Plasmaschwingungen für das bei kleinen Schichtdicken bzw. 
Partikelgrößen auftretende Absorptionsmaximum verantwortlich zu machen. Bei 
Silber und Quecksilber liegt dieses Maximum zwischen 4800 und 3400 A) auf der 
langwelligen Seite verfälscht durch die sogenannte klassische Absorption des 
Elektronengases. Der Wiederanstieg der Absorption beim Silber im Ultraviolett 
dürfte wohl der Band-Band-Absorption zuzuschreiben sein. Jedenfalls stimmt die 
Lage dieser Absorption mit der ‘berein, die man bisher bei Ag den Band-Band- 
Übergängen zugeschrieben hat. Es hat also damit den Anschein, daß tatsächlich 
das Spektrum eines Metalles durch das Zusammenwirken dreier Mechanismen be- 
stimmt wird, von denen der -dritte, bisher nicht beachtete, offenbar die Plasma- 
schwingungen sind. wi 

Auf Grund der theoretischen Überlegungen in I und II können wir die Ab- 
sorptionsstärke dieser Plasmaschwingungen abschätzen. Es ergibt sich, daß bei 
Schichten normaler Dicke die Absorption nicht zustande kommt, weil eine Impuls- 
behinderung auftritt. Diese Impulsbehinderung wird dadurch hervorgerufen, daß 
die Plasmaschwingungen, verglichen mit den elektromagnetischen Wellen, bei 
gleicher Frequenz immer einen kleineren Impuls haben als die elektromagnetischen 
Wellen und daß daher die gleichzeitige Erfüllung eines Energie- und Impuls- 
erhaltungssatzes unmöglich ist. Es zeigt sich jedoch, daß der Impulserhaltungssatz 
mit um so geringerer Genauigkeit zu gelten hat, je kleiner die lineare Ausdelmung 
der metallischen Schichten ist. Während die Stärke der Absorption dieser Plasma- 
schwingungen nur !/,,—!/» der Band-Band-Absorption der Metalle beträgt, wenn 
die Schichtdicken 10° cm überschreiten, kann die Plasmaabsorption die Band- 
Band-Absorptionsstärke erreichen, wenn die Ausdehnung des Metalls beträchtlich 
kleiner Als 10-5 cm ist. Ein Blick auf die Abb. 1 zeigt, daß diese Verhältnisse in 
der Natur gerade realisiert sind®). 


spapnungen ein entsprechendes Verhalten. Bei Silber liegt die Resonanzspannung bei 


3,7 Volt und bei Quecksilber bei 4,86 Volt. lea - 


| 
| a 
| 
ae _ ) Bei dieser Deutung der Versuchsergebnisse erscheint uns besonders befriedigend 4 
die Tatsache, daß die Band-Band-Absorption von Silber und Quecksilber hinsichtlich a 
der spektralen Lage beträchtliche Unterschiede zeigt. Es mag darauf hingewiesen werden, aq 
daß ja die Ionisierungsspannung der Atome von~Hg und Ag sich ebenfalls beträchtlich 4 - 
unterscheiden (7,54 Volt bei Ag und 10,4 Volt bei Hg). Auch zeigen die Resonanz- a 
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Auf die Anomalien der Absorption im Bereich der sichtbaren Strahlung im 
Zusammenhang mit den Absorptionsmessungen von Goos hat iibrigens auch 
David®) hingewiesen. Er macht bereits den Versuch, die vom bisherigen Stand- 
punkte der Festkörperphysik aus nicht zu erklärende Absorption auf gewisse 
Eigenschwingungen im sichtbaren und ultraroten Gebiet zurückzuführen. Seine 
Untersuchungen sind aber nicht so überzeugend wie die Rechnungen von Wolter, 
auf die wir schon hingewiesen haben und die den Effekt der kurzwelligen sichtbaren 
Absorption ganz eindeutig unter Beweis stellen. 

Die Untersuchungen von David haben dann H. Jensen”) veranlaßt, der 
Frage nach der Herkunft der von David geforderten EigenschWmgungen nach- 
zugehen. Jensen basiert bei seinen Untersuchungen auf einer Arbeit von Bloch"), 
der wohl erstmalig darauf hingewiesen hat, daß bereits ein normales Atom. durch 
die Möglichkeit der Schwankung seiner elektrischen Dichte zu Schwingungen be- 
Dr fähigt sein muß. Bloch benutzt diesen Gedanken, um die Bremsfähigkeit schwerer 


tragung auf diese Schwingungen bei der Bremsung von Partikeln berücksichtigt). 

Jensen untersucht das dynamische Gleichgewicht eines Fermischen Atoms, 
wobei er sich allerdings im wesentlichen darauf beschränkt, die Schwankungen 
der Dichte in einem solchen Gasplasma zu berechnen, während ihm die Existenz 
der tranWersalen Schwingungen entgeht. Er wendet seine Überlegungen ausdrück- 
lich auf die Schwingungen an, die bei kleinen Metallkiigelchen zu erwarten sind, 
wobei er sich auf den Fall konstanter räumlicher Elektronendichte beschränken 
muß. Er findet im wesentlichen eine Bande von für die Absorption in Frage 
kommenden Eigenschwingungen im Grenzgebiet zwischen dem Ultravioletten und 
Sichtbaren. 

Wie wir im Laufe der Arbeit zeigen wollen, wird gerade die langwellige Kante 
der Plasmaabsorption in Beziehung gesetzt werden müssen zur Eindringtiefe 
magnetischer Felder im supraleitenden Zustand@. Es würde also alles darauf 
ankommen, eine möglichst genaue Angabe über die Lage dieser langwelligen 
Absorptionskonstante auf Grund des experimentellen Materials zu machen. Das 
ist aber, wie ein Blick auf die Abb. 1 zeigt, heute noch nicht möglich, weil das 
Material noch zu spärlich ist. Auch macht die Überdeckung durch die Vorgänge 
im Ultraroten, die, wie es den Anschein hat, nur zum Teil auf der klassischen 
Absorption des Elektronengases beruhen, eine genauere Abschätzung der Plasma- 
absorptionskante sehr schwer. Dabei ist überhaupt zu beachten, daß Rückschlüsse 
aus Messungen des Produkts n kv!?) auf Absorptionsstellen sehr schwer sind, 
worauf bei der Besprechung der Band-Band-Absorption bereits von Fröhlich, 


10) H, Jensen, 2. Physik 106, 620 (1937). 
a 4) F. Bloch, Z. Physik 81, 363 (1933). 

22) In diesem Zusammenhang mag es interessant sein, darauf hinzuweisen, daß durch 
i den gleichen Gedankengang Langmuir seinerzeit auf die Entdeckungen der Plasma- 
schwingungen bei Gasentladungen gefiihrt wurde. Die normalen Bremswege von in eine 
Gasentladung geschossenen Elektronen waren stets beträchtlich kleiner, als auf Grund der 
damaligen theoretischen Ergebnisse zu erwarten war. Daraus schloß Langmuir, daß 
ein bisher noch nicht bekannter Prozeß der Energieübernahme vorhanden sein muß, als 
den er die Schwingungen der Plasmamasse des Gases erkannte. Langmuir konnte 
dann sehr bald experimentell nachweisen, daß diese Schwingungen auch Strahlung aus- 
zusenden imstande sind, wobei es ihm gelang, die Richtigkeit der nach ihm benannten 

Formel für die Frequenz dieser Schwingungen zu verifizieren. , 


00.8) n und & sind Real- und Imaginärteil des komplexen Brechungsexponenten. 


x 9) E. David, Z. Physik 114, 389 (1939). 


Atome gegenüber Korpuskularstrahlung zu berechnen, indem er die Energieüber- , 
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Sommerfeld und Bethe hingewiesen wurde. Das dort Gesagte gilt genau ent- 
sprechend für die Beurteilung der Lage der Plasmaabsorption. Das vorliegende 
experimentelle Material könnte etwa ausgedeutet werden durch die Annahme, daß 
die langwellige Absorptionskante der Plasmaschwingungen je nach Art des Metalles 
zwischen 4000 Ä und 5000 Ä liegt. 


Wenn man nun diesen Wert der Absorptionskante mit denjenigen Werten 
vergleicht, die sich aus der Langmuirschen Formel (I, 5) ergeben, so findet man 
für alle Metalle deutlich zu kleine Werte, für Ag und Au, die beide fast die gleiche 
Elektronendichte haben, Werte zwischen 1300 und 1400Ä. Es liegen danach 
nach unserer Ansicht die wirklich auftretenden Plasmaabsorptionskanten um den 
Faktor 3,5 nach längeren Wellen verschoben, um den gleichen Faktor also, um 
den sich die Eindringtiefen magnetischer Felder nach den neuesten Messungen 
von dem von London berechneten Wert unterscheiden. 


Es hat damit also den Anschein, als ob die von uns aus theoretischen Gründen 
postulierten Terme der Plasmaschwingungen, die zwangsläufig zu einer Theorie 
der Supraleitung führten, nicht nur vorhanden sind, sondern sich auch optisch 
im sichtbaren Gebiet der Metallabsorption bemerkbar machen. Damit würde also 
das Spektrum der Metalle um einen charakteristischen Anteil vermehrt werden, 
der bisher nicht beachtet worden ist. Dieses Ergebnis erinnert an einen Deutungs- 
versuch der Supraleitung durch -F. London, der darauf aufmerksam gemacht hat, 
daß die Supraleitung verstanden werden kann, wenn man die Existenz’ gewisser 
neuer Elektronenterme annimmt, über deren Eigenschaften er genaue Angaben 
machen konnte. Es ist nicht von der Hand zu weisen, daß die von uns gefundenen 
neuen Terme die von London postulierten Eigenschaften haben. Wenn das der 
Fall wäre, so würde ein tieferer Zusammenhang zwischen dem Londonschen 
und unserem Deutungsversuch zur Supraleitung vorhanden sein. 


2. Eindringtiefe magnetischer Felder und Absorptionskante 


Die Existenz der Plasmawechselwirkung hat, wie wir in (II) zeigen konnten, 
eine Relation zwischen Strom und Vektorpotential zur Folge, die der bekannten 


wir ähnlich ist: 
(IL, 5) 


2 


Der Unterschied gegenüber der Londonschen Relation besteht nur darin, daß 
die Elektronendichte n, nicht mehr räumlich konstant ist, sondern der Gleichung 


(II, 7) 
No 

genügt. Hier ist @ das elektrostatische Potential des Kristalls und Z, die Null- 

punktsenergie als Funktion der Dichte n,, wobei man im einfachsten Falle für Z, 

den Sommerfeldschen Wert für die Nullpunktsenergie eines Elektronengases 

wählen kann. Die Gleichung (II, 5) unterscheidet sich von (I, 3) nur dadurch, 

daß in (1,3) n eine räumliche Konstante ist. 


Ebenso, wie man in (1) mit Hilfe von (I, 3) zur Gleichung (I, 4) übergeht, die 
dort gleichzeitig nicht nur für den Strom j, sondern auch das Vektorpotential U 
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gilt, kann man mit Hilfe von (II, 5) zu einer Differentialgleichung übergehen, die 


jetzt wegen der räumlichen Inkonstanz von n, nur noch für X gilt: 


Fragt man nach zeitlich periodischen Lösungen, so gelangt man zu der Eigen- 
wertgleichung 
v2 Ng (x yz)e* 

Af + (3 0. (2) 
Der kleinste dieser Eigenwerte ist dann, abgesehen vom Faktor 42, die Reziproke 
des Wellenlängenquadrates der gerade noch von Kristall absorbierten Strahlung. 
Ist n, konstant und der mittleren räumlichen Dichte der Elektronen im Mittel 
gleich, so erhalten wir für die Grenzfrequenz genau die Langmuirsche Formel, 
die also für Ag einen Wert von 1350 Ä für die Grenzwellenlänge ergibt. 

Experimentell erhält man aber einen um etwa 3,5mal größeren Wert für die 
Grenzwellenlange. Daß man mit der Langmuirschen Formel den kleinsten 
Eigenwert von Gleichung (3) nur größenordnungsmäßig wiedergeben kann, darf 
uns aber nicht verwundern. Wir können umgekehrt, solange wir über den Verlauf 
von n, gar nichts wissen, die experimentellen Ergebnisse als eine Vervollständigung 
der Gleichung (2) ansehen. ‚Die Langmuirsche Formel stimmt deswegen mit 
der Erfahrung nur ungenau überein, weil die Voraussetzung, die zu ihrer Ableitung 
gemacht wurde, nicht mehr zutrifft: die Elektronen im Metall sind nicht mehr 
als frei zu betrachten. Die Abweichung zwischen der Langmuirschen Formel 
und der wirklich vorliegenden @renzfrequenz kann also ein Maß für die Gebunden- 
heit der Elektronen im Metall angesehen werden. 

Übrigens haben wir bisher stillschweigend angenommen, daß die Dielektri- 
zitätskonstante ¢ infolge des Vorhandenseins der stärker gebundenen Elektronen, 
die nicht zum Plasma gehören, gleich 1 gesetzt werden darf. Das ist, wie bereits 
Jensen!) bemerkt hat, nicht unbedingt richtig. Man würde wohl eine Dielektri- 
zitätskonstante zu erwarten haben, die größer als 1 ist. Ihr Vorhandensein würde 


eine Verkleinerung der kleinsten Eigenfrequenz », um den Faktor .. hervorrufen, 


Ve 
wie man sich unschwer iiberlegen kann. 

“ Wenn wir nun die Eindringtiefe magnetischer Felder in einen Supraleiter be- 
rechnen wollen, so haben wir Gleichung (1) mit zeitlich konstantem A zugrunde 
zu legen. Wir können also die Eindringtiefe aus folgender Gleichung berechnen: 

2 
AM — _ 9, 
C 
Nehmen wir an, daß die Oberfläche des Supraleiters eine Ebene ist und daß die 


x-Achse senkrecht zu dieser Ebene steht (x = 0 sci diese Ebene, x > 0 das Gebiet 
des Supraleiters), so können wir diese Gleichung mit dem "Ansatz 


wo ß-! also die Eindringtiefe bedeutet, in die folgende iiberfiihren: = 


Ab — 2p = +(B — 0, 
4) H. Jensen, l. e. N 
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Wenn wir nun noch annehmen, daß ® eine schwach variierende Funktion ist, so 
daß das Glied B= = klein ist gegen 6? ®, was wegen der Kleinheit der Eindringtiefe 
p-* fast immer erfüllt sein dürfte, so geht Gleichung (3) offenbar in Gleichung (2) 
über und wir finden für die Eindringtiefe 


Durch diese Gleichung wird ein ganz neuer Zusammenhang festgestellt zwischen 
zwei Größen, die bisher in keiner Beziehung zueinander standen: 


Eindringtiefe magnetischer Felder bei Supraleitung und Grenzwellen- 
länge der kurzwelligen Absorption im Sichtbaren des betreffenden 
Metalles sind zueinander proportional. Der Proportionalitätsfaktor 
ist —— 

Wenn man nun versuchen will, aus der Gleichung (II, 5) eine „Londonsche 
Relation‘ zu machen, so kann dies nur durch Mittelwertbildung geschehen. Man 
würde dann eine Gleichung zwischen dem makroskopisch meßbaren Suprastrom % 
und dem Vektorpotential X erwarten von der Form 

wo A die Londonsche Konstante, also eine räumlich konstante Größe, ist. Sie 
steht nach der phénomenologischen Theorie von London und v. Laue mit der 
reziproken Eindringtiefe 8 in folgender Beziehung: 


1 

‘In unserer Theorie kommt der Londonschen Konstanten keine unmittelbare Be- 

deutung zu. Auf dem Umweg über die Eindringtiefe aber können wir doch die 

Londonsche Konstante mit der Grenzfrequenz der Ultraviolettabsorption eines 
Metalles in Verbindung bringen: 4-1 = 2%. 

Die Experimente von Shoenberg und Pippard schließen eindeutig das 


‘ Zutreffen einer Eindringtiefe der Größe - wie sie sich aus dem ursprüng- 


cm 

at nt’ 
lichen Londonschen Wert m/ne? fiir A pate? al wiirde, als zu klein aus. Die ge- 
messenen Werte sind um den Faktor 3 bis 4 gréBer als die Werte der Londonschen 
Theorie. Die gemessenen Werte sind um genau denselben Faktor größer wie die 
wirklichen Grenzwellenlängen der Plasmaabsorption der Metalle, verglichen mit 
den Werten, die man erhält, wenn man die Langmuirsche Formel für ,,freie“ 
Elektronen zugrunde legt. Die Ergebnisse stimmen gut mit unserer Behauptung 


von dem Zusammenhang von Ultraviolettabsorption und Supraleitung überein. 


Im Hinblick auf dieses Ergebnis scheint es doch schon angebracht zu sein, 
einmal auszurechnen, wie groß die Wahrscheinlichkeit der Absorption von Licht- 
quanten durch ein Elektronenplasma ist. Um dies zu tun, brauchen wir einen 


\ 


3. Die Ubergangswahrscheinlichkeiten der Plasmaabsorption Gr 
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pre ne für die Energie der Wechselwirkung zwischen Plasma und elektro- 


magnetischem Feld. Stellen wir das elektromagnetische Feld durch sein Vektor- 
potential Yj, dar, das so gewählt werden soll, daß seine Divergenz verschwindet, 
so können wir also ansetzen, wenn wir die x-Achse in die Richtung der Fortpflan- 
zung der ebenen Welle, die y-Achse bzw. z-Achse je nach der Polarisation mit der 
Richtung von zusammenfallen lassen, 


d,. € + dye (6) 
Als Ene rgiedie hte machen wir mit Hilfe des durch die Gleichung (II, 5) definierten 
Stromes j den Ansatz: BE 
: e ech 


Indem wir weiter beachten, daß die Eigenschwingungen des ‚inneren‘ Vektor- 
potentiales A der Gleichung (2) genügen, können wir für j folgenden Fourier- 
Ansatz machen: 


Hier sind die p, komplexe Impulsvektoren harmonischer Oszillatoren, die zwar 


nicht mehr Funktionen der Raumkoordinaten 2, y, z, wohl aber der Zeit t sind; 
und zwar ist ihre zeitliche Abhängigkeit durch e*+?”irt gegeben. Nunmehr er- 
halten wir für V,,: 


vz vz 


Mit diesem ergibt sich der Ausdruck fiir die Wechselwirkungsenergie, indem wir 
über den ganzen Raum integrieren 


mc 

(8) 
Bevor wir die Integrale, die hier auftreten, vereinfachen kénnen, miissen wir 
noch über die Funktionen f, genauere Kenntnisse haben. Dazu ist es nur notwen- 
dig, darauf hinzuweisen, daß die Differentialgleichung (2) in der Funktion n, (2 yz) 
eine im Gitter des Metalles räumlich periodische Funktion enthält. Formal stimmt 
daher die Differentialgleichung (2) mit der Schrödingerschen Differential- 
gleichung der Eigenfunktionen eines Kristalles mit periodischem Potential überein, 
nur mit dem Unterschied, das der Eigenwert hier nicht mehr die Frequenz der 
Welle, sondern das Quadrat der Frequenz der Welle ergibt. Insbesondere dürfen 

wir schließen, daß die Funktionen f, die folgende Form haben: 


— F, (x y 2), (9) 


wo F, eine im Gitter des Kristalles periodische Funktion bedeutet. Unter §t ver- 
stehen wir den ,,reduzierten“ Ausbreitungsvektor, der immer so gewählt werden 
kann, daß seine Projektionen auf die Kristallachsen a, a, a, zwischen —} und +} 
liegen. 
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tae Uber die Verteilung der möglichen Eigenwerte der Gleichung (2) darf man 
Dur die gleichen Schlüsse ziehen, wie im Falle der Eigenwerte und Eigenfunktionen 
‚det in einem periodischen Potential. Die Eigenwerte liegen also auch hier kontinuier- 


am: lich und es wechseln erlaubte und verbotene Gebiete ab. Also auch hier zeigt 
ee sich eine Bandstruktur fiir die möglichen Werte von » der Gleichung (2). 
Man kann jedoch leicht überlegen, daß im Gegensatz zu den Ergebnissen für ? 
die Eigenwerte und Eigenfunktionen im periodischen Potentialfeld hier bei unseren 
(6) Plasmafunktionen nur ein einziges, und zwar das unterste Band!®) eine physikalisch * 
Realität haben kann. Dies sieht man so ein: Die eine Kante des untersten Bandes 
wird ein Eigenwert derjenigen Eigenfunktion sein, für die in Gleichung (9) & 
gleich Null gesetzt werden muß. Diese Funktion ist also im ganzen Gitter peri- 
odisch. Wir können sie dadurch charakterisieren, daß wir nach einer Lösung der 
7 Gleichung (2) fragen, die so beschaffen ist, daß sie innerhalb einer geeignet ge- 
(4) wählten Elementarzelle niemals verschwindet und außerdem auf der Oberfläche 
-tor- der Elementarzelle der Randbedingung = genügt. Wenn wir uns diese Funktion 
wil? in jeder Elementarzelle reproduziert dutta, erhalten wir im ganzen Gitter eine 
periodische Funktion, die dank der Randbedingung sich stetig differenzierbar 
von einer Elementarzelle zur anderen fortsetzt. x 
Die andere Kante des Bandes gehört im allgemeinen zu einer Funktion, bei 
der die in (9) auftretende e-Funktion in zwei unmittelbar benachbarten Elementar- 
War zellen sich nur noch durch das Vorzeichen unterscheidet. Eine solche Funktion 
ind; erhält man offenbar, indem man eine Lösung von (2) verlangt, die innerhalb einer 


Elementarzelle nicht verschwindet, aber auf der Oberfläche der Elementar- 
zelle der Randbedingung f = 0 genügt. Sorgen wir dann dafür, daß in zwei 
aneinander anstoßenden Elementarzellen diese Funktionen mit verschiedenen 
Vorzeichen aneinander gesetzt werden, so haben wir damit eine Funktion im 
gesamten Kristall, die die Form (9) besitzt. Offenbar ist dies eine Funktion, 
wir deren Wellenlänge von der Größe des Abstandes zweier benachbarter Elektronen 
wird. Also ist dies die kleinste Wellenlänge der Plasmaschwingungen, bei der unsere 
kontinuierliche Beschreibung des Plasmas trotz der atomistischen Struktur der 


it. Elektronen allenfalls noch einen Sinn haben kénnte. Die Eigenfunktionen, die 
zu höheren Bändern gehören, haben Wellenlängen, die kleiner sind als dem mitt- 
(8) leren Abstand zweier Elektronen entspricht, und sind daher nicht mehr zu berück- 
wir sichtigen. Daher ist nur das unterste Band der Plasmaeigenfrequenzen in der * 
ven- Natur verwirklicht. 

YZ) Wenn wir (9) in die Integrale von Gl. (8) einsetzen, so erhalten wir, wenn wir 
nmt bedenken, daß die Werte von r an zwei entsprechenden Stellen verschiedener ° , 
tial- Kristallzellen sich um den Vektor (rect 

ein, 

dee Ny Ay + Nz Az + Nz Ag 

rfen mit ganzzahligen n,, unterscheiden: - 


f no fy dz — [m F, Se , 


(9) 
ver- wobei die Integration nunmehr nur noch iiber eine einzige Elementarzelle zu er- 
den strecken ist. Die & hingegen ist über alle Elementarzellen auszudehnen. 
n 
5) Dieses Band kann allerdings u. U. in aaenean voneinander getrennte Teil- 
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Nun ist aber die Summe iiber die n, eine zwar endliche Summe, die über sämt- 
liche Elementarzellen des Kristallgitters zu erstrecken ist, bei der die Zahl der 
Summanden jedoch sehr groß ist. Unter ihnen werden im allgemeinen immer je 
zwei Glieder sein, die praktisch den gleichen Wert, aber das entgegengesetzte Vor- 
zeichen haben. Die Summe wird daher bei sehr großer Summenzahl praktisch 
gleich Null zu setzen sein, mit alleiniger Ausnahme des Falles, daß alle Glieder 
gleich +1 sind. Dies tritt dann und nur dann ein, wenn 


(10) 
wo g,, ein „ganzzahliger‘‘ Vektor des reziproken Gitters bedeutet: a8, 


Gm = mıb, + m,b, + m, 
mit 
und so fort, wobei wir für b, und b, entsprechende Gleichungen durch zyklische 
Vertauschungen der Indizes 1, 2, 3 erhalten. Bezeichnen wir mit Y„ das über 
eine Kristallelementarzelle erstreckte Integral ug: 
Yu = f (11) 
(Zelle) 


so geht die Wechselwirkungsenergie E, aus Gleichung (8) über in 


a N: Vn fir +f= Im 
(12) 
fir — f= g,, 


m 


[Pv U. N- y 
\P, 


vz 


Hier bedeutet N die Anzahl der im Kristall vorhandenen Zellen. 

Wir wollen nun diskutieren, inwieweit die Gleichung (10) eine einschränkende 
Bedingung für den Prozeß der Absorption elektromagnetischer Strahlung durch 
ein Elektronengas ist. Hierzu müssen wir die Größenordnung von §t, dem Aus- 
breitungsvektor der Plasmaschwingungen, kennen. § kann nach der Definition 
als „‚reduzierter‘‘ Ausbreitungsvektor zwar nie größer als etwa 10® werden, jedoch 
könnte er natürlich weit unter dieser Schranke bleiben. Der Ausbreitungsvektor 
¥ als Funktion der Frequenz » läßt sich allerdings nur für den Fall eines „freien“ 
Plasmas berechnen, d. h. eines solchen, das wir in Arbeit I zugrunde gelegt haben. 
Obgleich wir der Ansicht sind, daß die dort gemachten Voraussetzungen zu grob 
sind, wollen wir doch annehmen, daß die Dispersion der Plasmawellen angenähert 
wenigstens durch die Dispersion eines freien „Plasmas“ wie in I beschrieben 
werden kann. a 

Aus der Gl. (I, 4) entnimmt man, wenn man für j den Ansatz 


macht, das Dispersionsgesetz der Plasmawellen = 
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so daß mit dieser Gleichung (10) übergeht in 
VR — 


c 


Es ist also stets In] < |tl bezogen auf gleiche Frequenz. 
Der Vektor g,, wird von drei Vektoren b, aufgespannt, von denen jeder einen 
absoluten Betrag der Größe 10® hat. Ist also in g,, auch nur eine der drei Zahlen 
des Tripels m,, m,, m; von Null verschieden, so ist jedenfalls |g,,] > |f| und damit 
auch gegen |§t|, da |f| bei Strahlung im Sichtbaren und Ultravioletten höchstens 
von der Größenordnung 10° ist. Daher kann Gl. (10) bzw. Gl. (14) nur für g,, = 0 
erfüllt werden. 
Wir werden ferner zu n bestichsichtigen haben, daß neben (14) auch noch de 
Energiesatz 


hy=E, E,, : 


zu gelten hat, wo die E, Energieniveaus harmonischer Oszillatoren sind. Da die 
Differenz der Energieniveaus E,,—£,, beim harmonischen Oszillator, zwischen 
denen Übergänge méglich sind, gerade gleich hy ist, so ergibt der Energiesatz 
»—=». Wenn man dies in Gl. (14) einsetzt, so erkennt man, daß sie jedenfalls 
für Strahlung im sichtbaren und ultravioletten Gebiet nicht zu erfüllen ist. 
Andererseits dürfen wir nicht vergessen, daß Gl. (10) nur annäherungs- 
weise erfüllt zu sein braucht. Es liegt dis daran, daß die endliche Summe 


S— Inga, - (15) 
- n e 
= 
die wir in GI. (12) mit dem Buchstaben N zusammengefaßt haben, auch dann 
noch fast den Wert N ergibt, wenn $ + f nicht gleich Null oder g,,, sondern gleich 
QJ, + 6 ist, wo 6 nur hinreichend klein zu sein braucht. Es ist nur nötig, dafür zu 
sorgen, daß das Produkt aus 6 mit dem Vektor = N, A, mit dem größtmöglichen 


Wert von n, immer noch klein gegen 1 ist. 


Die Größenordnung der Länge der Vektoren a, ist 10. Der Wert von if 
ist für eine Strahlung der Wellenlänge 5000 A gerade gleich 2 - 10°. Nehmen wir 
an, daß wir uns in der Nähe der Absorptionskante der Plasmaschwingungen 
befinden, so ist |] noch sehr klein und |§t —£| von der Größenordnung von ||, 
d.h. gleich 2 - 10*. Nehmen wir weiter an, daß das Metall aus einer Schicht von 
T Atomlagen bestehe, und daß eine Strahlung senkrecht zu dieser Schicht, also 
in Richtung der X-Achse, diese Schicht durchdringe. Der Einfachheit halber 
wollen wir dann ferner annehmen, daß der kristalline Aufbau der metallischen 
Schicht kubisch ist und daß einer der Vektoren dp, z. B. a,, die untereinander vom 
gleichen Betrage sind und aufeinander senkrecht stehen, in die Richtung der 
x-Achse falle. Unter den gemachten Voraussetzungen ist also nur (® — f) a, von 
Null verschieden, während beide Produkte (® — £) a, und ( — f£) a, verschwinden. 
Daher geht die Summe Gleichung 16 in ein Produkt von 3 Summen über, von 
denen zwei Summen nur aus Summanden der Größe ‚1 bestehen. Das Produkt 
dieser beiden Summen ergibt also offenbar gerade die Zahl der Atome, die sich 
in einer monoatomaren Schicht des durchstrahlten Metalles befinden. Nennen wir 


die Zahl der Atome längs der y-Achse U, längs der z-Achse V, so ist das Produkt 
der genannten Summen gerade gleich U-V. 


23* 


4 
der 
‘or- (14) 
sch { 
der 
(10) 
__| 
| 
4 
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vee Die dritte Summe nimmt die Form 


N, 


an. Wenn wir die Zahl der atomaren Schichten in der Richtung der 2-Achse mit 7 
bezeichnen, so ist der Wert der Summe, abgesehen von einem komplexen Faktor 
„sin — Ha] 
vom Betrage 1, gegeben durch sin Ha] 


und das Produkt der drei Summen 
geht über in 2 > | 
, sin [Ta(8 — ha] \ 


sin — Ha] ’ 


(16) 
a 
wobei wir 7'-U-V, die Gesamtzahl der Teilchen der Schicht, wieder mit N 
bezeichnet haben. 
Die in der Gleichung (16) auftretende Funktion ist in der Optik wohl bekannt. 
Sie hat bei hinreichend großem 7' im allgemeinen sehr kleiné Werte und nimmt 
nur den Wert 1 an, wenn die Gleichung 


(@ —t)a, = +2,43... 


erfüllt ist. In dieser Größe S äußert sich also der Einfluß der Schichtdicke. Die 
Größe S ist für die hier interessierenden Größenordnungen von $ und £ gleich 
N, solange die Schichtdicke etwa 20 Atomlagen nicht übersteigt. Bei etwa 
1000 Atomlagen ist S nur noch ein knappes Zehntel von N. 

Ist also die Schichtdicke größer als 1000 Atomlagen, so besteht wenig Aus- 
sicht, die Existenz der Plasmaschwingungen durch optische Prozesse nachzu- 
weisen. Anders ist es allerdings, wenn wir metallische Schichten kleiner Schicht- 
dicke oder Metallkolloide nicht zu großer Teilchen benutzen. Dann kann man 
erwarten, Absorptionsprozesse nachzuweisen, die intensitätsmäßig der Ultra- 
violettabsorption der Metalle in nichts nachstehen. 

Hier muß noch nachgetragen werden, inwieweit der Ansatz (9) für die Eigen- 
schwingungen des Plasmas in Einklang steht mit der Annahme endlicher Kristalle, 
über die die Summe S in (15) ausgeführt werden soll. Die in Gleichung (9) einge- 
führten Funktionen entsprechen offenbar quasiebenen durch den Kristall durch- 
laufenden Wellen, die in einem endlichen Kristall gar nicht existieren können. 
Um diesen Widerspruch zu überwinden, denken wir uns eine unendliche Zahl 
endlicher Kristalle der in Wirklichkeit vorliegenden Form zu einem unendlichen 
Kristall vereinigt. Dann fordern wir, daß die in (9) eingeführten Funktionen sich 
in allen diesen Kristallen periodisch immer wieder reproduzieren. Dies bedeutet 
also, daß wir nicht nur für die Funktionen F in (9) Periodizität für die einzelnen 
Elementarzellen zu verlangen haben, sondern daß dazu noch eine Forderung der 
Periodizität der Funktion e?7‘(*® in den einzelnen Makrokristallen gestellt wird. 
Von diesen Makrokristallen werden wir also die Annahme machen, daß die Zahl 
der Elementarzellen in der Richtung der a,-Achse gleich 7’, in Richtung der a;- 
Achse gleich U und in Richtung der a;-Achse gleich V sein soll. Durch eine solche 
Periodizitätsforderung schränken wir offenbar die Wahl des Ausbreitungsvek- 
tors ® noch weiter ein, indem wir für & folgende Bedingung finden: 

8-7 (ky, ha, = 0,41, +2...) 
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und die Beschränkung, die wir für den reduzierten Ausbreitungsvektor eingeführt 
haben, hat, wie man sich leicht überzeugt, zur Folge, daß die Zahl der nunmehr 
diskret liegenden Energieniveaus mit der Gesamtzahl der Elektronen N = 7'-U -V 
übereinstimmt. Da zwei linear unabhängige transversale Plasmaschwingungen 

tT immer zu je einer Frequenz existieren, haben wir also 2N transversale Plasma- 
tor schwingungen, die zu dem einen existierenden Frequenzband gehören. Ent- 
nats sprechend werden wir von den longitudinalen Plasmaschwingungen, die wir 
bislang nicht behandelten, noch einmal N Eigenschwingungen erhalten, so daB die 
Gesamtzahl der Eigenschwingungen, wie in (I) gefordert, wirklich gleich 3N, 
16) d.h. gleich der Zahl der Freiheitsgrade von N Elektronen ist. 
| Die Eigenenergie des Plasmas wird durch eine Summe von Energien harmo- 
y -nischer Oszillatoren mit Frequenzen vy >», dargestellt: 
| i 2 2 2 2 
E = (Pox + Poy + Prz) + 2% + Yr 
ımt 
Zur Schrödinger-Gleichung für das durch die einfallende elektromagnetische 
Welle Gleichung (6) gestörte Plasmasystem gelangen wir, indem wir zu E die 
Wechselwirkungsenergie E,, Gl. (9) hinzu addieren und dann p,,, py, Py: durch 
Die bzw. ersetzen. Man sieht sofort, daß das Matrixelement 
sich fiir die Ubergangswahrscheinlichkeit des Uberganges Plasma ohne Quant zu Plasma 
wa mit Quant der Frequenz » gegeben ist durch 
zu- Wer = J Ze (17) 
- <vjol 
nan 
tra- Hier sind die [&,]oı [%»Joı Matrixelemente der Geschwindigkeitskomponenten 
der Indizes 0, 1 harmonischer Oszillatoren von der Frequenz v. Ihre absoluten 
jen- Beträge sind also untereinander gleich. Es gilt ae 
lle, » 
len 
ahl Das als Faktor auftretende, über eine einzige Elementarzelle zu erstreckende 
hen Integral hat eine einfache Bedeutung. Da e-n,(x yz) die elektrische Dichte ist, 
ich kann das Integral zwanglos gedeutet werden als derjenige Teil der elektrischen , 
ıtet Ladung einer Elementarzelle, der an der Absorption eines Quantes der Frequenz v 
nen beteiligt ist. 
der An dieser Stelle muß noch etwas über die Normierung von F, gesagt werden. 
ird. Da die durch Gl. (9) mit F, zusammenhängenden f, zu einer Fourierentwicklung 
ahl von j benutzt werden, müssen die f, normiert sein; 
Pi 
Schreiben wir also an Stelle von Gl. (9) die folgende: Fa 4 
1 ezrian.F,, 


PL] 
A 
x 
- 
« 
+ 
a 
= 
> 


€ 
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wo wiederum N die Zahl der Elementarzellen im Kristall ist, so haben wir erreicht, 
daß F, jetzt in der Elementarzelle normiert ist 


(Zelle) 
Daraus sehen wir, daß nicht e f n) F, dr gleich der zur Frequenzy gehörenden 
(Zelle) ; 
Ladung gesetzt werden darf, sondern daß erst - . 4 
ey = n, dt =e VN f dt 


(Zelle) (Zelle) 


diese Bedeutung haben kann. Beachten wir noch, daß die in Gl. (17) auftretende 


Summe S (Gleichung 15) den Wert 7:N hat, wn= = eine reine Zahl kleiner 


N 
als 1 ist, so können wir W,, aus Gleichung (17) schreiben re 
| 
Wa=n-N-e- Im = (18) 


Das Quadrat dieser Größe gibt die Absorption an. Diese ist also proportional 
dem Quadrat der zeitlichen Änderung eines el. Dipolmomentes, gebildet mit einer 
„reduzierten‘‘ Ladung e, und proportional der Zahl der Elektronen des Kristalles. 
Die mangelnde Impulsübereinstimmung zwischen Lichtquant und Plasmaquant 
hat eine scheinbare Reduktion der Zahl der bei der Absorption beteiligten Elek- 
tronen zur Folge, die durch 7? gegeben Ast. 


| 
4. Vergleich mit dem Experiment 


In der Formel (18) ist die einzige Größe, die nicht genau bekannt ist, die Größe ey. 
Sie ist sicher kleiner als die elektrische Ladung e eines Elektrons, im allgemeinen 
aber wohl noch von der gleichen Größenordnung. Von der Größe 7? wissen wir, 
daß sie unter Voraussetzung der Gültigkeit des für die Plasmaschwingungen ange- 
nommenen Dispersionsgesetzes (13) etwa 10”? für kompaktes Material beträgt 
und auf den Wert 1 heraufgeht, wenn die Dicke der durchstrahlten Schicht kleiner 
als etwa 20 Atomlagen ist. Dann könnte man erwarten, daß die Absorption der 
Plasmaschwingungen auch von der gleichen Größenordnung wird wie die Band- 
Band-Absorption. 


In der Abb. 1 haben wir eine Abbildung aus einer Arbeit von Wolter umge- 
zeichnet. Sie stellt im wesentlichen die Größe n ky für Silber als Funktion der 
Wellenlänge dar. Wolter selbst hatte nicht nk», sondern nkvd dargestellt, 
wo d die Dicke der Schicht bedeutet. Da er aber die Dicke der Schicht auf 
Grund einer Schätzung angeben konnte, war es möglich, in seinen Ergebnissen 
die Schichtdicke noch zu eliminieren. Dabei fiel zunächst auf, daß die Band- 
Band-Absorption, die etwa bei 3300 Ä beginnt, von der Schichtdicke völlig unab- 
hängig ist. Kleine Diskrepanzen, die sich bei der Ultraviolettabsorption noch 
ergaben, haben wir benutzt, um die Angaben über die Schichtdicken ein wenig 
zu korrigieren. So glaubten wir, die Schichtdicke der dünnsten Schicht von 17 A, 
wie Wolter schätzt, auf 14 Ä reduzieren zu können, um volle Übereinstimmung 
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im Ultraviolett zu erhalten. Aus dem gleichen Grunde würden wir die dickste 
Schicht nicht, wie bei Wolter, mit 180 A, sondern mit 155A angeben. 


Man sieht zunächst, daß das Maximum im kurzwelligen Sichtbaren, das nach 
unserer Auffassung auf die Plasmaschwingungen zurückgeführt werden muß, 
ungefähr bei 4500 A liegt und daß die Absorption dort von der gleichen Stärke 
ist, wie die Absorption im Ultravioletten. Sie ist etwa 20mal so stark, wie bei 
der Absorptipn des kompakten Materials. Vermutlich wird bei der Absorption 
des kompakten Materials noch ein beträchtlicher Anteil der langwelligen 
„klassischen“ Absorption zuzuschreiben sein, so daß ein Verhältnis der Absorp- 
tionsstärken wie 1:100 ent- : 


sprechend unserer Ab- 

schätzung keineswegs zuden 6 

Tatsachen im Widerspruch 5 | 
steht. Man sieht 4 > 
keine Abnahme der Absorp- 3 va 7 / N oe 
zu Schichten von 50 oder 2 \ r 7 7 x 

60 was ebenfalls in ! I 
Übereinstimmung mit un- 

seren Abschätzungen liegt. do 

Erst bei 155 Ä dicken 7 | 
Schichten macht sich eine 
beträchtliche Abnahme der 7, N 

bar, was wiederum mit 03 
unserer Theorie in Überein- 02 L N 
stimmung ist. Außerdem "4 N 
macht sich aber hier noch 

eine Verlagerung des Maxi an 0300 500 300 500 700 5300 300 7000 


mums bemerkbar, die deut- ae: 

ich davon herrihrt, daß im der von 
roten bzw. ultraroten Gebiet Wolter. Auf der Abszisse ist die Wellenlänge in 
ein weiteres Absorptions- my, auf der Ordinate die für die Absorption charak- 
maximum als Folge der ver- _teristische Größe nkv mit dem Faktor 10”? multi- 
größerten Schichtdicke auf- pliziert aufgetragen 

getreten ist. Diese konnten 

ebenfalls von Wolter?) am Ag, Au, Pt, K und Na auf Grund des vorliegenden ex- 
perimentellen. Materials. nachgewiesen werden. Diese Maxima sind im allgemeinen 
flacher als die Maxima im kurzwelligen Sichtbaren, die nach unserer Auffassung den 
Plasmaschwingungen zugeschrieben werden müssen. Auch unterscheiden sie sich 
von diesen dadurch, daß sie sich weder bei kompaktem Material, noch bei dünnen 
Schichten bemerkbar machen, sondern offenbar an die Existenz mittlere® Schichten 
gebunden sind. Besonders bei Pt konnte Wolter diese Erscheinung sehr deutlich 
nachweisen. Wolter führt die Existenz dieser Absorptionsstellen auf Verringe- 


16) Die geringfügige Vergrößerung der Absorption bei einer Schicht von 50 Ä 
Dicke ist vermutlich auf die Erhöhung des ultraroten Anteiles der Absorption zurück- 
zuführen. 

7) H. Wolter, Z. Physik 113, 547 (1939); 115, 696 (1940). 
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rung der Dichte im Vergleich zum kompakten Material verbunden mit einer ört- 
lichen Dichteschwankung zurück. 


Besonders schön sieht man die von uns gemeint Absorptionsmaxima an den 
in der Einleitung schon genannten Versuchsergebnissen von R. Havemann, die 
wirin Abb. 2 wiedergeben. Man erkennt hier deutlich bei Silber und Quecksilber das 
Maximum im kurzwelligen sichtbaren bzw. langwelligen ultravioletten Gebiet. 
Durch Intensitätsvergleichung zwischen dem ultravioletten Maximum pei Silber, das 
von der Band-Band-Absorption herrührt und auch genau an der zu erwartenden 
Stelle liegt, und den Ergebnissen von Wolter an dünnen Schichten würde man die 
Teilchendurchmesser des Silbersols auf etwa 100 bis 150 A schätzen, was man wohl 
auch als plausibel ansehen darf. Infolge des Fehlens der Band-Band-Absorption bei 

Quecksilber kommt dort die Absorptioy 


a im kurzwelligen Teil des Sichtbaren be- 
= sonders gut zum Vorschein, obgleich 
Be \ die Absorption im Vergleich zum Silber 
\ , hier etwas kleiner ist. Auch sieht man 

deutlich, daß das Maximum kurzwelliger 

06 zz liegt als das Maximum beim Silber. Be- 

05 I \ achtet man noch, da8 auf der langwelligen 

04 Seite die Absorption von der lang- 

03 l welligen und „klassischen“ Absorption 
En H vg in Benzol des Elektronengases herrithren muB, so 
ist die Annahme, daß die langwellige 
: el etwa A hegt, mit den ave- 

| 4 | ne mannschen Ergebnissen durchaus ver- 
300400500 500 700 einbar. 4500 Ä wäre aber gerade die 


Wellenlänge, die sich aus dem von 
Silber- und Quecksilberkolloiden nach Pip pard gewonnenen Wert für die Ein- 
Messungen von R. Havemann. Abszisse dringtiefe nach unserer Theorie für die 
und Ordinate haben die gleiche Bedeutung Grenzwellenlänge der vom Elektronen- 
wie in Abb. 1 plasma absorbierten Strahlung ergeben 
‘ müßte. Wie Havemann weiter mitteilt, 
liegen die Plasmaabsorptionsstellen bei Cadmium und Zink noch etwas kurzwelliger 
als bei Quecksilber und zwar bei Zn: noch kurzwelliger als -bei Cd. Die Elek- 
tronendichten von Quecksilber, Cadmium und Zink wachsen in der Reihenfolge 
der Aufzählung der drei Metalle, wobei alle drei genannten Metalle im Gegensatz 
zu Ag als zweiwertig angesehen werden. Es könnte danaab den Anschein haben, 
als ob die Langmuirsche Formel zwar absolut genommen nicht die richtigen 
Werte für die Grenzwellenlängen der Plasmaabsorption liefert, daß sie aber 
wenigstens qualitativ die Reihenfolge für das Einsetzen der Plasmaabsorption bei 
den verschiedenen Metallen richtig vorauszusagen gestattet. 


Abb. 2. Darstellung der Absorption von 


Es mag nun noch kurz darauf hingewiesen werden, welche Rolle die longitudi- 
nalen Plasmaschwingungen bei der Absorption von Licht spielen kénnen. Bisher 
hatten wir ja nur die Absorptionsintensitat der transversalen Plasmaschwingungen 
berechnet. Man müßte also, um die Absorptionsintensität der longitudinalen 
_ Plasmaschwingungen auszurechnen, die Wechselwirkungsenergiedichte zwischen 
den Plasmaschwingungen longitudinalen Charakters und den einfallenden elektro- 
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magnetischen Wellen, die transversal sind, berechnen. Als Ausdruck der Wechsel- 
wirkung könnte man das Raumintegral über die Größe en, 3 € zugrunde legen, 
wobei en,3 das Dipolmoment der Plasmaschwingungen und € die elektrische 
Feldstärke der einfallenden Wellen bedeutet. Aus der Forderung der Longitudinali- 
tät der Plasmaschwingungen folgt jetzt, daß der Ausbreitungsvektor § der Plasma- 
schwingungen in der Richtung des Vektors n, 8 liegt, während f, der Ausbreitungs- 
vektor der elektromagnetischen Wellen, nach wie vor senkrecht zu € ist. Man 
sieht daraus schon, daß die Änderung des Impulses, die im wesentlichen durch die 
Vektordifferenz § — f gegeben ist, einen Vektor darstellt, der weder in die Rich- 
tung von $ noch in die Richtung von € zeigt. Dies hat zur Folge, daß in den Aus- 
sin — f),a,] 
sin —f),a,] ’ 


als Faktor eingeht. Man könnte 


druck der Übergangswahrscheinlichkeit nicht nur die Funktion 
sin [U n(8 — f), a] 
sin [2 (St — f), ao] 
daher schon bei Schichten endlicher Ausdehnung nicht mehr eine nennenswerte 
Absorption durch longitudinale Schwingungen erwarteny Höchstens im Falle 
von kolloidalen Lösungen, bei denen die Teilchengrößen kleiner als 20 Atomlagen 
im Durchmesser betragen, kénnte man eine Absorption durch longitudinale 
Schwingungen bemerken. Aber auch diese kénnte das Bild des Absorptionsvor- 
ganges nicht wesentlich ändern, da die Grenzen der Absorption die gleichen sind 
wie bei den transversalen Schwingungen. Es geniigt daher, die letzteren allein 
zu betrachten. 


«sondern auch z. B. die Funktion 


Berlin-Buch, Institut für Festkörperforschung der Deutschen Akademie der 
Wissenschaften und Berlin, Institut für theoretische Physik der Universität sowie 
Berlin, II. Physikalisches Institut der Universität. 


(Bei der Redaktion eingegangen am 22. Dezember 1948.) 
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Anwendung der Vierpoltheorie auf die Probleme der Ban 
Reflexionsminderung, Reflexionsverstärkung 
und der Interferenzfilter 


Von Kurt Schuster 
Be 
(Mit 2 Abbildungen) 


Inhaltsübersicht 


Die bei der Behandlung von Leitungsvorgängen der Hochfrequenztechnik und 
Akustik seit langem bewährte Vierpoltheorie kann auch auf die Lichtausbreitung 
durch eine Folge planparalleler Schichten Anwendung finden. Für senkrechten 
Lichteinfall werden die Grundgleichungen angegeben. Als Beispiel wird das 
Problem u. Entspiegelung durch eine Doppelschicht dargestellt. 


Durch Aufbringen geeigneter dünner Schichten läßt gich die störende Reflexion 
an Linsen weitgehend beseitigen!) und das Reflexionsvermögen von Metallspiegeln 
stark erhöhen). Durch Aneinanderreihen einer Anzahl dünner Schichten aus ver- 
schiedenen Stoffen können Spektralfilter von sehr geringer Durchlaßbreite kon- 
struiert werden®). Die erzielten Wirkungen beruhen auf der Lichtinterferenz an 
dünnen Schichten. Zur theoretischen Behandlung der genannten Probleme kann 
mit Vorteil das aus der Hochfrequenztechnik und Akustik bekannte Verfahren 
der Vierpoltheorie*) herangezogen werden. Besonders einfach liegen die Ver- 
hältnisse bei senkrechtem Lichteinfall. Im folgenden soll daher dieses Verfahren 
in seiner Nutzanwendung auf die Optik kurz beschrieben und seine Brauchbar- 
keit an einem speziellen Beispiel dargelegt werden. 


Die Wirkung jeder Schicht ist durch die drei Größen Brechungsindex n, Ab- 
sorptionskoeffizient x und Dicke d gegeben. Die Kreisfrequenz werde mit w, die 
Lichtgeschwindigkeit im Vakuum mit c bezeichnet; A und v sollen Wellenlänge und 
Lichtgeschwindigkeit des betreffenden Mediums bedeuten. Die reelle Wellenzahl 
ist dann gegeben durch 


9 
4) Vgl. A. Smakula, Glastechn. Ber. 19, 377 (1941). 
8) Vgl. R. Meßner, Optik 2, 228 (1947). \ : Ge 
3) Vgl. W. Geffeken, Angew. Ch. A 60, 1 (1948). Paar 


. 4) Für die Akustik, sind die entsprechenden Gedankengänge z.B. in dem Artikel 
des Verfassers ,,Die Messung mechanischer und akustischer Widerstände‘ (Ergebnisse 
der exakten Naturwissenschaften 21, 313 (1945) kurz dargestellt. Im übrigen sei auf die 
ienschlägige Lieratur der Schwachstromstechnik und der Hochfrequenztechnik verwiesen. 
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Die für die Lichtausbreitung maßgebenden Größen sind der komplexe Brechungs- 
index g und die komplexe Wellenzahl p. Für q und p gelten die Gleichungen: 


q=n-(1—jx) 
© 
p=k-(l—jx) 9. 


Die Lichtstrahlrichtung, die senkrecht zu den Flächen der Schichtenfolge ange- 
nommen wird, bilde die x-Achse eines kartesischen Koordinatensystems. Die 
Lichtwelle sei in Richtung der z-Achse linear polarisiert, eine Annahme, die bei 
isotropen Stoffen und senkrechtem Einfall keine Einschränkung der Allgemeinheit 
bedeutet. Dann wird in jeder Schicht der ebene elektromagnetische Wellenvor- 
gang durch die beiden Feldkomponenten E, und H, beschrieben. In dem die 
Schichtenfolge abschließenden Medium, das als unendlich ausgedehnt angenommen 
wird, besteht der Lichtvorgang aus einer fortschreitenden Welle. In den einzelnen 
dünnen Schichten und in dem der Schichtenfolge vorgelagerten Medium setzt 
sich der Vorgang im allgemeinen aus je einer hinlaufenden und einer zurück- 
kommenden Welle zusammen. Das Verhältnis Z,/H, in einer rein fortschreiten- 
den ebenen Welle ergibt sich aus der Maxwellschen Theorie bei Benutzung des 
egs-Systems (E, in elektrostatischen, H, in elektromagnetischen Einheiten) zu 


v 
aa Im’ bei Benutzung des internationalen Systems (E, in Volt/em, H, in 
Ampere/cm) zu my I = x) Ohm. Die Größe E,/H, in einer rein fortschreitenden 


ebenen Welle soll als der „optische Wellenwiderstand‘“ Z des betreffenden 
Mediums bezeichnet werden. Der optische Wellenwiderstand ist in nicht absor- 
bierenden Stoffen reell, in absorbierenden komplex. Wird das Verhältnis E,/H, 
innerhalb eines beliebigen ebenen Wellenvorgangs gebildet, der sich aus einem hin- 
laufenden und einem zurückkommenden Anteil zusammensetzt, so soll allgemein 
vom „optischen Widerstand‘ 3 gesprochen werden. 

In einer bestimmten Schicht 1 mit dem komplexen Brechungsindex g, und 
der komplexen Wellenzahl p, sei nun der Wellenvorgang unter Benutzung des 
egs-Systems in folgender Weise dargestellt: 

EYP = A, B,etipıa (1) 
Die Koordinate x, soll ihren Nullpunkt in der der Lichtquelle zugewandten Grenz- 
fläche der Schicht haben. Werden die Werte der Feldkomponenten in dieser 
Grenzfläche kurz mit E, und H,, die Werte der Feldkomponenten in der der Licht- 
quelle abgewandten Grenzfläche der Schicht mit Z, und H, bezeichnet, so wird 


E, = A, +B, E, A, -e—imd + 
A, = q, — H,=q -e— — B, 
Dabei bedeutet d, die Dicke der Schicht. Durch Eliminierung der komplexen 
Amplituden %, und ®; erhält man einen unmittelbaren Zusammenhang zwischen 


den Eingangs-Feldgrößen und den Ausgangs-Feldgrößen der Schicht in Form der 
sogenannten. Vierpolgleichungen : 


(2) 


E, = cos p,d,-E,+2 snp,d,:-H, 
sen. H, = jq, sin p, d, - FE, + cosp,d, - Hy. 
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Durch Division der beiden Gleichungen ergibt sich ein Zusammenhang zwi- 
schen den Werten 3, und 3, des optischen Widerstandes in der Eingangs- und 
Ausgangsgrenzfläche der Schicht: 


cos p, d, + -sin p, d, 
= j qı 32 sin p, d, + cos dy 
SE Es ist zu beachten, daß sich die Indices bei p,, 9,, d, auf die Schicht selbst, 
bei den optischen Widerständen 3, und 3, auf die Grenzflächen der Schicht be- 

‘ ziehen. Die Vereinfachung der rechnerischen Handhabung der Probleme, die 
die Vierpoltheorie mit sich bringt, beruht letzten Endes auf der neuen Formu- 

lierung der Bedingung der Stetigkeit beim Durchgang durch eine Grenzfläche. 

Ss Die beiden Stetigkeitsbedingungen für die Feldkomponenten E, und H,, wie 
2 sie bei der üblichen Behandlung des Problems angesetzt werden, können dureh 
die Vierpoltheorie zu der einen Stetigkeitsbedingung für den optischen Widerstand 
zusammengefaßt werden. Der optische Widerstand in einer Grenzfläche ist je- 
weils durch die dahinter angeordneten Teile der Schichtenfolge gegeben. Gl. (4) 


3 ermöglicht es, das System der Schichtenfolge bis zur Eingangsgrenzfläche der 
£ ersten Schicht rückwärts durchzurechnen; sie ist dementsprechend als die Glei- 


chung der Widerstandstransformation zu bezeichnen. Für die optischen Auf- 
gabestellungen ist man nicht. unmittelbar am optischen Eingangswiderstand, 
sondern am Reflexionsvermögen interessiert. Das komplexe Amplituden-Re- 
flexionsvermögen St = B,/A, an der ersten Grenzfläche ist in bekannter Weise 


gegeben durch 
R— 3:—% (5) 
; B wobei Z, den optischen Wellenwiderstand des der Schichtenfolge vorgelagerten 
h Mediums und 3, den optischen Widerstand in der ersten Grenzfläche bedeutet. 
3 Das Energie-Reflexionsvermögen R ergibt sich daraus in einfacher Weise zu 3 
R=|KP=—R- R*, 3 
2 wobei R* den zu R konjugiert komplexen Wert bezeichnet. 
he Die Anwendung der Gl. (4) und (5) auf die Probleme der Reflexionsvermin- 
ie. derung, Reflexionserhöhung und der Interferenzfilter ermöglicht die Lösung der 
4 gestellten Aufgaben in besonders übersichtlicher Weise. Als Beispiel soll im folgen- 
Er. den die Entspiegelung durch eine Doppelschicht dargelegt werden. In der vor- 
# handenen Literatur®) wird die Dicke jeder Einzelchicht zu einer Viertelwellen- 
3 länge angenommen, und es werden die Wertepaare n,, n, berechnet, die zu einer 
B- A vollkommenen Entspiegelung führen. Die folgende Berechnung nach der Vierpol- 
7 theorie wird dariiber hinaus zeigen, daB in gewissen Bereichen ein beliebiges 
Wertepaar n,, zur Realisierung der vollkommenen Entspiegelung herangezogen 
We werden kann, wobei dann die zu wählenden Schichtdicken im allgemeinen von der 


Viertelwellenlänge abweichen. 

In Abb. 1 ist die zu untersuchende Anordnung schematisch dargestellt. Die 
beiden Schichten besitzen die Brechungsindices n, und n, und die Schichtdicken 
d, und d,. Der Brechungsindex des vorgelagerten Mediums hat den Wert n, 
der des abschließenden Mediums den Wert n,. Alle Stoffe sind nicht absorbierend 


5) R. L. Mooney, J. Opt. Soc. Am. 35, 574—583 (1945); Ref. Optik 2, 164 AT); 
L. Hiesinger, Optik 3, 485 (1948). 
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(x = 0,q=n,p=k). Der Wellenwiderstand des vorgelagerten Mediums beträgt 
also 1/n,, der des abschließenden Mediums 1/n,. Entsprechend der Begriffsbildung 


_ der Vierpoltheorie besteht die Aufgabe darin, für eine bestimmte Frequenz die 
/ Wellenwiderstände der beiden Medien durch einen Transformator einander ,,anzu- 


passen“, um dadurch die Reflexion zu beseitigen. Die optischen Widerstände 
an den drei Grenzflächen werden mit 3,, 3», 


3; bezeichnet. Die stufenweise erfolgende 
“Berechnung der Widerstände ergibt n\m|m|m 
3 
“ ing‘ 83° sin k,d, + cos kody (6) 
2° 3s 2 
cosk,d, + —-sink,d, Abb. 1. Schema der Doppelschicht 
jn, sin id, cosk,d, 


Im Falle der Reflexionsfreiheit i = 0 muß nach (5) 3, = Z, = 1/n, gemacht 
werden. Durch Einsetzen dieses Wertes in die letzte Gleichung nimmt die Be- « 


dingung der Reflexionsfreiheit die Form an: 


Durch Gleichsetzen mit der,zweiten Gleichung von (6) erhält man nach einigen 
einfachen 


tg hy dy te ky dy = 
” Ny 


Hieraus ergeben sich sofort die Gleichungen zur Bestimmung der beiden 2 


Schichtdicken: N 
ag 


é I ny No Na No N, 


Ne Na 
Ns ny (2-25). (= —*) 
Die Vorzeichen der Wurzeln sind im Einklang mit den Gleichungen (7) zu 
wihlen. Die Diskussion der Gleichungen (8) soll an Hand der Abb. 2 erfolgen. 


Hier sind die Brechungsindices n, und n, als kartesische Koordinaten aufgetragen. 
Die Brechungsindices n, und n, sind als fest vorgegebene Größen angenommen. Die 


. 
= 
2 
- und a 
‘ 1 
4 
_| 


Np 14 Ny 
Die gleichen Klammerausdrücke bilden auch Teilfaktoren der Radikanden 
in den Gleichungen (8). In den von den drei Geraden begrenzten schraffierten Ge- 
bieten ist das Vorzeichen der beiden Radikanden (8) positiv und damit die Auf- 
gabe lösbar. Die schräge Gerade stellt 
_ den von Mooney und Hiesinger be- 
handelten Spezialfall dar. Die vor- 
. . a liegende Rechnung zeigt, daß zur Er- 


ey reichung vollkommener Entspiegelung 


die Wertepaare n,, n, nicht auf der 
schrägen Geraden zu liegen brauchen, 
sondern irgendwo in’ den schraffierten 
Gebieten liegen können. Nur sind dann 
die zugehörigen Schichtdicken nicht 
My =V ehr gleich einer Viertelwellenlänge, 


m 
N NN sondern entsprechend den Gleichungen (8) 
um zu wählen. 
N 75 Fr Bei der Lösung komplizierterer Auf- 


gaben (Reflexionsverstärkung von Metall- 
spiegeln durch mehrere Schichten, Inter- 
_ Abb. 2. Bereiche, in denen die Aufgabe ferenzfilter) wird man die numerische 
der Entspiegelung lösbar ist Berechnung vorziehen. Man kann dann 
nach dem Schema der Gleichungen (6) 

für eine Anzahl herausgegriffener Wellenlängen die optischen Widerstände stufen- 
weise rückwärts numerisch berechnen und zuletzt mit Hilfe von Gleichung (5) 
den Übergang zum Reflexionsvermögen vollziehem. Allgemein ist zu erwarten, 
5 ge. daß die Anwendung der Vierpoltheorie auf optische Probleme ebenso wie in der 
ie -Hochfrequenztechnik und Akustik neue Gesichtspunkte für die Konstruktion 
Wellenfiltern liefern wird. 
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Jena, Carl Zeiss-Werke. 


7 
(Bei der Redaktion eingegangen am 28, Dezember 1948. 
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Ne no) n 
Ns 
tle 
die 
i und 
Sy: 
m 
tri 
kla 
D 
| 
| 
3 
. 


canden 
en Ge- 
e Auf- 
stellt 
be- 
‘vor- 
ir Er- 
gelung 
if der 
uchen, 
ierten 
dann 
nicht 
länge, 
zen (8) 


r Auf- 
[etall- 
Inter- 
rische 
dann 
n (6) 
tufen- 
ig (5) 
arten, 
n der 
ıktion 


I. Termaufspaltung und 
Ann. Physik 


Berichtigung zur Arbeit v von K.H. Hellwege: 
Elektronenterme und Strahlung von Atomen in Kristallen 


elektrische Dipolstrahlung 
(6) 4, 95 (1948) 


In die genannte Arbeit haben sich leider einige Schreibfehler eingeschlichen, 


die wie folgt zu korrigieren sind. 


In Tabelle'2 muß es 8. 110/111 und §. 113/114 fiir die Symmetrieklassen Dog 


und Ds, richtig so heißen: 


1 2 3 4 5 
Sym- Termfamilien Voll- 
me- | nach Anzahl der Terme in jeder Familie — 
t G bei be J . 
D. 
‘ 1 ( J+2] , [J—2 
(14 
1 lk 2 jk) 
)(-1* 
Zir J+2 , [J—2 
> 
=) 
2 ix) 
4 
1 21 [J—2 
1 
z)) 
1 » (le - 5k) 
J+1 


= 
. “ i 
ins 
a 
. 
. 
J 
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2 | 3 4 5 
| Termfamilien ‘ Voll- 
nachd. oii Anzahl der Terme in jeder Familie _| stand. 

i b Aufsp. 
| penth| Kristallquantenzahlen bei gegebenem J rm, 
D 
J+3] [J—3 
ix) 
+ (14(—1)% )(-1)% 
14(—1)* ) 
al 
( (tk - 5k) 
= (—1)* 
 [J+3] 
1 
(1-(-1)* ‘)(1+2[2]) 
Lik 2' - ix) 
J+3 
1+(—1)* 
a1 
1—(—1)* )(1+2[2]) 
2 (le - ik) 
i—(—3)5 
J+1 J—1] 
(1+| sur ) 
zı J+2]  [J—2]\ 
1-93 
L | 
Zh J+21 fy—2 
J+1] 
1-9 c, 
Ferner heißt auf 8. 113 die letzte Zeile von Dg richtig: 
9 
| » 


Pei . Auf S. 99 ist hinter Gl. (14) AN Satz „wobei m eine beliebige ganze Zahl 


ist“ zu streichen. 
Auf S. 104 in der Zeile vor Gl. (43): u = ®/, statt u = !/, und in der Zeile 
nach Gl. (43): Z’ statt Z,. 2 
Auf S. 105 letzte Zeile: Z’ statt Z,. 
Göttingen, II. Physikalisches Institut. 


(Bei der Redaktion eingegangen am 20. Januar 1949.) 
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